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Automates et systèmes à événements discrets 


Avertissement 


Ces notes ont été élaborées depuis 92 dans le cadre d’un cours donné selon les années en 
plusieurs endroits (troisième année à l’ENSTA, option automatique à l ’École des Mines, DEA 
Automatique et Traitement du Signal d'Orsay, DEA Automatique-productique, ENSIEG). 

En 21 heures, avec un public n'ayant jamais entendu parler de systèmes à événements 
discrets, les quatre premiers chapitres sont en général traités. 

Hormis la théorie de la commande supervisée, initiée par Ramadge et Wonham, ce cours 
s'inspire largement du travail commun au sein du petit groupe de travail sur l’algèbre max- 
plus, à l’INRIA Rocquencourt, auquel ont participé ou participent toujours, en dehors de 
l’auteur, Guy Cohen, Jean-Pierre Quadrat, Michel Viot, Pierre Moller, Ramine Nikoukhah, 
Marianne Akian. D'autres chercheurs et d’autres écoles ont contribué fortement par ailleurs 
au développement de cette théorie, voir la petite bibliographie raisonnée à la fin du premier 
chapitre. 

Depuis sa première rédaction, ce cours a subi des améliorations essentiellement 
cosmétiques (le plan étant inchangé). Notre intention initiale était de présenter en détail 
la théorie spectrale max-plus, qui est à notre avis le cœur du sujet. On obtient comme corol- 
laire du théorème spectral des résultats importants, par exemple le théorème de périodicité 
énoncé par Chretienne pour les graphes d'événements, ou des théorèmes asymptotiques en 
commande optimale déterministe (à espace d’état et de commande finis) qui n'apparaissent 
pas en général avec ce degré de précision dans la litérature. 

La théorie des séries rationnelles à coefficients sur le semi-anneau max-plus, qui offre 
un point de vue alternatif (de puissance équivalente) à la théorie spectrale, est seulement 
esquissée. La théorie de la résiduation, qui fournit un cadre algébrique au calcul des dates 
au plus tard, n’est pas abordée. En attendant la nouvelle version de ce cours, nous renvoyons 
le lecteur à “Synchronisation et Linéarité”, op. cit. 

Enfin, avec un peu de recul, ils nous semble que le progrès important de ces dernière 
années, et qui n'apparaît qu'en filigrane dans la version présente du cours écrit, et 
qui apparaît espère-t-on mieux dans le cours oral, est venu de l’augmentation du pou- 
voir de modélisation, grâce d’une part aux généralisation non-commutatives des graphes 
d'événements (automates max-plus, empilements de pièces), qui permettent de modéliser des 
phénomènes de concurrence (e.g. partage de ressources), et d’autres part aux systèmes dyna- 
miques monotones homogènes, et en particulier aux “fonctions min-max” et “réseaux de Pe- 
tri fluides”, qui constituent des classes maintenant assez mieux comprises et algébriquement 
analysables de systèmes. 


Chapitre 1 


Introduction 


1.1 Petite revue des Systèmes à Événements Discrets 


1.1.1 Problèmes typiques et outils 


Ce cours a pour but de présenter et d’illustrer les techniques d’analyse des systèmes dits 
“à événements discrets”, pour lesquels des “événements” (arrivée d’un client, achèvement 
d’une tâche, envoi d’un signal) donnent lieu à des phénomènes de synchronisation, saturation 
ou concurrence. On trouvera de tels phénomènes dans les systèmes de production (ateliers 
flexibles) les systèmes multiprocesseur, les réseaux informatiques, réseaux de transport ..…. 
Les systèmes à événements discrets peuvent se définir négativement par rapport aux systèmes 
considérés en automatique classique. On peut dire que l’automatique classique étudie les 
systèmes de la forme x’ — f(x,u) où f est une fonction régulière (C®, voire linéaire) de 
l’état x et d’un contrôle w, ainsi que leurs analogues discrets. Il est clair que les phénomènes 
de saturation et de synchronisation propres aux systèmes à événements discrets se laissent 
difficilement modéliser dans ce cadre (concédons que certains systèmes à événements dis- 
crets sont réguliers par morceaux, mais toute la problématique réside précisément en l’étude 
des ‘“‘changements de morceaux”). Les problèmes qui nous préoccupent sont typiquement 


— la spécification de systèmes à événements discrets (description d’un comportement 
logique ou temporel admissible), 


— la conception (le ‘“‘design”) de ces systèmes 


— la détermination qualitative et quantitative du comportement des systèmes à 
événements discret, soit plus particulièrement 


— l’étude du comportement asymptotique et la mise en évidence d’un régime sta- 
tionnaire ou périodique simple éventuel, 


— les questions d’évaluation de performance (typiquement, calcul du taux de pro- 
duction) 


— les questions de stabilité : les stocks de pièces ou le nombre de données en 
mémoire restent-ils finis, comment le système réagit-il à d'éventuelles pertur- 
bations ? 


— certains problèmes d’optimisation : par exemple, minimiser le nombre de ressources 
(machines, palettes, processeurs) nécessaires pour réaliser un taux de production 
donné. 
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En l'absence de théorie standard des systèmes à événements discrets, de multiples techniques 
ont pu être appliquées, certaines classiques, d’autres plus exotiques. Nous mentionnons en 
particulier : 


— Les approches “expérimentales” qui simulent ces systèmes sur ordinateur. Nous n’in- 
sisterons pas ici sur l’importance pratique de la simulation, qui constitue souvent le 
point de passage obligé avant la mise en place de tout modèle. Il y a tout un cor- 
pus théorique autour de la simulation, dont l’objet est de fabriquer de bons esti- 
mateurs (économes en temps de calcul), surtout pour des quantités intéressantes et 
difficiles (coûteuses) à calculer avec une précision raisonnable (gradients, fonctions 
d’événements rares). Nous mentionnerons particulièrement les approches de type ‘“ana- 
lyse de perturbation” [66, 39] qui visent à déterminer l’effet d’une perturbation sur le 
système sans avoir à le re-simuler, ce qui permet en particulier de mettre en œuvre des 
algorithmes d’optimisation de type gradient stochastique. 


— Les modélisations probabilistes classiques (par exemple par des réseaux de Files d’at- 
tente ou des chaînes de Markov). L’on peut mettre en évidence des bonnes classes de 
systèmes pour lesquels on a certaines formules analytiques (par exemple, les réseaux 
de Jackson) ou pour des systèmes dont la description totale est trop complexe, on 
cherchera des résultats simples asymptotiques, un peu dans l’esprit de la mécanique 
statistique qui met en évidence des lois limites pour des systèmes de grande taille. 


— Réseaux de Petri (temporisés) : il s’agit d’un langage graphique propre à la description 
et l’analyse des systèmes à événements discrets. 


— Systèmes de transitions finis (automates finis). Cette approche traite principalement 
du comportement logique de ces systèmes : on cherche à vérifier certaines propriétés 
comme l’absence de blocage ou l’équité entre les différentes tâches à servir. L’on peut 
exprimer les propriétés requises pour les trajectoires du système dans certaines lo- 
giques dites temporelles! Ou bien, comme dans la théorie de Ramadge & Wonham, 
on exprimera les spécifications du système par l’appartenance à un langage légal (re- 
connu par un automate fini dans les cas les plus simples). Il y a ensuite toute une 
algorithmique qui permet de vérifier des sous classes de propriétés, que l’on ramènera 
par exemple à des problèmes d’accessibilité dans des graphes (en général gros). Pour 
résister un peu à la? “Np-complète dérision”, on ne représente pas toujours l’automate 
en extension (par sa table), mais on code parfois sa fonction de transition par une for- 
mule Booleenne, qui peut se représenter de manière expérimentalement compacte sous 
forme de BDD (diagramme de décision binaire). 


— Systèmes dynamiques sur des corps finis. On code les événements et l’état du système 
par des vecteurs à coefficients dans Z/pZ (corps des entiers modulo p pour p premier) 
et l’on décrit la dynamique du système par des équations récurrentes polynômiales. Le 
pouvoir d’expression est identique aux systèmes de transition finis, mais la nouveauté 
consiste a voir les trajectoires comme des solutions d’un système dynamique poly- 
nomial, ce qui permet de ramener des problèmes dynamiques (e.g. calculs de points 
fixes, calcul d’ensembles invariants modulo l’action des commandes) à des problèmes 
de variétés algébriques ou d’idéaux, qui sont bien répertoriés mais assez coûteux en 
temps de calcul. 


ÎLe terme “temporel” dit simplement que l’on parle des trajectoires (il s’agit d’un temps logique, à ne pas 
confondre avec le temps physique). Pour introduction détaillée, voir A. Arnold [2]. 

2Nous avons emprunté l'expression au poème de A. Lentin qui fait prélude à [51]. Le terme “curse of dimen- 
sionality” est sans doute plus classique. 
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— Systèmes dynamiques sur des algèbres exotiques (de type max, + ou min, +). 
Dans ce cours, nous utiliserons essentiellement trois de ces outils : 


— les réseaux de Petri (temporisés), qui sont à la base de la spécification, représentation 
et simulation de ces systèmes. 


— l’approche “systèmes dynamiques linéaires”, qui est analogue à la théorie usuelle clas- 
sique des systèmes linéaires, mais après un changement de structure algébrique. Ce 
cadre ne s’applique qu’à une sous classe de systèmes à événements discrets. 


— l’approche à la Ramadge & Wonham via la théorie des automates. Celle-ci a un pou- 
voir de modélisation plus grand —on peut traiter certains phénomènes de concurrence. 
On privilégie les problèmes de nature logique” (par exemple absence de blocage). Ici 
encore, le point de vue est celui de l’automaticien qui voit le système à événement dis- 
cret comme un système dynamique et lui étend les notions familières (commandabilité, 
observabilité) à l’aide des outils standards de théorie des langages formels. 


La littérature sur les réseaux de Petri est fort riche, et comprend en particulier de nombreux 
résultats analytiques. Ici, nous n’utilisons les réseaux de Petri (en fait, principalement, les 
graphes d'événements) qu’à des fins de modélisation. On décrit les systèmes en termes de 
réseau de Petri, puis on les analyse au vu des équations linéaires associées. C’est dire que 
ce cours se limite aux techniques élémentaires de nature algébrique utiles à l’étude des 
systèmes à événements discrets. Nous pensons qu’il s’agit là de la meilleure introduction 
à ces systèmes, dans la mesure où ces approches algébriques (via la théorie des automates et 
les systèmes (max, +)-linéaires) produisent les résultats les plus simples relatifs à des sous- 
classes de systèmes à événements discrets. Il faut bien voir ici que les systèmes à événements 
discrets forment une classe trop large pour que l’on puisse prétendre donner des résultats 
précis généraux. Il ne s’agit pas ici de tout résoudre, mais de distinguer ce qui se résoud et ce 
qui ne se résoud pas par la théorie et de comprendre la pertinence physique des sous-classes 
mises en évidence (essentiellement les graphes d'événements et les automates). En outre, 
ces sous-classes ont valeur de modèles pour passer aux systèmes à événements discrets plus 
généraux, pour lesquels on aura des difficultés supplémentaires. Par exemple, la théorie des 
graphes d'événements stochastiques [4] se comprend mieux si l’on sait déjà ce qui se passe 
dans le cas déterministe. De même, la théorie des réseaux de Petri (non temporisés) est moins 
aisée que la théorie des automates de base dans la mesure où les langages qui décrivent les 
comportement des réseaux de Petri sont plus généraux que les langages dits réguliers (ou 
rationnels) qui sont reconnus par des automates finis [54, 61]. 


1.1.2 Quelques exemples de Systèmes à Événements Discrets 


EXEMPLE 1.1.2.1 (Assemblage). On considère le processus suivant. Dans un atelier, on as- 
semble une pièce de type A et une pièce de type B ce qui prend une durée +. Si uA(t) (resp. 
ug(t)) dénote le nombre de pièces de type À (resp. B) arrivées jusqu’à l’instant fr, et y(f) 
représente le nombre de pièces produites à l’instant f, on a l’équation : 


y) = min(ux( —T),ug(t —T)) . (1.1) 


3Il est possible d'intégrer le temps dans ce formalisme au prix d’une complexité plus grande. On peut comme 
Brandin et Wonham, rajouter une nouvelle lettre r dont l’occurrence représente un top d’horloge [15]. On peut 
aussi à la suite d’Alur, Courcoubetis et Dill [1] introduire des automates temporisés, c’est-à-dire des automates 
dont les transitions sont conditionnées par les valeurs de compteurs. Ce dernier formalisme est plus puissant, en 
ce qu’il autorise la modélisation d’un temps continu (réel) et non discret. On peut aussi utiliser des automates à 
multiplicités sur le semi-anneau (max,+) [37]. 
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Les applications u1,u8, y qui mesurent des quantités méritent le nom de compteurs. Duale- 
ment, on peut introduire les fonctions dateur : u’,(n) = date d’arrivée de la n-ième pièce de 
type À, u',(n) = ..., y'(n) = date de production de la n ième pièce. On a alors 


y'(n) = T + max(u,(n),u%(n)) . (1.2) 


EXEMPLE 1.1.2.2 (Contrainte de débit). Considérons une unique machine pouvant traiter 
une seule pièce à la fois pendant t unités de temps. Soit u(t) le nombre de pièces brutes 
arrivées à l’instant r et y(f) le nombre de pièces achevées au même instant. On peut écrire : 


y) = min(u(t —t),y({(—-Tt)+1). (1.3) 


Dualement, en notant w’(n) la date d’arrivée de la n ième pièce brute et y’(n) sa date 
d’achèvement, on a 


y'(n) = max(u'(n)+t,r +y(n—1)). (1.4) 


EXEMPLE 1.1.2.3 (Cantonnement). Dans une version simplifiée, une portion de voie ferrée 
fonctionne comme suit : un train ne peut rentrer que si le train précédent est déjà sorti (sinon 
le feu est rouge). Les trains sont supposés traverser un canton en un temps constant tr. Soient 
u(t) le nombre de trains arrivés à l’entrée du canton à l’instant fr, x, (1) le nombre de trains 
rentrés, et x2(f) le nombre de trains sortis. On a 


X1() = min(x2(f) + 1, u(t)) 
H()=x(—-T). 


En éliminant x, on obtient x2(f) = min(x2(f —t)+1,u(t —Tt)), soit une contrainte de débit 
identique à (1.3) ci-dessus, ce qu’on aurait pu voir directement comme suit. Du point de vue 
des comportements temporels, il y a isomorphisme entre une machine et un canton : le train 
joue le role d’une pièce, et la contrainte d’exclusion mutuelle des trains sur un canton revient 
à dire que la machine ne peut traiter qu’une pièce à la fois. 


EXEMPLE 1.1.2.4 (Communications Inter-Processeur). Un processeur À adresse des mes- 
sages vers un processeur B de la manière suivante. Les messages prêts à l’envoi sont placés 
en attente dans un tampon de capacité suffisante. À envoie le premier message à B ce qui 
prend un temps t (que l’on pourra supposer constant si la taille des messages est fixe). À 
attend d’avoir reçu un accusé de réception de B (qui met t’ unités de temps pour parvenir de 
B à À) pour en envoyer un nouveau. Soit #(f) le nombre de messages arrivés dans le tampon 
jusqu’à l’instant f, et soit x1(r) le nombre de messages envoyés par À jusqu’à l’instant f. On 
a 


XA() = max(xa(t—T—-T)+1lu(t)), 


où le lecteur, sans doute ecœuré, reconnait encore une variante de l’équation (1.3). Le même 
système dynamique apparait ainsi déguisé dans des applications différentes, d’où l’intérêt de 
développer une théorie des sytèmes à événements discrets qui parle de manière unifiée de ces 
systèmes (faisons un peu de prosélytisme). 


EXEMPLE 1.1.2.5 (Le professeur avisé). Un enseignant arrive à l’heure w’ pour faire son 
cours d’1h. Il a n élèves inscrits. L'élève à arrive à l’heure x’. Calculons l’heure y’ de fin 
de cours pour quelques politiques naturelles. Si notre enseignant attend tous les élèves pour 
commencer, le cours finira à 
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y'=1+max(u’, max x) , 
EEE 
où l’on reconnait une équation avec les lois max assez similaire aux équations écrites plus 
haut pour les dateurs, par ex. (1.2). Les choses étant ce qu’elles sont, une telle politique n’est 
pas toujours raisonnable, et notre enseignant peut choisir de ne jamais attendre plus de 1/4 
d’heure, quitte à faire cours devant une classe vide. Le cours finira alors à : 


y = 1 + min(u’ + 1/4, max(u’, ns KE (1.5) 


Une politique plus avisée peut être de se fixer un quorum g < n d’élèves, et de commencer 
dès que les g premièrs élèves sont présents. On trouve alors 


y'=1 + max(u, ( min max x!) ; (1.6) 
IC{1,..,n} iel 
card 1=q 


Ce que l’on pourra retenir de cet exemple, c’est que les dates des événements, dans des 
systèmes à synchronisation limitée (quand on attend, mais pas trop) peuvent “souvent” 
s’écrire en hybridant avec la loi min les équations donnant les dates des systèmes purement 
synchronisés, lesquelles font seulement intervenir les lois max, et +. 


EXEMPLE 1.1.2.6 (Magasin de taille finie). On considère un stock pouvant contenir au plus 
3 pièces. On a les deux événements 


a une pièce entre dans le stock 
b une pièce est déstockée. 


Ce système est représenté par le graphe de la Figure 1.1. Il s’agit d’un automate (détermi- 
niste). Il y a toute une théorie des automates (dont on introduit les éléments au chapitre 5). 
Ici, il nous suffit de savoir que de manière élémentaire, un automate est un graphe dont les 
sommets représentent les éfats possibles du systèmes et les arcs les changements d’états sous 
l'effet des événements (les lettres). L'état initial est distingué par une flèche rentrante. Ici, le 
sommet 0 représente l’état “0 pièce en stock” (initial), le sommet 1 une pièce en stock, etc. Par 
exemple, de l’état 1, les deux événements a, b sont admissibles (et conduisent respectivement 
à l’état 2 ou à l’état O0, comme les flèches l’indiquent). On cherche l’ensemble des séquences 
d'événements a et b admissibles (1.e. compatibles avec la capacité du magasin). Pour cela, à 
chaque chemin du graphe, on associe le mot lu de droite à gauche formé des lettres valuant 


les arcs. Aïnsi, au chemin 0 ES Lt 1 correspond le mot baa. Par contre, le mot 
ab ne correspond pas à un comportement possible de l’automate (l’événement “b” arrivant à 
l’état O fait vendre une pièce que l’on n’a pas en stock). Plus généralement, il est clair que 
l’ensemble des séquences d'événements admissibles correspond à l’ensemble des chemins du 
graphe partant de l’état initial 0. 


a a a 


RS CT a «7e 
0 HN 2 RE : uk) 
b D ou 


FIG. 1.1: Magasin de trois unités 
De manière plus précise, on spécifie l’automate par un alphabet Y (en l’occurrence Y = 


{a, b}), un ensemble d’états Q (ici O — {1,2,3}), un état initial go (ici go — 0) et une 
fonction de transition partielle (i.e. non nécessairement partout définie) 
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6:2xQ—0Q 


q' = ô(u,q) (1.7) 


est le nouvel état si l’on lit la lettre w à partir de l’état qg (ex : 6(a, 2) — 3, 6(b,2) = 1). On 
peut donc voir les automates comme une version finie des systèmes en temps discret familiers 
à l’automaticien, qui sont de la forme 


Xn+1 — (Un: Xn), Xn € R?, Un € R‘ 


clairement analogue à (1.7). 


EXEMPLE 1.1.2.7 (Chat et souris). On considère un prédateur et une proie se déplaçant (de 
manière imprévisible et asynchrone*) dans un ensemble de pièces communiquant par cer- 
taines trappes. Chaque trappe est réservé pour l’usage exclusif d’un seul animal, ï.e. seul le 
chat peut emprunter les trappes c1, c2 et seule la souris peut emprunter les trappes m1, m2. 
Notons “c,” l’événement “le chat passe par la trappe c2”, et ainsi de suite pour les autres 


P P: 
15 ARR 
KE 2 
Lei 
CI m] A+ 
m2 
P; 


FIG. 1.2: Chat et souris dans un labyrinthe et automate correspondant 


trappes. On se convainc que l’ensemble des comportements possibles du chat et de la souris 
est représenté par l’automate à droite de la Figure 1.2. Les sommets du graphe sont repérés 
par des couples “(i, j)” qui signifient “chat dans la pièce à et souris dans la pièce j”. Les 
arcs représentent les changements d’état relatifs aux mouvements d’un seul des deux ani- 
maux, notés par des lettres c; ou m; (les mouvements sont supposés immédiats, et effectuent 
à des instants arbitraires, on exclut en outre le cas où le chat et la souris changent de pièce 
exactement au même instant). Les états diagonaux de l’automate (1.e. les étas (5, i), pour 
i = 1,...,3, où le chat peut manger la souris) peuvent être qualifiés d’indésirables. Le but 
de la théorie exposée au chapitre 5 est de synthétiser des contrôles (par exemple de bloquer 
certaines trappes en fonction de la position du chat et de la souris) de manière à garantir un 
bon comportement du système, particulièrement en présence de contraintes (impossibilité de 
verouiller certaines trappes), voire même en observation partielle (lorsque l’on ne peut voir 
si certaines pièces sont occuppées). 


EXEMPLE 1.1.2.8 (Jeu de Tetris et systèmes à ressources partagées). Un “jeu de Tetris” ou 
système à ressources partagées peut être formalisé de la manière suivante. On se donne 


Le terme asynchrone signifie que le prédateur et la proie peuvent effectuer leurs mouvements à des instants 
physiques indépendants. On parle au contraire de système synchrone (par exemple pour un microprocesseur) 
lorsque toutes les actions ont lieu à des temps physiques donnés par les tops d’une horloge commune. 
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un ensemble de positions (en abscisse) ou de ressources R (e.g. R = {1,... ,n}). Nous 
appellerons pièce (ou tâche) a un bloc rigide éventuellement non connecté, représenté 
géométriquement par un ensemble de positions (ressources) occupées R(a) € R, un contour 
bas (dates de début d’exécution sur les ressources) £(a) : R(a) — R, un contour haut (dates 
de libération des ressources) h(a) : R(a) — KR, avec bien-sûr : Va € R(a), h(a) > £(a). Une 
pièce sera représenté par la région du plan R XR: P, = {(r,y) e R(a) XR | £(a), < y < 
h(a),}, ce que l’on interprétera comme suit : la tâche a requiert l’ensemble des ressources 
(machines, processeurs, main d'œuvre) R(a), et chaque ressource r € R(a) est employée 
de l’instant £(a), à l’instant h(a),. Une pièce P, peut être translatée verticalement d’une 
constante À, ce qui donne la nouvelle région £’(a) = À + £(a), h'(a) = À + h(a). Physique- 
ment, cela correspond à exécuter la tâche À unités de temps plus tard (ou —À unités de temps 
plus tôt si À < 0). Un sol ou une condition initiale est un vecteur ligne g € R? : la ressource 
r devient disponible initialement à l’instant g,. Si nous laissons tomber les k pièces a1 ...a, 
dans cette ordre, sur le sol g (on laisse tomber les pièces selon les lois de la gravité, en interdi- 
sant les translations horizontales et les rotations, exactement comme dans le vrai jeu de Tetris, 
voir la figure 1.3), nous obtenons un fas de pièces. Le contour haut x (w) du tas w — a; ...ax 
est le vecteur ligne dans R®, dont la composante r est égale à la hauteur de crête du tas au 
droit de la ressource r, 1.e. à la date de libération de la ressource r à l’issue de l’exécution de 
la suite de tâches w. La hauteur du tas est par définition y(w) — max,er x(w),. Physique- 
ment, y(w) donne ce que la littérature de l’ordonnancement appelle le makespan, c’est-à-dire 
le temps d’exécution de la suite de tâches, ou de manière équivalente, la date de libération 
de la dernière ressource occupée. On établit aisément la récurrence suivante qui détermine 
complètement la dynamique au plus tôt du système (_ dénote le tas vide) : 


X( = &, Vr e R(a), x(wa), — ne re — {(a)s +x(w)s] , 
Vr g R(a), x(wa), — x(w),. (1.8) 


CCE 
+ no 


R(c) = {2,4}, €(c) =[:,0,.,01, Ac) = [:,2,:,2] 


= 


Ê R(b) = {1,2}, €(b) = [0,0,.,-], h(b) = [2,2,:,.] 
R(a) = {1,2,3}, £(a) = [0,0,0,:.], h(a) =[1,1,3,:] 


FIG. 1.3: Un jeu de Tetris à trois types de pièces et quatre ressources 
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Arrivé à ce point, le lecteur a mérité une petite devinette : quel est le point commun entre 
les systèmes traités ci-dessus ? Cette section donne une première réponse : si l’on exclut 
les politiques trop avisées du professeur (Eq. (1.5) et (1.6)), nous allons voir que tous ces 
systèmes se représentent par des équations linéaires dans certaines structures algébriques 
exotiques. Nous verrons pour conclure que même l’exemple du professeur avisé, qui n’est 
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pas linaire, partage avec le cas linéaire certaines propriétés fondamentales (monotonie, et 
homogénéité). 


1.2.1 Systèmes (min,+)-linéaires 


Considérons l’ensemble RU{+c0}, muni du min (noté additivement, par exemple 263 — 
2 — min(2, 3)) et de l’addition usuelle notée multiplicativement (2 @ 3 — 5,0 @ 1 — 1). Ces 
notations sont justifiées par le fait que les propriétés usuelles sont vérifiées : 


(aBb)Bc—aæ(beæc) = min(min(a, b),c) = min(a, b, c) 
ab=b@a = min(a, b) 

a@b)@c—=a@(b@c)=a+b+c 
(aBb)@c—=a@cBb@c=min(a,b)+c=min(a+c,b+c). 


(1.9) 


Nous appellerons Raï cette structure algébrique. Posons £& = +00. On a 
EBx=x—miIn(+o0,x), 68 QX = 8 = +00 + x — +00 
ce qui montre que & est l’élément neutre pour la somme et qu’il est absorbant pour le produit 


(e joue le rôle du “zéro”). 0 est l’élément neutre pour le produit. L’équation (1.2) relative au 
processus d’assemblage se réécrit : 


PC) =uAG —-T)Eus(t—T). 
L’équation du limitateur de débit (1.3) devient 
PG)=1@y(—-T)Eu(t—7T) 
(noter que le 1 ne peut pas s’oublier). Sous cette forme, il est clair que ces systèmes sont 


linéaires (i.e. que si (#1, y1) et (#2, y2) sont des couples de solutions, (41 @ u2, y1 ® ÿ2) et 
pour À scalaire, (Au, Ày1) dont également des couples de solutions). 


1.2.2 Systèmes (max,+)-linéaires 


Considérons maintenant l’ensemble RU{—co} muni du max (noté additivement, 263 — 
3 — max(2, 3)) et de l’addition notée multiplicativement (2 @ 4 — 6). Nous noterons Ruax 
cette structure. Soit & — —00.Onas x —xetes @ x, 1e. e est le zéro de R;,. L’équation 
aux dateurs du processus d’assemblage s’écrit maintenant dans Rx : 


y'(n) = T @ (Win) Euz(n)) . 
De même, pour le limitateur de débit, 
J'(n) = rt Qu'(n) ET ® y'(n — 1). 
Ici encore, la linéarité est claire. Les deux structures algébriques Ruax et IRmin Sont des cas 


particuliers de dioïde (structure que nous définirons plus loin et qui fournit le cadre approprié 
pour traiter une classe de systèmes à événements discrets). 
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1.2.3 Description des langages de chemins par des récurrences linéaires 
On considère à nouveau l’exemple du magasin à trois unités de stock 1.1. On note 


{a, b}* _ {_,a,b,aa, ab, ba, bb,aaa,... ,} 

l’ensemble des mots formés seulement des lettres a et b. Nous avons admis le mot vide 
(contenant 0 occurrences de a et b), que nous avons noté “_”. À chaque mot correspond une 
séquence d'événements. On définit le produit de deux mots par simple juxtaposition (par 
exemple, ab @ b = abb). Ce produit est associatif. Le mot vide est élément neutre pour le 
produit. On considère maintenant P({a, b}*) (ensemble des parties de {a, b}*. Un élément de 
P({a, b}*°) est un ensemble de mots (on dit aussi un langage). On munit P({a, b}*) de l’union 
(notée additivement par “@”) et du produit des langages défini par 


L@L'={@l'|leL,l'eL'}. 


Par exemple, {a, b} @ {bb} — {abb, bbb}. Posons & — Ÿ et e — {-}. On vérifie aisément 
que les propriétés de type associativité, commutativité de l’addition et distributivité comme 
en (1.9) sont encore vraies, ce qui justifie les notations & et @. Le langage vide £ est neutre 
pour la somme et le langage e réduit au seul mot vide est neutre pour le produit. On cherche à 
calculer l’ensemble des événements possibles, c’est à dire l’ensemble des chemins du graphe 
1.1 allant du sommet 0 à un sommet quelconque. On introduit à cet effet 


X;; — ensemble des chemins de j à i de longueur n . 


Ona XŸ = e, X, = Ÿ pour i £ j. Au vu du graphe, il est clair que 


X. = X@{a) 

x — X, @ {a} 

Xi = X @{a]® X4 @{b} (1.10) 
XF! = X, @ {a} @ X} ® {b} 


(par exemple, la première équation exprime qu’un chemin qui part de 0 commence forcément 
par a). Autrement dit, l’ensemble des chemins du graphe est donné par un système fini 
d'équations récurrentes linéaires dans la structure algébrique (P({a, b}*), LU, .). En poursui- 
vant les calculs dans (1.10), on obtient entre autres 


Xh = 0 XX, = (a, X), = (IX = x, @{al= (ba), X5, = {bb}. 


X60 = Xa, @ {a} = (Xta} ® X5{b}){a} = {bbaa, baba} , 


ce qui corrobore l'intuition graphique. 


1.2.4 Semi-anneaux et dioïdes 


Le lecteur aura noté que les structures algébriques utilisées ci-dessus vérifient toutes les 
propriétés combinatoires usuelle de la somme et du produit (associativités, commutativité de 
la somme, distributivité), et en outre la propriété plus inhabituelle a & a — a (min(a, a) = 
a,max(a, a) = a,a Ua = a). Cela motive les définitions suivantes. 
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DÉFINITION 1.2.4.1 (SEMI-ANNEAU). On appelle semi-anneau un ensemble D muni d’une 
loi associative (notée additivement “®”), commutative, d’élément neutre &, et d’une loi asso- 
ciative (notée multiplicativement “@”), d’élément neutre e, telles que : 


1. Va,b,ceD,(a@&b)@c—=(a@c) & (b@c) 
2. Va,b,ceD,c@(a@b)=(c@a) & (c@b) 
3. Vae D,a@e=Ee@a—es. 


On a là toutes les propriétés de structure des anneaux, sauf qu’il n’est pas requis que @ 
soit une loi de groupe. 
Lorsque @ vérifie 


a@a—a, (1.11) 


on dit que D est un semi-anneau idempotent, ou dioïde. Si (D \ {&}, @) est un groupe, on 
qualifie (D, &, @) de semi-corps. Un semicorps dont l’addition vérifie (1.11) sera dit idem- 
potent. Un semi-anneau ou semicorps sera dit commutatif lorsque le produit est commutatif. 
On notera ordinairement le produit a @ b par a.b ou ab. 


EXEMPLE 1.2.4.2 (Dioïde des booléiens). L'ensemble {&,e}, muni des lois @ et & 


définies ci dessous est un dioïde commutatif, dit dioïde des booléiens. 


Cette table revient à interpréter & comme “faux”, e comme “vrai”, @ comme ‘ou logique”, 
et @ comme “et logique”. On notera que tout dioïde admet B comme sous-dioïde. 


EXEMPLE 1.2.4.3 (Algèbre (max, +)). La structure R. définie plus haut est un dioïde (et 
même un semicorps idempotent). On l’appelle aussi parfois “algèbre (max, +)”. 


EXEMPLE 1.2.4.4 (Algèbre (min, +)). La structure R,à définie plus haut (traditionnelle- 
ment appelée “algèbre (min, +)”) est également un dioïde. L'application x > —x est un 
isomorphisme de Ruax SUT Rain: 

EXEMPLE 1.2.4.5 (Dioïde des parties de R). L'ensemble des parties de R, muni de l’union 
et de la somme vectorielle est un dioïde. 


EXEMPLE 1.2.4.6 (Dioïde des langages sur {a, b}). La structure (P({a, b}*), U, .) définie 
en $1.2.3 est un dioïde. 


AVERTISSEMENT 1.2.4.7. Semi-anneau est une francisation du terme “semiring” qui est 
bien établi dans la littérature de langue anglaise. Il faut aussi savoir qu’une certaine tradition 
française, due sans doute à Paul Dubreil (l’inventeur des demi-groupes, ou semi-groupes 
comme nous écrirons ici) appelait demi-anneaux les semi-anneaux et réservait le terme de 
semi-anneau à un demi-anneau symmétrisable (1.e. plongeable dans un anneau, par exemple 
N plongeable dans Z est un semi-anneau au sens ancien du terme). Cette tradition semble 
aujourd’hui tomber en déshérence à cause du conflit avec le terme anglais. Nous avons opté 
ici pour l’emploi systématique du terme semi-anneau pour des structures non nécessairement 
symmétrisables. 


EXERCICE 1.2.4.8 (Formule du binôme dans les dioïdes). Donner une formule pour (a & 
b)" dans un dioïde commutatif. Donner une formule plus simple dans le dioïde Rx. 


EXERCICE 1.2.4.9 (suite). (sans l’indication) Soit D un dioïde commutatif dont le produit 
est simplifiable (1.e. ax = ay et a £ € entraîne x = y). Montrer que 
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(a@b}"=a"æb". 


[Indication : on pourra calculer (a @ b)(a & b)° et (a  b)(a? Eb?)] 


Dans un semi-anneau, on peut parler de matrices. Etant données A, B e D'ÂetCe 
Dk*?, on pose 


; k 
ABB); = AB (ASC); = DAx BC. 
1=1 


L'ensemble des matrices carrées de taille n est lui même un semi-anneau. 


He | dans l’algèbre (max, +). On a 


EXEMPLE 1.2.4.10. Soit À — | > 4 


42 | max(G+3,7+2) max(3+7,7+4) | [9 11 
7 | max(2+3,4+2) max(2+7,4+4+4) | | 6 9 


EXEMPLE 1.2.4.11 (Tetris). Pour l’exemple du jeu de Tetris (ou des systèmes à ressources 
partagées), 1l est immédiat de reconnaître dans (1.8) un produit de matrices dans les semi- 
anneau (max, +) 

x(wa) = x(w)M(a) , 


où le coefficient M(a),. vaut h(a), — £(a), sis,r € R(a), M(a),, = 0 quand r & R(a),et 
où les autres coefficients sont —co. 


Présentons, pour l’agrément du lecteur, quelques autres exemples plus ou moins clas- 
siques de semianneaux. 


EXEMPLE 1.2.4.12 (Semianneau tropical et autres semi-anneaux max-plusiens). 
L'ensemble N U {+co}, muni des lois min et +, est un dioïde, appelé semi-anneau tropical, 
en l’honneur de son inventeur, IL. Simon [65]. On notera cette structure Ni. C’est bien-sûr 
un sous-dioïde” de Rain. Plus généralement, à n’importe quel sous-monoïde M de (R, +) 
sont associés les dioïdes Min = (M U {+co}, min, +) et Max = (M U {—co}, max, +). 
AinSi, Zmin = (Z U {+00}, min, +). 

EXEMPLE 1.2.4.13 (Semi-corps de Maslov et Pap). Pour tout réel non-nul p, définissons la 
loi +, : R° x RT — R* : 


a+, b= (a? +bP) | (1.12) 
Évidemment, R; + (R*, +»; X) est un semi-corps, isomorphe à Rj = (R*,+, x), via 
l'application x > x? pour toute valeur de p. Cependant, 


lim a+,b=—max(a,b), 
P—+co 


ce qui nous permet de voir le semicorps idempotent “max-fois” Rf, — (R*, max, x), iso- 


morphe à Rnax (Via x > log x) comme la limite des semi-corps R;. Cette exemple exprime 
donc formellement que l’algèbre max-plus est la limite de l’algèbre usuelle. 


SLes notions de sous-semianneau, sous-dioïde, morphisme de semi-anneau ou de dioïde, sont définies comme 
en algèbre habituelle. Un sous-semi-anneau S’ d’un semi-anneau S est un sous-ensemble de S, contenant les 
éléments neutres € et e, et stable par somme et produit. Un morphisme d’un semi-anneau (S,@, @) de zero & 
et d’unité e vers un semi-anneau (5, @, ®) de zero &’ et d’unité e” est une application f : $ — S’ telle que 
f(@a@æb)= f(a)æ! f(b), f(a &b) = f(a) ® (b), f(e) = e', f(e) = €’. Comme il se doit, un semi-anneau S 
est donc sous-semianneau de S’ s’il s’injecte dans S” et si l’injection est un morphisme de semi-anneau. 
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EXEMPLE 1.2.4.14 (Semi-anneau goulot ou flou). L'ensemble [0, 1], muni des lois max et 
min est un semi-anneau parfois appelé “algèbre goulot”, où algèbre floue. C’est en fait un cas 
très particulier de treillis (voir Chap. 2). Nous parlerons peu de ces structures ici. 


EXERCICE 1.2.4.15. Un élément x d’un dioïde peut il avoir un symétrique (1.e. un élément 
x'tq.x x" —Ee)? 

EXERCICE 1.2.4.16. Cet exercice a été communiqué à l’auteur par Gérard Duchamp. 1. 
Montrer que pour tout p entier positif impair, R muni de la somme +, définie dans 
l’exemple 1.2.4.13, et du produit usuel est un corps, isomorphe au corps des réels. 2. Monter 
que pour tout a,b € KR, la limite a +, b na limz_, 4 +241 b existe, et donner une cha- 
ractérisation simple de la loi +4. 3. Montrer que l’ensemble R, muni de la somme +, et du 
produit usuel satisfait tous les axiomes de structure des corps, sauf un. Lequel ? 


EXERCICE 1.2.4.17 (Compacts convexes). 1. Montrer que l’ensemble des compacts 
convexes de R”', avec pour somme 


A @ B = conv(A U B) 


(l'enveloppe convexe de l’union de À et B) et pour produit 8: A@B—=A+B—{a+b| 
a € À, b € B}, est un dioïde. 2. ({) Montrer que le produit de ce dioïde est simplifiable (1.e. 
AQC=B@CetC Ze — À = B). Indication : on pourra associer à un compact-convexe 
X sa jauge jx : R' — KR, jx(p) = sup,.-x(p, x), et montrer que jxg@y = jx + jy pour tous 
compact-convexes X, Y.3. L'ensemble des convexes (non nécessairement compacts) munis 
des mêmes lois est il encore simplifiable ? 


EXERCICE 1.2.4.18 (Théorème de Cayley-Hamilton). (D Soit D un semi-anneau et une 
matrice À € D"*". On définit le déterminant positif de À par 


det* À = ee) Sd Àio(i) 


o permutation paire à 


et det À de manière analogue pour les permutations impaires, de sorte que l’on a “formel- 
lement” det A — det* À — det” A. On appellera mineur principal positif d’ordre k le det* 
d’une k x k matrice carrée extraite, les indices de lignes étant les mêmes que ceux des co- 
lonnes. Soit A} la somme des C* mineurs positifs principaux d’ordre k et A} la somme des 
C* mineurs négatifs principaux de même ordre. Montrer que, pour n pair° 


AA A" I@AÏA"?6@...@AtId=A}A"! @A;A"?@...@A;Id. (1.13) 


Cette formule généralise le théorème de Cayley-Hamilton. On a classiquement, quand la loi 
@ = + admet une loi opposée, notée —, 


Pa(X) = det(A — XId) = Œ(-e) (AÏ — Ar)X* 
k=0 
et la formule (1.13) exprime l'identité P4(A) = 0 sans signe moins. 


1.2.5 Un exemple de système (min, +)-linéaire en temps continu : le limitateur 
de débit 


Nous montrons maintenant que les phénomènes de saturation et de synchronisation 
aperçus plus haut et représentés par des équations (max, +) où (min, +) linéaires ne sont 


6La formule pour 7 impair s’obtient en changeant A; et A de coté 
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pas spécifiques au cas discret, mais apparaissent également dans des systèmes continus. De 
manière un peu polémique, on pourrait dire que l’important, dans les systèmes à événements 
discrets, ce n’est pas le caractère discret, mais c’est l’algèbre sous-jacente. On aboutit alors 
à une nouvelle notion de linéarité. 

Considérons en effet le système mono-entrée mono-sortie S : u > y représenté à 
gauche de la Figure 1.4. Un fluide circule à travers un long tuyau vers un premier réservoir 


temps de séjour d 


FIG. 1.4: Système continu et analogue discret 


(vide à la date 0). L'entrée w(r) représente la quantité cumulée de liquide entrée dans le tuyau 
jusqu’à l'instant f (la fonction f + u(t) est donc croissante, et l’on a u(t) = 0 pour ft < 0). Le 
liquide met un temps d à parcourir le tuyau. Le liquide passe du premier au second réservoir 
à travers une ouverture qui limite le débit instantané à la valeur maximale de B > 0. On 
note y(r) le volume de fluide dans le second réservoir à l’instant f. On a une quantité initiale 
y(0) — c. Le débit d’entrée dans le second réservoir étant limité à B, on a : 


Vi,WO > 0, y(t +0) < y(t) + B0 . (1.14) 
D'autre part, 
Vt, yG)<ut—-d)+c. (1.15) 
Il en résulte immédiatement que Yf et VA > O, 
YO) <yG—0)+8B0 <u(t—-d—-6)+c+B0, 
d’où pour tout f, 
VE) < Pt —4=0)#c+pOl= nt PERD (1.16) 
Posons 


t<d ; 
Lee POP (1.17) 
c+B(t—-d) sinon. 


et définissons y par : 


50% inf[u(é — +) + K(T)] (1.18) 
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Toute sortie y vérifie y < y. En outre : 
YU+0) — inf [ue (s +0—7Tt)+k(r)] 
TE 
= inf [eee —s)+k(s +0)] 
€ 


S 


inf [ue — s) +k(s)] + BO 
=  y(t) + 0 . 


d’où il résulte que y est la solution maximale de (1.14). C’est cette solution y qui décrira le 
comportement physique du système si l’on suppose que le liquide s’écoule aussi rapidement 
que possible. On a ainsi montré que la sortie au plus tôt du système est représentée par 
une inf-convolution de l’entrée avec la fonction k. De cela, il résulte immédiatement que le 
système u > S(u) est (min, +) linéaire, i.e. vérifie les deux propriétés suivantes : 


min-superposition  Vu,u'e R° S(min(u, v)) = min(S(u), S(v)) 

invariance additive VieR SA +u) = À + S(u). 
REMARQUE 1.2.5.1. On peut voir le calcul qui précède comme la résolution (en fait le 
calcul de la solution maximale) de l’équation fonctionnelle (min, +) : 


y) — min (inf[ y (4 —0)+ BO],u(t —- d)+c). 


REMARQUE 1.2.5.2. Le système à droite de la Figure 1.4 est l’analogue discret du limi- 
tateur de débit. On prend Bh machines travaillant en parallèle et traitant chacun une pièce 
pendant h unités de temps (ce qui fait bien un débit de B). Les pièces doivent attendre un 
temps de préparation d dans la place amont avant d’être traitées. Enfin, un stock initial de c 
pièces (présentes depuis l’instant —oo) est disponible. Il est clair que la quantité de pièces 
servies au plus tôt y(f) vérifie : 


»G) = c+min(Bh + y(t—h),u(t —d)) . 
En substituant y(f — h), et en poursuivant de même, il vient 


YO = € + inf{Bhi + u(t — d — hi)] . (1.19) 


Ceci est une inf-convolution discrète, analogue à (1.18). 


1.2.6 En résumé : systèmes linéaires sur un dioïde et systèmes monotone ho- 
mogènes 


On peut dire que les exemples de systèmes à événements discrets donnés plus haut, à 
l’exception du cas du professeur avisé (Eq. (1.5),(1.6)), relèvent toutes de la notion suivante 
de linéarité. 

DÉFINITION 1.2.6.1. Soit D un dioïde. T — {0}, N, Z ou R le temps (quand 7 — {0}, 
on parlera de système purement statique, dans les autres cas, de système dynamique). Un 
système linéaire mono-entrée mono-sortie (discret ou continu) sur D est une application 
S: Dpt 
qui vérifie les axiomes d’additivité et d’homogénéité : 
Vu,veDT, VAeD, 
(À) SuBv) = S(u)ES(v), 
(H) SA@Su) = ÀA8S(4). 


1.3 Supplément : quelques autres problèmes où apparaissent les dioïdes 19 


La généralisation au cas à ñn entrées et p sorties est immédiate (on a alors un système de 
(D")T dans (D’)7). Le fil conducteur de ce cours est l’étude de ces systèmes “linéaires”. 

On pourrait relaxer l’axiome d’additivité (A), en le remplaçant par l’axiome de monoto- 
nie : 


Vu,ve DT, 
(M) u<v —=  S(u) <S(v), 


où, en anticipant un peu sur le chapitre suivant, on munit D de l’ordre a <b & ab — 
b (quand D = Rx, a < b ssi max(a,b) — b ssi a < b, c’est l’ordre usuel). On voit 
facilement que l’axiome d’additivité (A) entraîne l’axiome de monotonie (M). L'exemple de 
l'enseignant attendant un peu mais pas trop ses élèves (Eq. (1.5) et (1.6)) est justiciable de 
ce cadre : les applications (w’,x'") > y’ ainsi décrites vérifient en effet les axiomes (M) et 
(H). Ainsi que nous le verrons dans la partie du cours consacrée à la théorie spectrale, une 
partie notable des résultats structurels peut s’établir simplement seulement à partir de ces 
deux axiomes, qui semblent bien établir la frontière entre les “bon” systèmes à événements 
discrets (1.e. ceux qui sont analytiquement solubles), et les mauvais. 


1.3 Supplément : quelques autres problèmes où apparaissent les 
dioïdes 


Cette section ne relève pas de la théorie des systèmes à événements discrets proprement 
dite et peut être sautée en première lecture. 


Nous suggérons ici d’autres classes de problèmes pour lesquels des structures de type 
dioïde sont utiles. Ne serait-ce que d’un point de vue historique, il faut rappeler que l’intérêt 
pour ces structures algébriques est venu d’une part de la communauté de la recherche 
opérationnelle (cf. par example [45]) à partir des problèmes de programmation dynamique 
(dont un exemple générique est présenté ci-dessous en $1.3.1), et d’autre part à la suite 
de Maslov, de l’étude de certains phénomènes asymptotiques (inspirés par l’approximation 
quasi-classique de l’équation de Schrüdinger) où des opérateurs linéaires sur le semi-anneau 
(min, +) apparaissent comme limite d’opérateurs linéaires classiques. Il ne saurait être ques- 
tion de présenter ces travaux 1ci (cf. [52]) et l’on s’est contenté de montrer quelques exemples 
où l’algèbre max, + traduit des phénomènes asymptotiques plus élémentaires. 


1.3.1 Problèmes d’optimisation Markoviens en horizon fini 


Soit un système (par exemple un stock) pouvant prendre n états notés 1,... ,n. On 
étudie le système sur N + 1 pas de temps 0,...,N. A chaque pas de temps k, on peut 
décider de modifier l’état à un prix qui dépend de l’instant, de l’état de départ et de l’état 
d’arrivée, soit 


. def À Li ph eee at 
a(k, i, j) = coût de transition de l’état i à l’état j à l'instant k. 


Notons que a(k, i,i) n’est pas forcément nul (ne rien faire pouvant avoir un certain coût). 
S'il n’est pas possible d’aller de i à j, on conviendra d’un coût infini, soit a(k,i, j) — 
+oo. En outre, on se donne un coût final c(ix) associé à l’état in auquel on arrive à la 
dernière itération. On considère le problème d’optimisation suivant : partant d’un état initial 
io, trouver N états consécutifs i1,... ,in qui réalisent le minimum du coût total, soit 
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vo) = min [a(l, io, t1) +a(2, ü,i2) +... + a(N,in_1, in) + CGn)] 
11...I1N 


Le coût optimal v vu comme fonction de l’état initial est aussi classiquement appelé fonc- 
tion valeur. Afin de calculer v au moyen d’une récurrence, on introduit le sous-problème 
d'optimisation de l'instant N — Kk à l'instant N, soit 
def _. ne + ; : 
vG) = min [a(&, i, ie) + a(k+ Ligier) +2. + (Ni in in) + CN) 
1k..lN 

On à vo = v. La fonction valeur du sous-problème vw, se calcule par récurrence rétrograde à 
l’aide de l’équation dite de la programmation dynamique suivante 


Vi, vG) = minfa(k, i, j)+wnOQ)L KZ) (1:20) 
j 
moyennant la condition finale 


Vi, un() = ci) . 


L’équation de la programmation dynamique (1.20) est presque évidente. Elle exprime que si 
Î, Je, Jes1, -.. , Jn eSt une trajectoire optimale, alors jg, je11,... , jn eSt aussi une trajectoire 
optimale pour le même problème sur un pas de moins avec le nouvel état initial j, (‘une sous 
trajectoire T’ d’une trajectoire optimale T est optimale”). On a donc 


ak +1, j, jen) +. +aCN, jn-1, nv) + CON) = HO) , 


et en outre, le choix de j minimise le coût de la trajectoire globale T, soit a(k, i, j)+w:1(j), 
ce qu’exprime (1.20). L’équation de la programmation dynamique (1.20) se réécrit 


() = ak i, j) & mr) (121) 
j 


dans le dioïde Rin, Soit en introduisant la matrice A(k) = (a(Kk, i, j))1<i,;<n et en voyant c 
et v, comme des vecteurs colonnes 


ve = AK) @ vi, UN =C. 
Il en résulte aussitôt que 
v = vo = A(1) Q A(2)@...@A(N)@c. 


Nous retiendrons que le coût optimal est donné par un produit de matrices dans l'algèbre 
min, +. Le contrôle optimal s’obtient en calculant par récurrence vw; par (1.21) et en 


mémorisant à chaque itération le (ou un, le contrôle optimal n’étant en général pas unique) iz 


de : à 2 def 
qui réalise le minimum. Dans le cas où le coût ne dépend pas du temps, on a A(k) = cst = À, 


et donc 
vu = AN @c 


de sorte que l’étude de l’asymptotique de la fonction valeur (lorsque l’horizon N tend vers 
l'infini) se ramène à celle des puissance itérées des matrices dans l’algèbres min, +. Cette 
théorie est abordée au Chapitre 4. 


EXERCICE 1.3.1.1 (Un taxi très déterministe). Considérons un chauffeur de taxi qui circule 
entre 3 villes et un aéroport, comme indiqué sur la Fig. 1.5. Nous supposerons, pour sim- 
plifier, que le chauffeur choisit ses clients et sa destination, les prix des courses étant in- 
diqués sur le graphe. Soit Vy(i) le gain optimal du chauffeur en N courses, sachant qu’il 
part de la ville i. Montrer qu’on Vy satisfait une récurrence sur le semi-anneau max-plus 
Vy = AVy_1, Vo = b, où À est une matrice carrée et b un vecteur colonne, tous deux à 
coefficients dans Rax, que l’on déterminera. Calculer V;,(1). 
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FIG. 1.5: Un taxi très déterministe 


1.3.2 Petits Problèmes Asymptotiques 


L’algèbre (max, +) apparaît naturellement dans certains phénomènes asymptotiques, par 
exemple 


et + else” MGM Jorsque t — +00. (1.22) 


Plus formellement, considérons $ le semi-anneau des combinaisons linéaires à coefficients 
positifs d’exponentielles, soit 


eo % £ 
SJ où & {re S'aexptbi)| &.b; ER) 


k>1 i=1 


et définissons l’exposant maximal d’une combinaison linéaire d’exponentielles s € S par 


def .. 1 
œ{s) — nus 7 HER . 
On a 


Q)  g(s +57) = max(p(s), p(s’)). 


(Gi) (ss) = @(s) + p(s’) (1:23) 


donc @ est un morphisme de semi-anneau de S dans Rx. De manière moins pédante, on 
peut dire que “l’algèbre max, + est celle des exposants” ou l’algèbre des “petits o”, dont le 
lecteur sait bien qu’elle est idempotente (“e + € — €” en analyse). Plus généralement, soit 
pour fixer les idées K € R" un compact et f : K — KR continue. La méthode de Laplace 
bien connue [30] pour l’estimation asymptotique des intégrales montre en particulier que 


lim nu . (| exp dx ) = sup f(x) (1.24) 
[0 K xek 


ce qui étend (1.23),(i) au cas des “sommes infinies”. 
APPLICATION 1.3.2.1. Soit A(r) une matrice dont les coefficients sont de la forme A(f);; — 


exp(ai;t), soit N € N, et posons B = (a;;). Exprimer en fonction de B la limite 
lim ie [AG] 
1—o f él 


(A(1)Ÿ désigne la puissance N ième de la matrice A(r) dans l’algèbre usuelle). 


De manière moins anecdotique, des asymptotiques de ce type sont centrales en physique 
statistique lorsque l’on fait tendre la température vers 0 [21, 20]. 
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Petite bibliographie générale raisonnée pour le cours 


Il y a sur les systèmes à événements discrets, ainsi que sur le semi-anneau (max, +) et ses 
applications quelques références de base, certaines très accessibles, d’autres qui demandent 
un effort sérieux, que le lecteur soucieux de prolonger ce cours pourra fort utilement lire. 

Une excellente référence à jour sur l’approche (max, +) linéaire des Systèmes à 
Événements Discrets (déterministes et stochastiques), ainsi que plus généralement sur le 
semi-anneau (max, +), est “Synchronisation et linéarité”[5], de Baccelli, Cohen, Olsder et 
Quadrat. Ce texte contient des parties très accessibles, d’autres plus techniques, et résume à 
peu près l’état de l’art en 92 sur les SED (max, +). Depuis, ce sont essentiellement les exten- 
sions de cette théorie au cas “non linéaire” (récurrences (min, max, +), systèmes dynamiques 
fluides (min, + x)) ainsi que les extensions aux automates qui ont progressé. On pourra se 
reporter à l’excellente collection d’articles “Idempotency” [48] éditée par J. Gunawardena, à 
paraître incessamment aux Cambridge University Press. 

Une autre référence plus générale sur les Systèmes à Événements Discrets stochastiques 
est le livre de Glasserman et Yao [401]. 

Sur la théorie de Ramadge et Wonham, on pourra lire le petit livre de Kumar et Garg [57] 
mais aussi les papiers originaux [59, 16, 69] qui sont très accessibles. 

Sur le semi-anneau (max, +) proprement dit, une référence classique est “Minimax al- 
gebra” [29] de Cuninghame-Green, qui contient énormément, mais n’est pas toujours facile 
à lire en raison des notations et du vocabulaire qui datent un peu. La théorie a aussi pro- 
gressé depuis. Le livre de U. Zimmermann [71] est très original, et plus orienté “optimisation 
combinatoire”. Le “Graphes et Algorithmes” [45] de Gondran et Minoux contient un cha- 
pitre sur les dioïdes et les algèbres de chemins dans les graphes, qui est lumineux (les autres 
chapitres sont aussi hautement recommandables, c’est vraiment une excellente référence de 
base pour le lecteur intéressé par la théorie des graphes et ses applications). Une collection 
d'articles sur l’ Analyse Idempotente à été éditée récemment par Maslov et Samborskiï [52]. 
Elle donne une idée assez à jour des très importants travaux de l’École “russe”. On y trouvera 
particulièrement des application du semi-anneau (max, +) à l’étude des EDP de Hamilton- 
Jacobi, qui apparaissent en commande optimale, ainsi qu’en analyse asymptotique (méthode 
dite WKB pour la détermination de la phase quasi-classique de l’équation de Schrôdinger). 
La théorie spectrale en dimension infinie (qui étend la théorie de Perron-Frobenius que l’on 
verra dans ce cours) y est traitée. 

Sur les langages formels et la théorie des automates, deux références excellentes mais 
peut-être un peu trop ambitieuses pour un début (ne pas commencer par la sauf à avoir une 
âme d’algébriste) sont le manuel d’Eilenberg [33], ainsi que le livre de Berstel et Reute- 
nauer [9] sur les séries rationnelles. Le livre de Salomaa et Soittola [63] est sans doute plus 
accessible. Nous conseillerons pour un début les survols sur le sujet (ex. Perrin) dans le 
Handbook of Theoretical Computer Science [67]. Le manuel d’Autebert [3] sur les langages 
algébriques est sans doute assez accessible. Pour aller plus loin, voir le livre de Berstel [8]. Un 
merveilleux petit livre, pas toujours très canonique mais vraiment très original, et profond, 
est celui de Conway [28]. 

D'autres références ad hoc seront livrées au fil du cours, et à la fin de chaque chapitre. 


Chapitre 2 


Equations linéaires dans les dioïdes 


Introduction 


Ce chapitre passe en revue la théorie algébrique de base des dioïdes. On s’est limité aux 
résultats utiles pour l’étude des systèmes à événements discrets. L'idée qui est à l’œuvre dans 
cette théorie est qu’à la différence des semi-anneaux généraux, les semi-anneaux idempo- 
tents (dioïdes) peuvent être munis d’une structure ordonnée canonique. On rappelle quelques 
points indispensables sur les structures ordonnées (semi-treillis, treillis). On étudie ensuite 
les équations et inéquations de type 


x=ax®b, x>-ax@b (2.1) 


qui apparaissent comme les équations implicites de graphes d'événements temporisés. Il est 
naturel d’utiliser une technique de point fixe : l’opération étoile (a* = e&a@ ad @...“— 
(1 — a) !”, étoile de Kleene, bien connue en théorie des langages et des séries formelles) est 
introduite à cette occasion. Disons que l’étoile “remplace” l’inverse dans certaines structures 
où l’inverse n’est pas défini. Le fait essentiel pour la suite est que l’étoile d’une matrice 
s’interprète simplement en termes de chemins. En gros, les coefficients de A* représentent 
l’ensemble des chemins du graphe associé à la matrice À : cela permet de relier un point de 
vue combinatoire (chemins) à un point de vue algébrique (l'étoile satisfait tout un formulaire 
utile, que l’on établit). Il y a une autre classe d’équations que l’on sait également étudier, 
celles de la forme 


ax =b . (2.2) 
Disons simplement qu’on montre que la sous-solution maximale, 
a\b = sup{x | ax < b} 


jouït de bonnes propriétés. “a\b” est appelé résiduel de b par a. Cette notion est utile pour 
des problèmes au plus tard, pour lesquels le calcul de la sous-solution maximale revient 
à déterminer les dates de lancement les plus tardives compatibles avec un objectif. Cette 
théorie, qui est une généralisation élémentaire des dualités usuelles et des correspondances de 
Galois est tout à fait instructive, mais non centrale dans ce cours. Nous renvoyons le lecteur 
intéressé à Blyth et Janowitz [12, 4] Signalons également qu’une théorie de type Cramer 
[4, 55, 56, 36], plus délicate que la théorie usuelle, existe pour les équations plus générales 
Ax b=Cx@d. 
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2.1 Dioïdes et structures ordonnées 


2.1.1 Rappels sur les structures ordonnées 


On donne ici le minimum vital sur les structures ordonnées. Cette section contient des 
définitions et des résultats très classiques ainsi que des exemples les illustrant. On pourra lire 
directement la suite, et se reporter à cette section en cas de besoin. 


DÉFINITION 2.1.1.1 (SEMI-TREILLIS). Un ensemble ordonné (E, <) est un sup-semi- 
treillis si toute partie à deux éléments de E admet une borne-sup, 1.e. 


Va, béE?, ixeE, (x-aetx-b)&xxc) 


On notera c = a V b. Définition duale pour un inf-semi-treillis. On notera a A b la borne inf 
de a et b. Si (E, <) est à la fois un sup et un inf semi-treillis, on dit que (E, <) est un treillis. 


EXEMPLE 2.1.1.2. R ordonné usuellement est un treillis. On a a V b = max(a, b), a A b — 
min(a, b). Plus généralement, tout ensemble totalement ordonné! est un treillis. 


EXEMPLE 2.1.1.3 (Treillis des Parties d’un Ensemble). Etant donné un ensemble F, l’en- 
semble des parties de F ordonnées par l’inclusion est un treillis, où les bornes sup et inf sont 
données respectivement par 


AVB=AUB, AnAB=ANB. 


EXEMPLE 2.1.1.4 (Treillis des sous-espaces vectoriels). L'ensemble des sous espaces vec- 
toriels d’un espace vectoriel F est un treillis, avec 


AVB=A+B, AAB—=ANB. 


EXEMPLE 2.1.1.5 (Treillis des compacts convexes). L'ensemble des compacts convexes de 
R?, muni des deux lois suivantes : 


anñhb=anNb, aVb=—conv(aUb) 


où conv dénote l’enveloppe convexe, est un treillis. 
NOTATION 2.1.1.6 (DIAGRAMME DE HASSE). Il est traditionnel de représenter les treillis 
par des diagrammes dits de Hasse semblables à ceux de la Figure 2.1. Dans ces diagrammes, 


[1,1,1] 


[1,00] LOU (0, i,1] 


a b M. his 


(A) (B) {0,1} 


FIG. 2.1: Diagrammes de Hasse 


on trace un arc ascendant de x à y six < y. Les arcs triviaux x + x et les arcs se déduisant 
par transitivité sont omis. Le Figure 2.1,(A) exprime donc que a Ab minore a et b et que a Vb 
majore a et b. On a donné sur la Figure (B) un diagramme de Hasse un peut plus intéressant, 
à savoir celui de l’ensemble des triplets d’éléments de {0, 1}, muni de l’ordre partiel usuel : 
u < vssi Vi € {1,2,3},u; < v;. 


lie. ou l’on a pour tous x et y ou bien x < y ou bien y < x 
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PROPRIÉTÉS 2.1.1.7. Soit T un treillis. Pour tous a, b,c e T,ona: 
1. a V a = a (idempotence) 

a Vb = bV a (commufativité) 

aV(bVc)= (a Vb)Vc(associativité) 

a V (b A a) = a (absorption) 

aV(bAc) <(aVb)A(a Ac) 

6 aA(bVc)> (a Ab)V(a Ac) 


M 8 & ND 


Démonstration. Ces propriétés se vérifient aisément. Montrons par exemple (3). On a 
rue x>->v)Sx>-uVu 
d’où x > a V (b V c)ssi x majore a, b et c. L’associativité en résulte. ( 


Un treillis où les deux dernières inégalités sont des égalités est dit distributif. Par exemple, 
le treillis des parties d’un ensemble (ordonnées par l’inclusion) est distributif. Le treillis des 
sous espaces d’un espace vectoriel F est un treillis qui n’est pas distributif. 


EXERCICE 2.1.1.8. Soit 7 un treillis. Montrer que les deux propriétés suivantes sont 
équivalentes 


1. Va,b,ceT, aV(bAc)=(aVb)A(avVoc) 
2. Va,b,ceT, aA(bVc)=(aAb)V(a Ac). 
Autrement dit, l’une des deux distributivités entraîne l’autre. 
CONTRE EXEMPLE 2.1.1.9. L'ensemble {a,b,c,d} muni de l’ordre donné par le dia- 


gramme de Hasse de la Figure 2.2 n’est pas un treillis. En effet, l’ensemble des majorants 


C d a<c, a <d 
b b<c, b<d 


FIG. 2.2: Ceci n’est pas un treillis 


communs à a et b est égal à {c, d}. Les deux majorants c et d sont minimaux et non compa- 
rables. a et b n’admettent donc pas de plus petit majorant (1.e. de borne sup). 


DÉFINITION 2.1.1.10 (ENSEMBLE ORDONNÉ COMPLET). On dit que l’ensemble ordonné 
(D, <) est complet si toute partie À de D admet une borne-sup, que l’on notera in- 
différemment 


\/A ou Na: 


Lorsque (D, V, A) est un treillis, on dira que D est un freillis complet lorsque toute partie 
admet une borne-sup et une borne-inf. 


On notera? T = \/D le plus grand élément de D. Lorsque qu’il s’agit de l’ordre naturel 
d’un dioïde, on notera la borne-sup @ au lieu de \/. En conséquence, pour une partie infinie 
X, on écrira 


ep pour “top”, terme consacré en anglais 


26 2. Équations linéaires dans les dioïdes 


x L sup X s 


xeX 


Voici sans doute le plus célèbre des résultats sur les treillis complets, d’une simpli- 
cité remarquable. Nous le mentionnons surtout pour la culture, car nous ne l’utiliserons pas 
vraîment dans ce cours. 


THÉORÈME 2.1.1.11 (POINT FIXE DE KNASTER-TARSKI). Une application croissante 
d’un treillis complet T dans lui même admet un point fixe. 


Par exemple, une application croissante de [0, 1] dans lui même admet un point fixe. 


Démonstration. L'idée de la preuve est de montrer que 

a inffteT|t> ft) (2.3) 
est un point fixe (1.e. f (a) — a). On montre d’abord que l’ensemble 

A={GEeTIt2z f()) 
dont on prend l’inf est non vide. Cela résulte du fait que 7 étant complet admet un plus grand 
élément T qui satisfait par définition T > f(T). 
Ensuite, par définition de a, on a 
t> fH)=tza 


donc par monotonie de f 


12 f(1) = f() Z f(a). 
Ainsite At > f(t) > f(a) doncinf{f | ft e A} > f(a), donc a = inf A > f(a) ce 
qui montre une inégalité. D'autre part, en utilisant encore la croissance de f, 

a > f(a)= f(a)z f(f(a)) , 

donc f (a) € À, donc f(a) > inf A = a. (| 
REMARQUE 2.1.1.12. La preuve ci-dessus montre en fait que a est le plus petit point fixe 
de f. On notera bien que l’on a seulement utilisé la moitié des hypothèses, ï.e. le fait que 
T admet un plus grand élément et que tout partie admet une borne inf. On peut montrer 
dualement (1.e. en renversant toutes les inégalités) que 

b sup{t ET | FO > 1} (2.4) 


est le plus grand point fixe de f, 


2.1.2 Dioïdes comme structures ordonnées 


Dans un dioïde (D, &, ©), on définit la relation d'ordre naturelle < par : 
a<bs&adb=b. (2.5) 


Comme l’addition est idempotente, on a a = a @ a > a, ce qui montre la réflexivité. L’an- 
tisymétrie et la transitivité sont claires. Cette relation d’ordre est compatible avec les lois de 
structure de TD), 1.e. : 


a<b=a&c<b@&c 
a <b=ac<bc. 


(D, <) est ainsi un semi-treillis dans lequel la borne-sup est donnée par @ (a @ b est le 
plus petit majorant de a, b) et où & est le plus petit élément. 
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EXEMPLE 2.1.2.1. Dans Ruxx (cf. 11.2.4.3), la relation < coïncide avec l’ordre coutumier. 
Dans Rhin (cf. 1.2.4.4), on a x < y ssi min(x, y) = y, et donc < est l’ordre dual de l’ordre 
usuel (par exemple, 2 > 3). Soit X un ensemble, dans (P(X), U), on a A <B & AUB = 
B & AC B. 


Nous introduisons maintenant la sous classe commode et importante des dioïdes com- 
plets. 


DÉFINITION 2.1.2.2 (DIOÏDE COMPLET). Le dioïde D est complet s’il est complet en tant 
qu’ensemble ordonné par (2.5) et s’il vérifie les deux propriétés suivantes, dites de “distribu- 
tivité infinie” : 


(BDuca a)b = Dix ab 
b(Buea a) = Dis ba . 


VACD, VbeT, 


Il en résulte immédiatement que 


ŒoBn= © «. (2.6) 


acA  beB (a,b)EAxB 


EXEMPLE 2.1.2.3 (Dioïde Complété de Rnax). IR U {oo}, muni du max et du +, avec la 
convention (+00) + (—00) — —0o, est un dioïde complet que l’on notera Ra. 


EXEMPLE 2.1.2.4. L'ensemble des applications de R dans R U {Ho}, muni du max point 
par point et du produit de sup-convolution défini par : 


(f ® 8)4) = suplf( — x) + g(x)] 


xeR 
(avec (+00) + (—00) — —co) est un dioïde complet que l’on notera Ra. L'élément neutre 
pour le produit est la fonction e donnée par e(f) = —co sit 0 et e(0) = 0. 


REMARQUE 2.1.2.5. On peut montrer qu’il y a incompatibilité entre la complétion et l’in- 
versibilité de la loi produit, soit le résultat suivant : un semi-corps idempotent non trivial n’a 
pas de plus grand élément, et en particulier n’est pas complet. Soit en effet oo plus grand 
élément d’un semi-corps idempotent D. On a 60.00 > 00.6 = co, et en simplifiant © = e. 
Ainsi, tout x non nul vérifie 


x <e (2.7) 
D'où en passant aux inverses x—! > e. L’autre inégalité s’obtient en appliquant (2.7) à x7!, 
d’où x —=eet D — {&,e} est trivial. 


NOTATION 2.1.2.6 (INTÉGRATION IDEMPOTENTE DE MASLOV). On notera, à la suite de 
Maslov [53, 50], 


def 
a # Da 
iel 


iel 


par analogie avec l’intégrale usuelle. On a en particulier la règle de Fubini : 


ALT 
Urer {IX JR) ielJ jeJ(i) 


qui n’est autre que l’associativité de la borne sup. L’infinie distributivité du produit se réécrit 
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ce(s a) =# cCQOa. 
iel iel 


La sup-convolution se réécrit comme suit : 
Cf @ 8)4) = suplf( — x) + g(x)] =$ fG —x) @ g(x) 
xeR xeR 


ce qui montre bien l’analogie avec la convolution usuelle. 
Le résultat suivant permet de munir un dioïde complet d’une borne inf. 


PROPOSITION 2.1.2.7. Soient D un dioïde complet, a, b € D. La borne inf de a et b existe 
et est donnée par 


anb=\/{xlx<aetx <b}. (2.8) 


Le seul point non trivial est que l’ensemble à droite de (2.8) est non vide : cela résulte de 
ce que D admet un plus petit élément égal à &. On notera que l’inf d’une partie quelconque 
X C D existe et se caractérise de manière analogue : 


Ax=\V6e?D] VrexX, y<x}. 
On a une autre classe importante de dioïdes où la borne inf est définie. 
PROPOSITION 2.1.2.8. Soit D un semi-corps idempotent. Toute partie finie admet alors une 
borne inf. En outre, la borne inf distribue par rapport au produit, i.e. 
(a Abjc—=acAbc, c(aAb)=caAcb. (2.9) 
Dans la littérature, (cf. Dubreil [31]), la structure (D \ {&}, @) s’appelle aussi groupe 
réticulé ou groupe ordonné en treillis. 


Démonstration de la Proposition. Supposons x, a,b non nuls. Les propositions suivantes 
sont équivalentes : 


x<a et x <b 
x That et xt <br! 
xt EG tot 
< @l@b 
< b(a@b) la. 
On en déduit que sia@b Æ£eona 
anñnb=b(a@b) la. (2.10) 


Sia = eetb = 6, trivialement a À b = &. La distributivité résulte de ce que x < ac et 
x < bc est équivalent à xc7! < (a A b). (| 


CONTRE EXEMPLE 2.1.2.9. Le dioïde B[X] des polynômes à coefficients booléiens en une 
indéterminée X est un exemple de dioïde où la borne inf est bien définie, et qui n’est pas 
complet (considérer Gex X! 4 B[X]), et où le produit est loin d’être inversible (il n’est 
même pas simplifiable, par exemple, (e & X)(e & X°) = (e & X)(e  X & X°)). Il faut 
bien voir que les conditions exhibées ci-dessus sont des conditions suffisantes commodes 
pour l’existence d’une borne inf (connaissant la borne sup), mais ne sont absolument pas 
nécessaires. 
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CONTRE EXEMPLE 2.1.2.10. Dans des dioïdes généraux, il est faux que l’inf (quand il 
existe) distribue par rapport au produit. Considérons par exemple le dioïde des parties de 
R, muni de l’union et de la somme vectorielle. La borne inf n’est autre que l’intersection. 
Soient À = {0}, B={1},C —=R Ona 


(AAB)@C=(ANB)+C=H+R=HZ(A@C)A(B@C)=RNR=R. 


Cela montre que la borne sup @ et la borne-inf A jouent en général des rôles asymétriques 
dans les dioïdes. 
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2.2.1 L’Équation x — ax @ b dans les Dioïdes Complets 


On considère l’équation x = ax @ b. En posant f(x) = ax @ b et en notant que f est 
croissante, on peut appliquer le théorème du point fixe de Tarski et en déduire l’existence 
d’une solution. Cela ne nous satisfait pas cependant, parce que les formules (2.3) et (2.4) 
donnant le plus petit et le plus grand point fixe sont peu effectives. Ici, comme il s’agit d’une 
équation linéaire, on peut caractériser plus précisément le plus petit point fixe à l’aide d’une 
nouvelle opération fondamentale : l’éfoile. 


THÉORÈME 2.2.1.1. L'inéquation x > ax @b dans un dioïde complet admet une plus petite 
solution, égale à a*b, où a* est définie? par : 


a* — Da . 
neN 
En outre, X — a*b réalise l’égalité (c’est-à-dire x — ax @ b). 
Démonstration. On a par une induction immédiate à partir de x > ax @ b que 
x>(eDa®...@a")b@a"tlx, (2.11) 


d'oùx >(eBa@...@a")b — Di a*b pour tout n. On a montré 


x>ax@b= xx a*b= ab : (2.12) 
kEN 


Réciproquement, le fait que x — a*b est solution de x — ax @ b résulte aussitôt de la 
propriété de distributivité infinie 2.1.2.2. [] 


EXERCICE 2.2.1.2. Donner une formule analogue pour la plus petite solution de x > axc@ 
b (on ne supposera pas le dioïde commutatif). 


NOTATION 2.2.1.3. L'opération “plus” dérivée de l’étoile sera utile : 


+ def 


at = ada D 04 6... —t4a =4'a; (2.13) 
On a 
eat =a*. (2.14) 


3 parfois appelée “étoile de Kleene”. Dans l’algèbre usuelle, x* n’est autre que (1 — x)7!, du moins dès que 
Ix| < 1. 
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EXEMPLE 2.2.1.4. Dans Rax, On à d* = +oo sia > Oeta* = 0 si a < 0. Dans Ruax, 
l’étoile de a n’est donc définie que si a < 0. Cela est reflété par le fait que si a > 0, l'équation 
x = max(a + x, 0) n’admet pas de solution finie. 


EXEMPLE 2.2.1.5. Considérons le dioïde des parties de R, (P(R), U, +). Soit a = [2,3]. 
Onaa?=a+a—T[4,6],...,d’où 


a* = {0} U[2, 3] U[4, 6] U[6, 9] UT8, 12]U...= {0} U[2, 3] U[4, +ool . 


Nous établissons ci-après les propriétés de base des opérations étoile et plus. 
PROPOSITION 2.2.1.6. Soit D un dioïde complet, a, b € D, et c un entier strictement positif. 
On a 

1. (a*)* = a*, 
(a*)* = a, 
a*=(a@...æat)(a}" 
(a @ b)* = (a*b)*a*, 
. (a b)* = b*(ab*)”, 
a* = a*a*. 
(ab*)* = a(a &b})* 
. (ab*)* =e&a(a&b} 


D N D Li À  N 


En outre, lorsque a et b commutent : 


9. (a @ b)* = a*b* 


Démonstration. (1).x > ax @e entraîne x > a(ax Be) Be > a?x Be et plus généralement 
x > a"x @e, soit en sommant x > a*x @ e. On en déduit que la plus petite solution de 
x > ax @ e, 1.e. a*, est plus grande que la plus petite solution de x > a*x @ e, 1.e. (a*)*. Par 
ailleurs, trivialement, a* < (a*)*, d’où l'égalité. 
(2) : a* < (a*}" < (a*)* = a* (par (1). 
(3) résulte du théorème de division euclidienne. En effet, en développant (a° 6... 
a°—1)(a°}*, on obtient la somme des a/+% avec0 <i<c—letjeN. 
(4),(5). Soit x = (a Bb)“. On a 

x =(aæb}x@e (Théorème 2.2.1.1) 

X=ax@bx@e 

x > a*(bX @e) (propriété (2.12)) 
x > (a*b)*a* (idem) 

donc (a @ b)* > (a*b)}*. Pour montrer l’autre inégalité, comme % est la plus petite solution 
de x > x(a @b)@e, il suffit de voir que x’ = (a*b)*a* satisfait cette dernière inéquation. Or 


(a*b)‘a*(a&b)@e — (a*b)'a® @(a*by'a'b&e 
— (a*b)*'at @(a*b)* &e 
(a*b)*a” @ (a*b}* = (a*b)'a* 


Arrivé à ce point, le lecteur peut s’interroger légitimement sur le caractère pénible et un 
peu miraculeux de cette seconde partie de preuve. Nous verrons à la fin de la section 2.2.2 
(paragraphe : “autres preuves du formulaire”) un raffinement qui permet d'économiser la 
vérification que x’ est solution de x’ > x'(a Bb) Be. 


2.2 Algèbres de Chemins 31 


La formule (5) s’obtient à partir de (4 par symétrie. 


(6) : par (4), a* = (a & a)* = (aa*)*a* = (a*)*a* = a*a* (par (2)). 

(7) Via (5), a(a @ b)* = ab*(ab*)* = (ab*})*. 

(8) : résulte de (7) et de la formule (2.14) ci-dessus. 

(9) : On a alors pour k > 1 : (a*b}* = (a*b* = a*bf (par (6)). D'où (a Bb)* = (a*b)*a* — 
(e & D, a*bl)a* = a* & D, a*bf = a*b*. {me 


EXERCICE 2.2.1.7. Parmi les formules de la proposition 2.2.1.6, dire celles qui restent 
vraies lorsque a, b sont des éléments d’une algèbre de Banach, avec [al] < 1, |b|| < 1, 
et [la + b|| < 1 (on définit toujours a* = a° + a + a? + :..). Pour une réponse sans calcul, 
voir la discussion à la fin de la section 2.2.2. 


EXERCICE 2.2.1.8. On se place dans le dioïde Rinax [IX] (séries formelles à coefficients 
dans Rmax). Calculer les séries s1 = (2 @ X})* et s> = (X? @ X°)*. 


EXERCICE 2.2.1.9. Montrer que Rx [[X]] n’est pas complet. Soit s — Ben SnX" € 
Rinax LLX]1 (série formelle à coefficients dans le dioïde R:,). Montrer que si so — &, s* 
existe. Etudier la réciproque. 


2.2.2 Cas particulier des dioïdes de langages et interprétation du formulaire 


Nous rappelons ici quelques définitions très classiques, déjà esquissées dans l’introduc- 
tion (cf. $1.2.3). Les dioïdes de langages sont des cas particuliers de dioïdes complets, qui 
fournissent une interprétation très naturelle du formulaire ci-dessus. 

Soit À un ensemble (traditionnellement appelé alphabet, et dont les éléments seront 
appelés lettres). Un mot de longueur k est un k-uplet w = (a,... ,a) € A*. On le notera 
par simple concaténation : w — a1...a4. Par exemple, abaa est un mot de longueur 4 sur 
l'alphabet A = {a, b, c}. A° désigne conventionnellement l’ensemble des mots de longueur 


66 °° 


0, comprenant le seul mot vide, que l’on notera e ou “_”. 
A AU AU AU. 


désigne l’ensemble des mots (de longueur arbitraire mais finie). Etant donné w € Æ*, on 
notera |w| la longueur de w (ï.e. l’unique k tel que w € Af). Etant donnés deux mots w = 
a...a € À et v = bo...b, € À’, on définit le produit (dit de concaténation) u = wv = 
a...@&bo...b, € AT. Par exemple, (ab)(abb) — ababb. Le mot vide e est élément neutre 
pour ce produit, qui munit 4* d’une structure de monoïde“, dit monoïde libre” sur À. 


DÉFINITION 2.2.2.1 (DIOÏDE DES LANGAGES DE 4*). Une partie de 4* (un ensemble de 
mots) sera appelée langage. Le produit de deux langages L et L' est défini par 


EL'=NEl EMeErSLT 
Par exemple, 


{a, b°?} @ {c} = {ac, bc} . 


ÆRappelons qu’un monoïde est un ensemble M muni d’une loi @ associative et admettant un élément neutre 

SD'un point de vue naïf et informel, le terme “libre” signifie que les éléments de 4* ne vérifient aucune 
“relation particulière” (du style ab? = ba?). De manière précise, posons À = {a1,... ,a)} (pour la simplicité, 
on suppose ici que À est fini). Alors 4* est défini de manière unique (à un isomorphisme de monoïde près) par 
la propriété universelle suivante : pour tout monoïde M et pour tous x1,...,xp € M, il existe un morphisme de 
monoïdes @ : A* — M tel que Vi € {1,...,p}, @(ai) = xi. 
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L'ensemble des langages$ P(4*), muni de l’union et du produit est un dioïde. Il a pour zéro 
l’ensemble vide Ÿ et pour unité le singleton {e} (e note le mot vide). 


On notera que l’étoile d’un langage L n’est autre que l’ensemble de tous les concaténés 
des mots de L, par exemple 


{a,ab}* = {a, ab, aa, abab, aab,aba,...}. 


En particulier, on retrouve que {a, b}* est l’ensemble des mots sur À — {a, b} ce qui est 
cohérent avec la notation 4* pour le monoïde libre. Nous ferons la convention suivante qui 
permet de simplifier les formules : 


$ Dans la suite, on identifiera le mot w € A* avec le singleton {w} € P(4*). Cela nous 
permet d'écrire par exemple a @ aba = {a} U {aba} = {a, aba). 


Application : autres preuves du formulaire On peut donner une preuve alternative bien 
facile du Formulaire 2.2.1.6 : il suffit de montrer que ces identités rationnelles sont valides 
dans le dioïde des langages sur {a, b}, autrement dit, “si c’est vrai pour les langages, c’est vrai 
dans n’importe quel dioïde complet”. Pour établir que ces identités sont vraies dans le dioïde 
des langages P({a, b}*), on peut user de deux techniques différentes. Prenons l’exemple de 
l'identité (a @ b)* = (a*b)*a*. 

La première technique de preuve est purement combinatoire : à gauche de d’égalité, on 
a la somme de tous les mots sur l’alphabet {a, b}. Aïnsi, cette égalité ne fait que dire qu’un 
mot quelconque en a, b se factorise sous la forme (a”!b)(a"2b)...(a"*b)a, qui est le terme 
générique du développement de (a*b)*a*. 

La seconde technique consiste à réecrire l’argument d’élimination de la preuve de la 
Proposition 2.2.1.6, cette fois ci dans P({a, b}*). Pour tout langage L € {a, b}* ne contenant 
pas le mot vide, et pour tout langage L' C {a,b}", il est aisé de montrer que l’équation 
x — Lx @ L' admet solution unique dans P({a, b}*) (à la main ou par un argument de 
contraction, voir l’annexe A.exercice A.8). Puisqu’il y a solution unique, on peut raisonner 
par équivalence : x = (a@b)* & x =(aBb)x®Be & x = a*(bxBe) & x = (a*b)*a*. 
On fait ainsi l’économie d’une moitié de la preuve. 

Le lecteur pourra s’exercer à re-prouver ainsi le reste du formulaire 2.2.1.6. La première 
technique (combinatoire) à l’avantage de révéler la différence de nature entre ces différentes 
formules. L'identité (a @ b)* = (a*b)*a*, sa duale, ainsi que a* = (a° &@..: @ af! )(a°)* 
sont des tautologies (on peut établir une bijection entre les termes des développements des 
deux membres de ces égalités). Ces trois identités, ainsi que (ab*)* — a(a &b)* et (ab*)* — 
e @ a(a @ b})* qui se déduisent de (a & b)* — (a*b)*a*, restent donc valides dans toute 
structure où l’on peut donner un sens raisonnable a l’étoile. Le lecteur pourra vérifier par 
exemple que ces formules restent vraies en posant x* — (1 — x) _!, lorsque a,b, c sont 
des éléments d’un corps, et que les inverses existent. Par contre, les preuves des formules 
a*a* = a*, (a*)* = a*, et dans le cas où ab — ba, (a & b)* — a*b* font intervenir de 
manière essentielle l’idempotence de D. 


6Dans tout ce texte, on note P(X) l’ensemble des parties d’un ensemble X 

TDe manière précise, si D est un dioïde complet et a’, b' € D, il existe un unique morphisme @ de dioïde 
complet (i.e un morphisme préservant les sommes infinies) de (P({a, b}*), U, -) tel que @(a) = a’, g@(b) = b’. 
(c’est dire que (P({a, b}*), U, -) est le dioïde complet libre sur {a, b}). Le morphisme envoie une identité 
rationnelle valide dans P({a, b}*) sur une identité rationnelle valide dans D. 
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2.2.3 Équations implicites matricielles 


On considère toujours l’équation 
x = Ax @b (2.15) 


A étant ici une matrice carrée et x, b étant des vecteurs, tous à coefficients dans un dioïde 
complet D. En se plaçant dans le dioïde des matrices carrées D"*", et en complétant b et x 
en des matrices carrée par des colonnes de &, on se ramène à appliquer le Théorème 2.2.1.1, 
et donc la plus petite solution de (2.15) est donnée par x — A*b et vérifie en outre l’égalité. 
Cependant, l'expression de x — A*b — @, A"b n’est pas effective. On cherche ici un 
algorithme permettant de calculer les étoiles de matrices dès que les étoiles des scalaires sont 
connues. Dans l’algèbre usuelle, c’est l’algorithme de Gauss qui fait le travail (on inverse 
la matrice si on sait inverser les pivots). L’algorithme de Gauss s’adapte dans les dioïdes. 
Le point clé est le Lemme d’inversion par bloc suivant, qui est une version ‘“‘dioïde” d’une 
formule bien connue dans le cas des corps. 


LEMME 2.2.3.1. Pour la matrice suivante partitionnée en quatre blocs, on a 


a | b | __| a* @ a*b(ca*b & d)*ca*  a*b(ca*b & dy” (2.16) 
LT 4 (ca*b @ d)*ca* (ca*b & d)* ‘ 
Démonstration. Soit 
= | X11 X12 | 
X21 X22 
A* est la plus petite solution de X > AX @ Id, soit 
Xn > 4x Bbhxu Be 
X12 > 4X12 ® bxn 2.17) 
Xi > Cx11 ® dx 
Xn2 > Cxy Edxn Be . 


La seconde équation permet d’éliminer x12 : 
X12 > 4“bxn . 
En substituant x, dans la dernière équation, il vient 
X22 > ca*bxn Edxn @e, 
d’où 
X22 > (d & ca*b)* ; 


On obtient ainsi le coefficient (2,2) de (2.16). Les autres coefficients s’obtiennent par un 
argument analogue. Jusqu'ici, on a seulement montré que A* est supérieur ou égal au second 
membre de (2.16). Pour montrer l'égalité, il faut reporter ce second membre dans X = AX@ 
Id et vérifier qu’il satisfait l’égalité, ce qui est une vérification de routine (instructive) laissée 
au lecteur (le lecteur paresseux pourra employer la seconde technique de preuve à la fin de la 
section 2.2.2 pour s’en abstenir). . Ü 
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COROLLAIRE 2.2.3.2. On a 
(ca*b @ d)* = d* & d'c(bd*c & a)*bd* . (2.18) 


Démonstration. On obtient en effet cette formule pour A4, si l’on élimine le système dans 
un ordre différent (l’expression à droite de (2.18) n’est autre que A, dans (2.16) à un chan- 
gement de lettres près). Ü 


Le lemme 2.2.3.1 appliqué de manière inductive permet de calculer l’étoile d’une ma- 
trice quelconque, et l’on peut formuler un algorithme analogue à l’algorithme de Gauss. Nous 
pourrions ici le faire directement, mais préférons faire appel à une interprétation combinatoire 
qui donne une vue très simple de ces questions. De même que le Formulaire (scalaire) 2.2.1.6 
de l’étoile devient immédiat si on l’interprète en termes de mots, de même, les formules de 
type inversion par blocs s’interpréteront très facilement en termes de graphe. 


2.2.4 Interprétation combinatoire de l’algorithme de Gauss 


DÉFINITION 2.2.4.1 (MATRICE ASSOCIÉE À GRAPHE). On rappelle qu’un graphe 
(orienté) est un couple de deux ensembles (S, A), avec À € S?. S est l’ensemble des 
sommets, À l’ensemble des arêtes. On pourra repérer chaque arête (i, j) par une lettre 
distincte d’un alphabet 4. A un graphe à » sommets, on associe la matrice M € P(4*)"”" 
telle que 

à » lettre valuant l'arc ji s'ilyaunarcjr i 


iÿ = j 
E sinon. 


FIG. 2.3: Graphe et matrice associée 


L'exemple d’une telle matrice et de son graphe est donné sur la Figure 2.3. 

On appelle chemin de longueur k de j à i une suite p = (ig,ix_1,... , 11) avec i, = j 
et i, — i telle que pour tout ! = 1,...,k —1,i, + ü11 Soit un arc. Sii — j, on parle de 
circuit. Le poids d’un chemin p est défini par 


w(p) = Mi: &...@ M;;i 


Quitte à identifier le mot w(p) et le singleton {w(p)}, on pourra voir M comme une matrice à 
coefficients dans P(A*). Le théorème suivant est fondamental : il traduit le fait que multiplier 
des matrices revient à concaténer des chemins. 


THÉORÈME 2.2.4.2. M° est égal à la somme des poids des chemins de longueur k allant 
de j à i. M est égal à la somme des poids des chemins du graphe de j à i (de longueur 
arbitraire). 


Démonstration. La propriété pour les M entraîne immédiatement la propriété pour M*. 
Nous montrons seulement la propriété pour M?, le cas général n’étant pas plus difficile. On a 


M5 = D MM; . 
k 
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MM; est non nul ssi il y a un arc de j à k et de k à i, 1.e. ssi (i, k, j) est un chemin de 
longueur 2 de j à i, auquel cas w(i, k, j) = Mix M; est le poids de ce chemin. CO 


s 


On notera qu’avec la définition de M, les mots se lisent de droite à gauche (par 
exemple, cb est un chemin 1 — 2). Il faut bien voir que l’intérêt de la définition ‘trans- 
posée” de M;; = lettre j + à (ij se lit de la droite vers la gauche) est d’être cohérente 
avec l’écriture usuelle des produits de matrices de droite à gauche. Par exemple, si X 
est une matrice telle que X;; représente des chemins j — à, alors la matrice X’ = MX 
rajoute les arcs de M à la fin des arcs de X, ce qui se révélera utile dans les équations 
matricielles standard de l’automatique, -typiquement x(n) — Ax(n — 1) - où les ma- 
trices opèrent à gauche. Avec la convention non renversée pour M, il faudrait faire des 
produits de matrice de gauche à droite, c’est d’ailleurs ce que font les probabilistes 
pour les chaînes de Markov où les mesures sont représentées par des vecteurs lignes 
(équation de Fokker-Planck), mais allez donc écrire une représentation d’état inversée 
x(n +1) =x(n)A u(n)B, y(n) = x(n)C devant un automaticien orthodoxe 


La variante suivante de l’algorithme de Gauss est classique (cf. Backhouse & Carre [7], 
Gondran & Minoux [45l) : 


ALGORITHME 2.2.4.3 (DE JORDAN). Soit À une matrice n xn à coefficients dans un dioïde 


complet D, et A, ..., AU) les matrices définies par : 
AO = À 
pilier, APE NAS NCA PA" (2.19) 
On a At) = At, 
On rappelle que A* A A7 @ A... — AA* = A*A et que A* = Id & At 


(les calculs de A* et de A* sont pratiquement équivalents). Le terme A sera qualifié de 
k-ième pivot par analogie avec l’algorithme de Gauss usuel. 


EXEMPLE 2.2.4.4. Appliquons l'algorithme de Jordan à un scalaire a. Il vient 
aO a ,aD = 0 GaV (a )aO = 4 Daa*a =apa Sa ®...=at. 


EXEMPLE 2.2.4.5 (Enumération des Chemins). On illustre l’algorithme de Jordan 2.2.4.3 
sur le graphe à 2 sommets représenté sur la Figure 2.4. Soit 


FIG. 2.4: Enumération des chemins 


dy 4 
= nu 40 | 
d21 An 
la matrice A* représente alors les chemins non triviaux du graphe. On a en appliquant 


2.2.4.3 : 


+ * * * * * * * 
x=at=| @ 4j ,di2(an ® ana; dn)'anai, Ajjdi(an ® ay4ï412) (2.20) 
7 A * * * , ’. 

(an ® ana\,a)"ana\, (an ® ana,an)t 
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< + 
et en développant 4;, 
+ 
A5 = 422 ® dnid12  dr1d11412 P Ar421Q12 D A21Q12Q22 P .… 


ce qui est bien l’ensemble des chemins de 2 à 2. Au passage, on a retrouvé la Formule (2.16). 


Démonstration. (de l’ Algorithme). On considère les coefficients de À comme des 
indéterminées, ce qui revient à se placer dans le dioïde des langages sur l’alphabet À = {a;; | 
1 <i,j < n}.lIl suffit de voir que le résultat est vrai “formellement” dans (P(.4*), U, -) 
pour qu’il soit vrai dans n’importe quel dioïde complet. La preuve donnée ici repose sur une 
interprétation des langages produits à chaque itération par l’algorithme de Jordan comme en- 
sembles de chemins du graphe associé à la matrice A. Nous renvoyons le lecteur à [7, 44] 
pour une autre preuve de nature plus algébrique de l’algorithme de Jordan. Appelons ici 
chemin de longueur k de j à i un mot de la forme 


P = dii,di,is +. dis ji 
avec i; = i etiy — j. Cela revient à identifier le chemin un chemin p — (iiisisig) au 
sens de 2.2.4.1 avec son poids w(p). On dira que le chemin passe exclusivement par Z si 
l,... ,ix-y € Î (mais tous les à € JZ ne sont pas forcément atteints). Le coefficient a;; 
s’interprète comme le chemin de longueur 1 allant de j à i. Le lemme suivant donne une 
interprétation combinatoire fort utile de l’algorithme de Jordan. 


LEMME 2.2.4.6. A) est égal à la somme des chemins de longueur au moins 1 de j à i 
passant exclusivement par les sommets 1,...,k. 


Démonstration. Dans (2.19), le second terme AË Pat Dikat ni 


; S’interprète comme la 
somme des chemins de j à k puis de k à k, puis de k à à, ces chemins passant par ailleurs 
exclusivement par 1,...,k — 1, soit la somme des chemins de j à i passant au moins une 
fois par k et exclusivement par 1,...,k — 1. Les chemins ne passant pas par k ont été déjà 


énumérés dans le premier terme de (2.19). La récurrence en résulte. ( 


Finalement, A est égal à la somme des chemins de longueur au moins 1 de j à à, 1.e. 
A = Aë. Cela achève la preuve de 2.2.4.3. Ü 


Ces résultats nous permettent de donner une propriété algébrique vérifiée par les chemins 
des graphes. Nous dirons qu’une partie de À de 4* est rationnelle ssi elle peut s’écrire à l’aide 
d’une formule finie faisant intervenir seulement les opérations &, ©, x, ainsi que des parties 
finies de À. De manière plus formalisée : 


DÉFINITION 2.2.4.7 (PARTIES RATIONNELLES). Le dioïde des parties rationnelles de 4* 
est le plus petit sous dioïde de P(.4*) contenant les parties finies et stable par @, @, *. 


Le lecteur notera l’analogie avec les fractions rationnelles classiques : les parties finies 
remplacent les polynômes, et l’étoile remplace l’inverse. Le résultat suivant forme une moitié 
du théorème de Kleene qui sera énoncé au chapitre S. 


COROLLAIRE 2.2.4.8 (FONDAMENTAL). L'ensemble des chemins entre deux sommets 


donnés d’un graphe fini est une partie rationnelle. 


Démonstration. Résulte par exemple de l’algorithme de Jordan 2.2.4.3 qui calcule AÏ à 
partir des coefficients de la matrice À à l’aide d’un nombre fini de sommes, produits, et 
étoiles. [] 


BNoter ici encore que le chemin se lit de la droite vers la gauche. 
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2.2.5 Chemins de poids extrémal 


On considère un graphe dont les arcs sont valués par des réels. Soit AÀ;; la valuation de 
l'arc j + à (noter la transposition). Le poids d’un chemin p = (i1,... ,i,) est ici défini par 


wWaA(p) = Ai +... + Aix | 


On cherche le poids minimum d’un chemin de j à à. Il suffit de penser à des problèmes de 
trajets routiers de longueur minimale, ou de durée minimale pour se convaincre de l’intérêt 
de se problème. Les problèmes de programmation dynamique du type évoqué en 1,$1.3.1 
rentrent également dans ce cadre lorsque le coût ne dépend pas du temps. 

On se place dans le dioïde Rain (.e. (RU {—00.+ 00}, min, +)), où le zéro est £& = +00. 
Lorsque j ++ i n’est pas un arc, on posera À;; — €. 


PROPOSITION 2.2.5.1. À; est égal au poids minimum d’un chemin de j à i. 


Démonstration. Résulte immédiatement du théorème 2.2.4.2 : prendre la somme des poids 
des chemins, c’est prendre le minimum de ces poids. Ü 


Dans le cas où À est à coefficients dans le dioïde Rain et non Riin (ce qui exclut la valeur 
—0o pour les coefficients et correspond au cas où les valuations des arcs du graphe de À sont 
finies), il est naturel de se demander à quelle condition A* appartient à R,, (ce qui signifie 
que le poids des chemins du graphe est inférieurement borné). 


PROPOSITION 2.2.5.2. Soit AE Ri*". L'étoile A* converge dans RŸ ssi il n’y a pas de 
circuit de poids négatif dans le graphe associé à À. On a alors 


A*=Id@A@...@A"! (2.21) 


Démonstration. Soit un chemin p de j à i. Si p est de longueur au moins ñn, p passe au 
moins deux fois par le même point. On peut donc factoriser p — p1C1p2... pxCxprr1 Où les 
c; sont des circuits élémentaires (éventuellement triviaux) et où p1p2 ... p.11 est un chemin 
élémentaire de j à i, comme dans le petit dessin suivant. 


On a, comme le poids ne dépend que de l’image commutative du chemin 


W(p) = W(Pp1P2... Per) @ W(C1...Cx) < W(P1P2... Pkr1) 


d’où il résulte que le min dans 4°: est atteint pour un chemin élémentaire (de longueur au 
plus égale au nombre de sommets du graphe diminué d’une unité, c’est-à-dire n — 1), d’où 
À; = Id;; &4;; E...@ A Réciproquement, si c est un circuit (disons le longueur k) de 
poids strictement négatif et passant par i, on a 


Aÿ < (w(c)}" + —00 


ce qui montre que A* ne converge pas. Ü 


En reprenant le premier point de la preuve ci-dessus, on montre qu’il existe un circuit 
de poids négatif ssi il existe un circuit de longueur au plus n de poids négatif, ce qui est très 
facile à vérifier. Par exemple : 
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LEMME 2.2.5.3. Il y a un circuit de poids négatif dans le graphe associé à A ssi tr(A @ 
... @ À") < 0, où la trace d’une matrice est définie comme de coûtume : t M — @, Mi. 


D'un point de vue pratique, on ne calcule pas Id & À & A? @...@ A"7!, mais plutôt la 
suite définie par la récurrence 


Xo=Id, Xe11 = AXx OU . 


On a X4 = Id &...@ Af, et l’on s’arrête dès que Xy41 = Xx (= A* nécessairement) 
où que trXx; — Œ;(Xx)x < 0 (circuit de poids négatif). Cette méthode exige en général 
un temps O(n*). On peut aussi calculer par élévations successives au carré (Id & A avec 
2 > n (temps O(n°logn)) ou appliquer l'algorithme de Jordan (temps O(n*)). Lorsque 
l’on veut résoudre x — Ax @ b pour un vecteur b fixé, il y a évidemment des algorithmes qui 
calculent directement le vecteur A*b sans calculer A*, ce qui serait maladroit ou infaisable 
(par exemple lorsque À et b sont de grande tailles et creuses, il est hors de question de calculer 
A* qui est pleine). Pour résoudre x = Ax@b, l'algorithme le plus simple est comme ci-dessus 
l’itération de point fixe xo = b,x1 — Axo @ b, etc. On voit facilement que, en notant n la 
dimension, x, — x1_1 — A*b si A*b n’a pas de composante —co, et que x, < x,_1 (1e. 
Xn > Xn_1) Sinon. On a ainsi un algorithme en temps O(n°) pour calculer A*b (ou O(nE) si 
À est codée de manière creuse et si E désigne le nombre d’arcs du graphe associé). Lorsque 
les coefficients de À sont négatifs ou nuls, un algorithme de type glouton, du a Dijkstra, 
permet de calculer une colonne de A* (c’est le problème du calcul des poids minimums d’un 
sommet i fixé à tous les autres sommets du graphe) en un temps O(n?) (ou O(n + E), i.e., 
en temps linéaire), à comparer avec le O(n*) du cas général. Il y a de multiples variantes de 
ces algorithmes de chemins, qui sont examinées en détail dans [45]. 


EXERCICE 2.2.5.4. Calculer l’étoile de la matrice suivante à coefficients dans le dioïde Rrin 


£ 2 3 
A=| —1 € —2 
5 3 72 


Notes et références bibliographiques 


Le lecteur intéressé par la théorie des treillis pourra se reporter aux livres de Birkhoff 
[10] et de Dubreil-Jacotin, Lesieur et Croisot [31]. Sur la résiduation, la référence de base 
est Blyth et Janowitz [12]. Les treillis sont étudiés en détail d’un point de vue topologique 
dans [38]. Sur les systèmes d’équations linéaires dans les dioïdes, leur lien avec la théorie 
des graphes, la meilleure référence est Gondran et Minoux [45]. Citons aussi Carre [19]. On 
pourra aussi se reporter à Cuninghame-Green [29], à Cao, Kim et Roush [17], à Zimmerman 
[71] ainsi que, pour une présentation plus récente et plus orientée Systèmes à Événements 
Discrets, à Baccelli et al. [4]. 


noter qu’il s’agit bien de n et plus de n — 1 


Chapitre 3 


Représentation linéaire des Graphes 
d’Evénements Temporisés 


Introduction 


On présente sommairement les réseaux de Petri, qui sont un langage graphique per- 
mettant de décrire les systèmes à événements discrets. L'emploi des réseaux de Petri pour 
la spécification, la simulation et l’analyse de ces systèmes remonte aux années 60 et s’est 
maintenant généralisé. La littérature sur les réseaux de Petri généraux est très abondante. On 
s’est limité ici aux préliminaires indispensables sur les réseaux de Petri, le but étant d’arri- 
ver aux graphes d'événements temporisés, qui sont une sous-classe intéressante de réseaux 
de Petri pour lesquels on a des résultats algébriques plus simples et plus précis. Les graphes 
d'événements se représentent en effet par des équations récurrentes (max, +) linéaires sur les 
dates ou (min, +)-linéaires sur les compteurs. On arrive alors à une équation vectorielle de la 
forme x,:1 = À @ x, dans l’algèbre (max, +) (ou dualement, dans l’algèbre (min, +)). Cela 
permet d’attaquer les problèmes de calcul du régime périodique et du taux de production via 
la théorie des matrices dans l’algèbre (max, +), qui sera présentée au chapitre suivant. 


3.1 Introduction aux Réseaux de Petri Temporisés 


3.1.1 Réseaux de Petri 


Rappelons qu’on spécifie un graphe (orienté) par la donnée d’un ensemble S (ensemble 
des sommets) et d’une ensemble À d’éléments de S? (ensemble des arcs ou des arêtes). Un 
réseau de Petri est un graphe avec deux sortes de sommets : les places (représentées par des 
cercles) et les transitions (représentées par des rectangles). Un exemple de réseau de Petri est 
fourni sur la Figure 3.1. Si P désigne l’ensemble des places et 7 l’ensemble des transitions, 
on a donc la partition $ = P U T. Les arcs relient les transitions aux places et les places 
au transitions (mais il n’y a pas d’arc entre deux sommets de même type). Un tel graphe est 
dit biparti. Etant donné un sommet s € S, on note l'(s) l’ensemble des successeurs de s et 
T-l(s) l’ensemble des prédécesseurs de s, soient 


r(s)={sesl(ss)eA} Tl(s)={s ES] (s,s)e A}. 


A chaque arc a € A est associé un naturel n,. Des jetons circulent dans le graphe selon la 
règle suivante. On se donne un marquage M c’est-à-dire une application P — N. M(P) 
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représente le nombre de jetons dans la place P. On dit que la transition T € T est tirable ou 
franchissable si pour toute place P en amont de T a au moins n(p 7) jetons dans la place, 1.e. 


VPeT "(T), M(P)>nçr 


Plusieurs transitions peuvent donc être simultanément franchissables. Si l’on décide de tirer 
une transition 7, on retire n(p r) jetons dans chaque place amont P et l’on rajoute n(r 0) 
jetons dans chaque place aval Q. On obtient donc le nouveau marquage M’ donné par 


si PET (T), P#T(T) M'(P)=M(P)-n@pr 

si PET(T), PTT) M'(P)=M(P)+nap 

si PEeT(T)NT UT) M'(P)=M(P)+nrp-n@T 
sinon M'(P) = M(P) 


Soit par exemple le réseau de Petri de la Figure 3.1. On a représenté les naturels n, sur les 
arcs. Par exemple, np,7r, — 2. On observe que les transitions 7, et T; sont simultanément 
ürables. On a représenté ici le tir de la transition 7:. 


FIG. 3.1: Tir de la transition T> 


$ Dans la suite, en l'absence d'indication contraire sur les dessins, on supposera que les 
valuations des arcs n, valent 1, ï.e. que l’on retire ou ajoute un unique jeton pour chaque arc. 
Nous donnons maintenant divers exemples de phénomènes représentables par des réseaux de 
Petri. 

3.1.1.1 (Contrainte conjonctive). Pour le réseau de Petri de la Figure 3.2, la transition 74 
peut être tirée à condition qu’il y ait un jeton dans chacune des transitions amont. On peut 
représenter de la sorte l’assemblage de 3 pièces en proportions identiques. 


Ti D T3 Ti D T3 
Ÿ Pi \ 7 
— D CS = RE À * 
SRE &ex LU 2 
MT, LT, 
PF À 
Ne, LU 


FIG. 3.2: Conjonction 


3.1.1.2 (Portion de voie ferrée). On reprend l’exemple simplifié de voie ferrée traité dans 
l'introduction. Soit le réseau de Petri de la Figure 3.3. Un jeton présent en P; représente un 
train prêt à entrer sur le tronçon. Le jeton en place P; représente un feu vert. Lorsque la 
transition 7, est tirée, le train rentre sur le tronçon (place P;), et comme il n’y a plus de jeton 
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en P; (feu rouge), aucun nouveau train ne peut rentrer. Lorsque le train sort du tronçon (tir de 


T3), un jeton est à nouveau produit en P; (le feu repasse au vert), soit la séquence représentée 
sur la figure. 


FIG. 3.3: Feu rouge contrôlant un tronçon 


3.1.1.3 (Partage d’une ressource). On a représenté sur la Figure 3.4 deux philosophes 
dînant à l’aide d’une seule fourchette. La fourchette est initialement sur la table (un jeton 
dans la place P;). Lorsque l’un des deux philosophes décide de manger, disons le premier 
philosophe, la transition 7; est tirée, et le jeton passe dans la place P, (qui correspond à l’état 
“le philosophe 1 mange”). Pendant ce temps là 7, n’est pas tirable. Le second philosophe 


peut manger. Lorsque le premier philosophe repose la fourchette sur la table (tir de 7:), on 
est revenu à l’état initial. 


philosophe 1 mange fourchette libre philosophe 2 mange 


FIG. 3.4: Deux philosophes partageant une fourchette 


3.1.1.4 (Producteur et consommateur). Un producteur fournit une par une des pièces qu’il 
dépose dans un stock. Le consommateur retire une par une au fur et à mesure de ses besoins 
ces pièces du stock lorsque celui n’est pas vide. On se concainc que ce système est représenté 
adéquatement par le réseau de Petri de la Figure 3.5. 


producteur consommateur 
FN stock DES 
ee 1 se 
Lé Es ra FF \ É La 
| a +. 
Nr) — Lo Y 


FIG. 3.5: Producteur et consommateur 
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3.1.1.5 (Blocage). Le réseau de Petri de la Figure 3.6 est en situation de deadlock (blo- 
cage). La transition 7 ne pourra pas être tirée à cause de l’absence de jeton dans le cycle 
(comprenant la place P:). 


FIG. 3.6: Deadlock 


3.1.1.6 (Représentation formelle). La définition mathématique standard des réseaux de Pe- 
tri est la suivante [14]. Un réseau de Petri est un quadruplet (P, T, pre, post). P est l’en- 
semble des places, T est l’ensemble des transitions (disjoint de P), pre est une application 
P x T — N dite d'incidence avant post est une application P x T — N dite d'incidence 
arrière. On obtient la représentation graphique ci-dessus à partir de cette représentation for- 
melle comme suit. Si pre(P,T) > 0, on trace un arc de P à T : la place P est en amont de 
T.Si post(P,T) > 0, on trace un arc de T à P : la place P est en aval de T. On pose alors 
n(p,r) = Pre(P,T)etn(r.p, = pOSt(P,T). La condition de tir d’une transition T s’écrit 


VPEP, M(P)>pre(P,T). 
Le franchissement de T donne le nouveau marquage M : 
VPEP, M'(P)=M(P)—-pre(P,T)+post(P,T). 


Ainsi, le tir de la transition T rajoute le vecteur —pre(-, T) + post(:, T) au vecteur des 
marquages (cf. [61]). 


EXERCICE 3.1.1.7 (Philosophes Dinant). On considère maintenant 3 philosophes devisant 
autour d’une table et mangeant avec des baguettes comme indiqué sur la Figure 3.7. Il y a 
trois baguettes sur la table (une entre deux philosophes). Un philosophe a besoin de deux 
baguettes pour manger et ne peut prendre que les baguettes situées à sa gauche ou à sa droite. 
Représenter ce système par un réseau de Petri. On pourra associer à chaque philosophe i — 
1,2,3 les deux places “M; : philosophe i mange” et “R; : philosophe à réfléchit” (les deux 
activités sont exclusives). 


_philosophe 1 
8 > 
| | 


| FN, | 
\ ® © / 
philoso phe 1 S ne philosophe 2 


FIG. 3.7: 3 philosophes dînant avec des baguettes 
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3.2 Graphes d’Événements Temporisés 


3.2.1 La Sous Classe des Graphes d’Événements 


Un graphe d'événements est un réseau de Petri tel qu’en chaque place, ne rentre qu’un 
seul arc a et ne sort qu’un unique arc b, et qu’en outre n, = ny, = 1. Ces conditions inter- 
disent les configurations de la figure 3.8,(a1,2,3). On appelle conflit structurel une situation 


ceux ci ne sont pas des graphes d'événements graphes d'événements 
l T k 1 T Lee 
AS 5 Ë 
| ! Ÿ” ! | VŸ _ 
) Ç CN = D) ! C) = D) © D A .) F l Ta 
1 1 je fe) 2 ee 3 | - PE ] 
) L Ve vi. Ÿ pe _/ 
cr 
@n Se T3 | 4) A on Ts Co) 
_ | Le | 1) É eo» 13 
concurrence à concurrence à décimation des multiplicateur de 
la production la consommation jetons (1 pour 2) jetons (2 pour 1). 


FIG. 3.8: Sous classe des graphes d'événements 


où un jeton dans une place peut a priori servir pour le tir de deux transitions. Il est clair que la 
condition définissant les graphes d’événements interdit de tels conflits. Par exemple, le réseau 
de la Figure 3.8,(a2) (qui n’est pas un graphe d’événements) donne un conflit structurel. On 
peut tirer au choix 7: ou 7; pour le marquage du dessin. Si l’on tire 7>, alors on ne peut plus 
ürer 13, et réciproquement. Une telle situation représente physiquement une “concurrence à 
la consommation”. Par contre, pour le réseau (b2), l’ensemble des transitions franchissables 
pour le marquage du dessin est {71, 1, T4}. On peut tirer T1, 12, T, toutes les trois dans un 
ordre arbitraire. Cette propriété permet de définir un comportement déterministe du réseau, 
où l’on passe au nouveau marquage M’ obtenu par tir de toutes les transitions franchissables 
pour les configuration donnée. 

Une propriété centrale des graphes d'événements est la conservation du nombre de jetons 
dans les circuits. 


PROPOSITION 3.2.1.1. Pour un graphe d'événements, la somme des marquages des places 
d’un circuit donné est constante. 


Démonstration. 11 suffit de le vérifier pour le tir d’une unique transition : c’est immédiat. [ 

On dira que le graphe est vivant ssi pour tout marquage accessible M et toute transition 
T, il existe une séquence de franchissements à partir de M contenant T. Il en résulte que T 
peut être franchie une infinité de fois. Le résultat suivant est bien classique. 


PROPOSITION 3.2.1.2. Un graphe d'événements est vivant ssi il ne contient pas de circuits 
sans jetons. 


Si le graphe a un circuit sans jetons, il résulte de (3.2.1.1) que les transitions de ce circuit 
ne pourront jamais être tirées, d’où l’on conclut à la non vivacité du graphe. Nous ne prouvons 
pas pour le moment la réciproque, mais nous produirons plus bas une preuve algébrique fort 
simple dans le cas temporisé. 


3.2.2 Temporisation des Graphes d’Événements 


Dans un graphe d'événements temporisé, on associe à chaque transition T une durée 
minimale de tir 0(T). Entre la mise à feu de la transition, où les jetons sont extraits des 
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places amont et la fin du tir où les jetons sont mis dans les places aval, il s’écoule un temps 
au moins égal à 0(T). On associe également à chaque place P un temps minimal de séjour 
0(P) : les jetons doivent séjourner dans la place P au moins 8(P) unités de temps avant de 
devenir disponibles pour le tir des transitions aval. La construction de la Figure 3.9 permet 
de se ramener au cas où seules les places sont temporisées (et où le tir des transitions est 
instantané). On peut donc sans perte de généralité émettre la restriction : 


S dans toute la suite, on supposera que les places seules sont temporisées. 


transitions temporisées ñ places temporisées seules 
a TT O(T1) = 0 
| LD 
Ti F O(T)=t : CL 8(P)=t 
2 
(D) = 0 


FIG. 3.9: Réduction d’une transition temporisée à une place temporisée 


3.2.2.1 (Tâche Périodique). Le réseau de Petri de la Figure 3.10 représente une tache 
répétitive. Soit une cellule de production fonctionnant comme suit : deux palettes (initiale- 
ment en position P;) sont disponibles pour charger les pièces. Les pièces, une fois palettisées, 
sont amenées à l’une ou l’autre de deux machines travaillant en parallèle (place P:), le trai- 
tement d’une pièce prend au moins 3 unités de temps, et à l’issue, la pièce finie est déposée 
dans un stock P; (supposé de capacité infinie) et la palette libérée revient en position initiale 
P; attendre une nouvelle pièce. Comme la place P; correspond au traitement des pièces par 
l’une ou l’autre des deux machines, elle est temporisée à 3 unités deux temps (temps minimal 
de fabrication). Les deux jetons initialement dans la place P;, représentent les deux palettes. 
Comme d’après la Proposition 3.2.1.1, le nombre de jetons dans le circuit est constant, on 
note qu’au plus deux pièces peuvent être simultanément traitées. Considérons le cas où une 
grande quantité de pièces (de jetons) sont disponibles à l’instant 0 dans la place PQ et où le 
système fonctionne au plus tôt ce qui signifie que chaque transition est mise à feu dès que 
possible. La transition 7; est alors tirée deux fois instantanément. Les deux jetons restent 3 
unités de temps dans la place P;, T; est alors tirée deux fois à l’instant 3, puis les deux jetons 
séjournent un temps nul en P; et l’on est revenu à la situation initiale. On constate sur cet 
exemple que le régime au plus tôt est périodique : deux pièces sont traitées toutes les trois 
unités de temps tant qu’il reste des jetons dans la place Po, soit un taux de production de +. 
Dans ce chapitre (ainsi que dans le suivant), nous montrerons que ce phénomène est général 
et caractériserons fort simplement le taux de production. 


3.2.3 Equations aux dateurs 


Il est temps de donner les équations décrivant le comportement temporel des graphes 
d'événements. A chaque transition À, on associe le dateur x;, qui est une application crois- 
sante Z — Ra. On suppose que les jetons sont présents dans les places depuis l'instant 
—0o, et qu’à un instant initial (qu’on peut prendre égal à 0), aucune transition n’a été tirée. 
x;(n) désigne la date à laquelle arrive le n-ième tir de la transition (depuis l’instant initial). 
Le cas générique est celui de la Figure 3.11, pour lequel on a l’inégalité portée sur la figure. 
Le lecteur aura remarqué que nous notons maintenant les durées minimales de séjour dans 
les places par un nombre de bâtonnets. En l’occurrence, la place x: > x3 est temporisée à 
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Re"! fe, fo ne, fe ur 
| PA. | | | | 

Fa [SN 0(P:) =3 ONCE NE CS LES 
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K 7-57 Ne 2 Ne 27 Vel. 7 
} P: 
. : FUN 3 TS FE TN (6e 
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t=0 0<1<3 1=3 3<1<6 


FIG. 3.10: Tâche cyclique 


3 unités de temps, et la place x, > x3 à 2 unités de temps. L'intérêt de cette notation est 
de faire jouer un rôle symétrique au temps et aux quantités!. Le terme x1(7 — 1) + 2 rend 
compte d’une part de ce que tout jeton doit patienter au moins deux unités de temps dans la 
place entre x, et x; avant d’autoriser un nouveau tir de x;, et d’autre part que le jeton initial 
introduit un décalage d’une unité entre les numéros de tir des transitions x; et x3. 


XI — —— À2 


y 
Ke 


) 230) > max(2 4x (n — 1), 3 +207 — 2) 


\ 
CT Ci 
CI) GI 
L'an Lee 


= X3 | temps de séjour dans la place 
(1 unité de temps) 


FIG. 3.11: Inéquations aux dateurs pour une transition 


On considère maintenant le graphe d'événements de la Figure 3.12. Les dateurs vérifient 


a | = temps de séjour dans la pl: 
_ [Le = temps de séjour dans la place 
f ( 


qi 7 (unité de temps) 
| KA 
xl c x2 
| er) + =marquage initial 


& Le 7 | (jetons) 
on 


(in) _{(e) 
u» FE 
= 


FIG. 3.12: Exemple de graphe d'événements temporisé 


les inégalités suivantes : 


x1(n) Z  max(ui(n) +3, x2(n — 1)) 

x2(n) Z  max(l +uz(n —1), 1 +xi(n)) G.D 
x3(n) 2° max(2+x3(n — 1), x1(2), x2(n) + 1) | 
JG) 2 max(x3(n) +3,x2(n — 1)) 


lOn verra plus loin que les équations satisfaites par les compteurs sont duales de celles vérifiées par les dateurs 
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Dans le dioïde Ryax (où a @ b — max(a, b) et a @ b — a + b), on obtient le système 
d’inéquations linéaires équivalent : 


X(n) Z 3u(n)®x(n— 1) 
X(n) > lu(n—1)6 lx(n) (32) 
Gr) > 2x3(n —-1)@x1(n) & Ix(n) | 
JU) > 3x3(n) Ex —1) 
Matriciellement, cela s’écrit : 
Xx1(n) £ £ € x1(n) £ e € xi(n — 1) 
Hfn) | >| 1 &e € H(n) |@| & € € xfn—1) |@ 
x3(n) e l € x3(n) £ € 2 x3(n — 1) 
3 
u1(n) 
£ | BD | (3.3) 
£ 
x1(n) 
— 1 
vm)z[e e 3]| mo [ele e Ne .. 64 
X3(n) ‘ . 


soit un système d’inégalités de la forme : 


x(n) > Aox(n) ® Aix(n — 1) @ Bou(n) & Biu(n — 1), y(n) > Cox(n) & Cix(n — 1) . 
(3.5) 


avec des notations évidentes. L'étude du comportement au plus tôt des graphes d'événements 
temporisés est donc équivalente à l’étude des solutions minimales des équations aux dateurs 
de type (3.5), et c’est l’objet de la fin de ce chapitre ainsi que du Chapitre suivant. 


3.2.4 Equations aux compteurs 


En associant à chaque transition x; (u;, y ... ) le compteur x;(t) (nombre de tirs de la 
transition i jusqu’à l’instant t), on peut écrire le système d’inéquations (min, +)-linéaires 
suivant : 


X1() < min +x2(1), u1(1 —3)) 
x (1) < min(xi( — 1), 1+ut — 1)) y(t) < min(1 +6), x3(t — 3)) , 
x3() < min(x1(r), 206 — 1), 1 + x3(1 — 2)) 

(3.6) 


Avec les notations du dioïde Rn (où a @ b = min(a,b) eta @b = a+beta > bssi 
a @b=assi a < b), on peut écrire : 


x1() X 1x2) E u1(t — 3) 
(6) > x —-1)@ lu — 1) JO) x IMG) x; —3), (3.7) 
2x3) > x1(6) E x2(t — 1) E 1x3(4 — 2) 


On obtient en continuant de la sorte un système s’écrivant comme (3.5), mais cette fois ci 
dans le dioïde Rs et relativement aux fonctions compteurs. Les dioïdes Rain et Rax four- 
nissent donc deux systèmes d’inégalités (duaux en un certain sens) qui permettent de manière 
équivalente d’étudier le système. Dans la suite, on utilisera l’un ou l’autre de ces formalismes. 
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3.2.5 Forme Explicite “ARMA” 


De manière analogue à (3.5) et en toute généralité, les fonctions dateur au plus tôt d’un 
graphe d’événements sont solution minimale dans le dioïde Rx d’un système implicite 
d’inéquations linéaires de la forme : 


a b 
x(n) z @ Aix — à) © D Biun — j) (3.8) 
i=0 j=0 


où x(n) et u(n) sont des vecteurs d’état et d’entrée. Via le Théorème 2.2.1.1 du Chapitre 2, 
la plus petite solution de (3.8) est donnée par 


a b 
x(n) = Aix —-i)e D B;un — j) (3.9) 
i=1 j=0 
où À, = Ati Bo = AfBo, ... Noter que en toute généralité, les coefficients des 


matrices sont dans Rinax et non dans Rx (il faut être dans un dioïde complet pour que 
l’étoile converge). Aïnsi, (3.9) peut très bien donner une valeur infinie de x(n). 

L’équation (3.9) est l’analogue exact d’un système ARMA?. Bien sur, on a un résultat 
identique dans le domaine des compteurs. Calculons par exemple la forme ARMA associée 
à (3.7). Dans la pratique, on se garde bien de calculer A, mais l’on effectue des subtitutions 
dans les équations (ce qui est en fait équivalent, nous laissons le lecteur méditer cela un 
moment). Par exemple, on substitue x2(f) = x1(f—1)@ 1u2(f — 1) dans la première équation 
de (3.7), et l’on obtient après avoir substitué x,(f) dans x3(f) : 


16) = Ixi(t — 1) E 2u2(t — 1) Bui(t —3) 
HG) =xi(t—-1)@ lu —1) 
2x3) = Ixi(t — 1) E 2u2(t — 1) Puit — 3) xt — 1) E Ix3(t — 2) 


70) = 1x2) @ x3(1 —3) . 


Le lecteur notera que comme la solution minimale de l’équation ax @ b = x vérifie l’égalité, 
le comportement au plus tôt réalise les égalités par exemple dans (3.2), ce qui n’était pas 
évident a priori. 

Nous sommes maintenant en mesure de caractériser les graphes d'événements tempo- 
risés non bloqués. Nous dirons que le graphe est vivant si, en présence d’entrées infinies (1.e. 
u — & non contraignante), chaque dateur au plus tôt x(n) vérifie x(n) < +oo (tout événement 
arrive en un temps fini, ou encore, chaque transition est tirée un nombre infini de fois). On a 
alors : 


THÉORÈME 3.2.5.1. Le graphe d'événements temporisé est vivant si et seulement il n'y a 
pas de circuit sans jeton et de temporisation strictement positive? 


Démonstration. On part de l’équation 3.8 dans le domaine des dateurs. Pour que x(n) < 
+00, il suffit que l’étoile de la matrice À, ait ses coefficients inférieurs à +co, donc, d’après 
la Proposition duale de 2,2.2.5.1, que Ao n’admette pas de circuit de poids positif, ce qui 
donne la condition de l’énoncé. La réciproque est claire, car en présence d’un circuit sans 
jeton de poids c > 0, on peut écrire pour un dateur x; d’une transition de ce circuit une 


inégalité de la forme x;(n) > x;(n) + c ce qui montre que x;(n) ne peut être fini. [] 
2j.e., auto-régressif à moyenne mobile 
FLa temporisation d’un circuit est égale à la somme des temporisations de ses places 
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Il y a ici une petite difficulté si l’on considère un circuit sans jetons et de temporisation 
nulle comme sur la figure 3.13, (a). L’inéquation aux dateurs s’écrit x(n) > x(n)Œu(n) 
et admet donc la solution minimale x = u : ce système se comporte comme l’iden- 
tité. Cependant, lorsque l’on considère ce graphe comme un graphe d'événements non 
temporisé, on est dans une situation de deadlock (cf. 3.1.1.5). On retrouve cependant la 
situation de blocage si l’on temporise positivement la place comme en (b). On obtient 
alors l’inéquation x(n) > max(1 + x(n), u(n)) qui, dès que u(n) + &, admet l’unique 
solution triviale x(n) — + (1.e. la date de l’événement n est +co : blocage). On 
retiendra que le modèle algèbrique des graphes d'événements temporisés est cohérent 
avec la théorie des graphes d'événements non temporisés sous jacents à ce genre de 
pathologies près qui disparaissent si l’on suppose par exemple que toutes les places ont 
des temporisations strictement positives. 


an sd 
temporisation nulle {7} 1. 
HESENTE À qd xx place temporisée à une unité 


\ ET 
er 4 J 


ee | KE 
(a) non bloqué (b) bloqué 


FIG. 3.13: Cas pathologique 


3.2.6 Forme d’état 


Par des manipulations combinatoires tout à fait classiques (en augmentant l’état), on peut 
ramener un système ARMA de type (3.9) à la forme dite d’état 


x(t) = Ax(t—1)@ But), y()=Cx(). (3.10) 


En l'occurrence, en prenant les vecteurs x (1) = [x1(r), x2(f), x3(t), x3(1—1), x3(1—2), x3(1— 
3)]! etu(r) = [uit — 3), u2(t — 1)], on se ramène à la forme standard 


x(t) = Ax(ft—1)@ But), y(t) = Cx(t) (3.11) 
avec 
1 € & €& € € e 2 
6e £ E£ £ € € | 
À = a De ; Cie. ee ei 
£ £ € £ € € £ € 
£ £ £& e E € £ € 
£ £ £& €£ e € £ € 


(3.12) 


Un cas particulier intéressant est celui d’une entrée non contraignante (4 = €). L’équation 
(3.10) se réduit alors à x(1) — A'x(0), de sorte que l’étude du comportement du graphe en 
régime autonome se ramène au l’étude des puissances de la matrice À. On a déjà noté sur 
l’exemple 3.2.2.1 que certains graphes d’événements étaient périodiques, i.e. que pour un 
certain taux de production À, on avait une relation de la forme 


Xr = TXÀ+X() 
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soit dans les notations du dioïde 
T 
Xe4r = À X(1) 


(toutes les T unités de temps, À x T tirs supplémentaires arrivent pour chaque transition). 
Au chapitre suivant, on va montrer 1/ que sous une hypothèse de forte connexité du graphe, 
le système est périodique à partir d’un certain temps, 2/ que le taux de production est égal à 
l’unique valeur propre de la matrice À, 3/ que cette valeur propre est simplement caractérisé 
comme le poids moyen maximal des circuits du graphe (et se calcule donc aisément). 


EXERCICE 3.2.6.1. 1/ (Facile). Ecrire l’équation d’état du graphe d'événements temporisé 
de la Figure 3.14. 2/ (Plus dur). Peut on trouver des matrices À’, B’, C’ de dimensions plus 
petites telle que la relation entrée sortie 4 > y soit identique ? 


RES 7 sn. 
x] re AS , GAONsAAOnE en toutes lettres 
= ans les places) 


FIG. 3.14: Quelles sont les équations de ce graphe d’événements ? 


3.3 Représentation entrée-sortie 
3.3.1 Représentation par des sup-convolutions 
Soient À e R!*", B € Rnat , C € Rhax . Considérons le système à p entrées et q sorties : 
x(k) = Ax(k — 1) Bu(k), y(k) = Cx(k) (3.13) 


dans l’algèbre (max, +). Il vient, en remplaçant x (k — 1) par sa valeur dans (3.13) et à l’issue 
d’une récurrence immédiate : 


1-1 


VIEN y&)=CAxk-1D@ Œ CA'Bu(k—i) (3.14) 
i=0 
Posons 
€ sik < 0 
h(k) = 
@) CAB sinon. 


Pour simplifier, nous supposerons dans un premier temps À scalaire (on considère un système 


mono-entrée mono-sortie, soient C € Rl*" et B e R’*!). 


PROPOSITION 3.3.1.1. La plus petite solution du système (3.13) est donnée par : 


y(n) = @A(p) & u(n — p) = supla(p) + u(n — p)] . (3.15) 
peZ pe 


50 3. Représentation linéaire des Graphes d’Événements Temporisés 


Démonstration. Au vu de (3.14), il est clair tout y vérifiant (3.13) est plus grand que la 
somme à droite de (3.15). Réciproquement, il faut montrer que cette somme est effectivement 
solution. Posons 


(ee) 
x) = @ A'Buk—i), FU) = CYR) = hxu(R) . 
i=0 
Il est immédiat que (x, y) est solution de (3.13) ce qui conclut la preuve. (| 


EXERCICE 3.3.1.2. Que devient la proposition pour un système à plusieurs entrées et plu- 
sieurs sorties ? 


L'expression à droite de (3.15) s’appelle sup-convolution de h par u. Cet dénomination 
est justifiée par l’analogie avec le produit usuel de convolution 


+00 
J x V0) = JATIEG = T)dr. 
On notera dans la suite 
UxV: uxV(n) = Suplu(n — p) + v(p)]. (3.16) 
pEZ 


On vérifie aisément les propriétés suivantes : 
PROPRIÉTÉS 3.3.1.3. 1. uxv=vxu 
2. ux(v®Ow)=uxv@uxw 
3. (u x v) x W = u *x (u x W) 


Soit d’autre part : 


On a 
exU—=uUkxe—=U 


pour tout u € R?,.. Autrement dit, e est l’unité pour le produit de sup-convolution. e est 


donc l’analogue de la fonction & de Dirac, et pour cela, on appellera e impulsion. On notera 


| PR TZ ; : 
que si u, v € RZ,,, on a en général u x v € R,, (il peut se faire que le sup dans (3.16) soit 


infini, i.e. u + v(n) = +00 € Ra). On résume ces observations par la proposition suivante. 


—Z 
PROPOSITION 3.3.1.4. KR, muni du max (point par point) et du produit de sup- 


convolution est un dioïde, ayant la fonction nulle k > & pour zéro et e pour unité. 


Il résulte de (3.15) et de h x e — h que la sortie correspondant à l’impulsion e est égale 
à h. Autrement dit, h s’interprète comme la réponse impulsionnelle du système. 

Plus généralement, on pourrait se demander quels sont les systèmes représentables par 
des sup-convolutions. Dans l’algèbre usuelle, on sait que les systèmes linéaires stationnaires 


se représentent par des convolutions (en un sens plus ou moins fort suivant la régularité). 
Ici, on a un résultat analogue. Soit S un système (max, +) linéaire sur R,.,, i.e. un système 


vérifiant 


max? 


S(u @ v) = S(u)BS(v), S(AU) = ÀS(u) 
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pour tous signaux w,v et pour tout scalaire À. Notons y l’opérateur de décalage en 
numérotation : 


yu(n) =u(n—T). 
Un système est dit stationnaire s’il commute avec l’opérateur de décalage, 1.e. 
Sy =yS . (3.17) 


Autrement dit, si (n+> u(n),n + y(n)) est un couple entrée-sortie admissible, alors, (n + 
u(n—1),n + y(n — 1)) l’est également. Par une induction immédiate et comme y est 
bijectif, la propriété (3.17) entraîne 


VkeZ, Sy =y"s. 


(si on décale l’entrée de k numéros, la sortie est décalée de même). Nous aurons besoin d’une 
hypothèse de régularité. Le système S est dit continu s’il préserve les sup infinis, 1.e. 


(Bu) = PS) 


iel iel 


; . : —Z 
pour toute famille de signaux {u;};er avec Vi el, u; e R.. 


rs : s ne 2 | | 2 3 
THÉORÈME 3.3.1.5. Soit S un système linéaire stationnaire continu sur R,... Îl existe une 


—Z 
application h € R,.., (dite réponse impulsionnelle) telle que 


max 


Vu, Su=hxu. 
Démonstration. On a pour toute entrée u, u = u *xe, soit 


u(n) = Buke(n —h) = Buky'e(n) 


keZ, keZ 


soit 


u = ŒBuky'e : 


k 


Par continuité, linéarité et stationnarité de S : 


Su = Buks(y'e) — MBuky'sSe ; 


keZ, keZ, 


1.e. 


(Su)(n) = Buk)(Se)(n — k) = (Se) + u)(n) . 


keZ, 
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3.3.2 Série de transfert 


Il y a une manière purement algébrique de retrouver le résultat de représentation par des 
sup-convolutions. Le système d’équations (3.13) s’écrit en utilisant l’opérateur y de décalage 
sur les quantités. 


x = Ayx @ Bu, y—=Cx. 


En appliquant le Théorème 2.2.1.1 du Chapitre 2 sur l’équation implicite x — ax @ b, on 
obtient la plus petite solution 


x = (47) Bu = Œ A‘By'u 
kEeN 
y =C(Ay) Bu . 


La série formelle 


H = C(Ay)B = ŒCA'By" ER? [II] 
kEeN 


est appelée série de transfert. On a 


y) = Hu(n) = @ CA‘ By u(n) = CA Bu(n —k) 
keN keN 


de sorte que l’on a retrouvé la représentation (3.15). 
EXEMPLE 3.3.2.1. Considérons le système d’ordre 1 : 
x(n)=3x(n—1)@lu(n), y(n) = 1x(n) . 


Onax = 3yx@ lu, y = 1x, d’où 


y =2(G37y)'u = P2r3u ; 
ieN 


La réponse impulsionnelle vaut donc 


£ sii <0 
h{(i) — | 3.18 
@ a sinon. ( ) 


3.3.3 Dioïde des dateurs 


Dans le cas des graphes d'événements, les dateurs sont des fonctions croissantes. Nous 
noterons Dat l’ensemble des applications croissantes Z — R:,,. Posons 


£ Sik <0 
5 e'(k) = 3.19 
HA [ ik > 0 Hi 


e' joue un rôle analogue au “Dirac” e, mais relativement aux applications croissantes : 


NET —Z : . 
PROPOSITION 3.3.3.1. Une application u € KR, est croissante ssi u x e' = u. 
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Démonstration. 


u xe'(k) = suplu(K — p) +e’(p)] = supu(k— p) = supu(l) 
p p>0 I<k 


d’où la conclusion. (E] 


COROLLAIRE 3.3.3.2. L'ensemble Dat des applications croissantes Z — Rx muni du 
max et du produit de sup-convolution est un dioïde, ayant e’ pour élément unité. 


Si un transfert h est croissant, alors on à h — h x e’, ce qui montre que h est la réponse 

à l’entrée e’. L’entrée e’ peut être qualifiée “d’impulsion croissante”. Elle s’interpète physi- 

quement comme suit. Au vu de (3.19), e’ signifie que les les événements numérotés 0,1,2,3 

. arrivent à l’instant 0 (et que les événements antérieurs se sont produits à l’instant —co). 
Autrement dit, le système commence à fonctionner (en régime maximal) à l’instant 0. 


3.3.4 Opérateurs sur les dateurs 


5 CARS 2 Sa es 
L'opérateur y est a priori défini sur R,,,,. On peut cependant considérer la restriction de 
y aux dateurs (1.e. aux applications croissantes). 


PROPOSITION 3.3.4.1. Pour toute application croissante u, on a 
QE hu = y Pa. 
Démonstration. On a 


(y"@ y )u(k) = max(u(k — n),u(k — p)) = u(k— min(n, p)) = y" Pu(k) | 


PROPOSITION 3.3.4.2. On a plus généralement la règle de simplification : 
t>setn<p —= VueDat, (ty"@sy?)u =ty'u . (3.20) 
Démonstration. 


@y"@sy?)u = (G(y"@7y?)@sy?)u (via 3.3.4.1) 
= (y @(6y?sy?)u 
— (ty"@ty')u (cart ®s =t) 
= ty" 
Cette règle admet une interprétation simple en termes graphiques. Par exemple, les graphes 
d’événements temporisés des figures 3.15,(a1) et (a2) sont équivalents du point de vue de la 


relation entrée sortie, 1.e. réalisent la même application au plus tôt u > y. De même pour les 
graphes (b1) et (b2) (on rappelle que les temps de séjour sont représentés par des bâtonnets. 


EXEMPLE 3.3.4.3. Soient deux machines travaillant en parallèle et traitant les pièces à rai- 
son d’une toutes les trois unités de temps. On a représenté sur la Figure 3.16,(a) le graphe 
d'événements temporisé, en (b) sa transcription opératorielle et les équations correspon- 
dantes. Il vient en substituant x: : 


22 > 3y/X%2 D 3, Y > X 


d’où la solution au plus tôt 


54 3. Représentation linéaire des Graphes d’Événements Temporisés 


Æ, u L u u 
RSS (.) () (D = ab 
Re ha 
(al) (2) (b1) (b2) 


FIG. 3.15: Interprétation graphique des règles de simplification 


dates 
—ù # — 76) 
ne [ $ \ 
( ) Ke) x > 40 y2x : 
2 x f = x El | 
|: + ( EN (3 x pa — 
DID) = v)U) 1 
| x2 > 3x1 + 
= » re #2 T—< 
‘a @) 1) 
Lo) V7 Y XX u( A RP NER RCE ARES événements 
LL y k — } \ 
! — 09 1 8" 
TL 75 7 2: k u(&) = e/(k) 
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FIG. 3.16: Graphe d’événements, opérateurs associés, équations, réponse impulsionnelle 


y = 3(G7y°)'u = Are | (3.21) 
keN 
Finalement 
y) = hxu(k) = sup[3 x ( +1) +u(k —2i)] . (3.22) 
ieN 


En prenant pour entrée l’impulsion croissante e’ (une infinité de jetons arrivent en w à l’instant 
0), on obtient pour sortie une fonction croissante en escalier à accroissements périodiques de 


3 
pente moyenne 3. 


EXEMPLE 3.3.4.4. Sur la planète Mars, il y a deux mathématiciens qui collaborent de la 
manière suivante : le premier produit trois conjectures simultanées tous les trois ans à partir 
de l’année zéro. Le second travaille exclusivement sur les conjectures du premier, met deux 
ans pour élucider une conjecture, et peut travailler sur au plus deux conjectures en même 
temps. On prétend qu'après un transitoire, une conjecture sera élucidée chaque année. Ce 
résultat peut sembler non intuitif, dans la mesure où l’on s’attend plutôt à un phénomène de 
ppem des cyclicités. La preuve algébrique qui suit montre que les cyclicités sont en fait reliés 
à certaines équations linéaires en nombre entiers, pour lesquelles via Bezout, on obtient le 
pgcd des cyclicités. 

On se convainc que le système représente par le graphe d'événements de la Figure 3.17 
avec l’entrée u — e’. Le premier mathématicien est représenté par la transition x. En fait, on 
a x = 3y°x De’, ce qui entraîne x1 = (3y°)*e’. On voit facilement que* 


351 sik>0 


X1(k) = 
1) —0o sik <0, 


4[x] note la partie entière de x 
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ce qui représente bien trois événements se produisant toutes les trois unités de temps à partir 
de l’instant 0. x3(n) représente la date à laquelle la n-ième conjecture est résolue par le second 
mathématicien. En éliminant l’inconnue intermédiaire x2, on obtient le dateur au plus tôt : 


xX1 > 3y9x] @u dates 


N\ 
( Ka TK AR Hs 
_ 


ee Nec tue à 7 22 > y?x3 Ex 
LS EX a 
| K-y | EE | TT | X3 > 2x2 
e X] X2 X3 + 


T 
événements 


FIG. 3.17: Deux producteurs placés en série 


x3 = 2027) Gy)"e . 
Or 
(37°) y?) = (æc x D) (æe x 1) = AT xi+2x jy ri, 
ieN jeN Lien 
On a 
0=3x0+2%x0, 
2=3x0+2%xI, 


3—-3xX1+2X0 
4—=3x0+2x2 


et l’on voit sans difficulté que tout entier supérieur à 2 s’écrit sous la forme 3i + 2j. On a 
donc 


37°)" Cy?)" =eB2y 37 D4y*@...—eæ2y (1y)*. 
Finalement, les dates de résolution des conjectures sont données par 
x3 = 2(eB2y(2y)")e' . 


Il s’agit de la fonction en escalier représentée à droit ede la Figure 3.17, à accroissements 
ultimement périodiques de période 1 (alors que les deux mathématiciens ont des périodicités 
de 3 et 2). C’est donc un cas où des pgcd interviennent pour le calcul de la période. 


3.3.5 Opérateurs de décalage sur les compteurs 


Pour les compteurs, la situation est complètement duale. Un compteur est une application 
croissante : Z — Run. Notons 6 l’opérateur de décalage : 


Su(t) ut —1). 


ô vérifie la règle suivante duale de y" @ y? = y"int@p) : 
(Do) NE Pur. (3.23) 
Cela résulte en effet de 
(8° D6')u(t) = min(u(t —t}),u(t —s)) = u(t — max(s,f)) . 


Nous laissons le lecteur se convaincre que toutes les manipulations précédentes en dateurs 
admettent des versions duales en compteurs. 
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3.3.6 Application : simplification de Graphe d’Événements 


Nous concluons ce chapitre en calculant la série de transfert du graphe d'événements 
de la Figure 3.12. Il s’agit d’un exemple spectaculaire où l’application de règles algébriques 
permet de simplifier les graphes d'événements. Nous pensons que le lecteur a intérêt à re- 
garder en détail la résolution (apparemment technique) des équations qui suivent. Une bonne 
partie des propriétés de de l’opérateur y y apparaît. On pourra aussi se référer à [27] pour la 
théorie générale. On a d’après (3.2) les équations au plus tôt relatives aux dateurs : 


X(n) = 3ui(n) Bxi(n — 1) 
x3(n) =  2x3(n — 1) Ex1(n) & lx(n) 
JG) =  3x3(n) ®x2(n — 1) 
Soit avec l’opérateur y. 
X1 = 3u1 D 7x2 
%X —= lyw@lx 
3.25 
3 —= 2x3 x © Lx A 
y = 3x6 yx 
En substituant x, dans la seconde équation, il vient 
%2 = lyu2 @ dus @ 1yx . 
La solution au plus tôt vérifie donc : 
%2= (y) (yu2 du) 
d'où x — 3u1 ®7y(1y)(lyu2 @ du) 
En notant que 
3®4y(y) =3e@æ1y(y))=3(1y) , 
on obtient 
x = 3(1y)'u @ 17*(1y)'u | 
En reportant les valeurs de x; et x: dans l’équation en x: : 
23 = 2yx3 © 3(1y)'u & 17 (y)'u  1(17) (yu2 D dur) : (3.26) 
Par application de (3.20), on peut écrire 
347) 57) = 5(y) 
ly Gr) @2y(y) = @y@17)y) =2y7(yY 
L’équation (3.26) se simplifie donc en 
x3 = 2yx3 @ 5(1y)*u1 B2y(ly)'u . (3.27) 


En prenant la plus petite solution de (3.27), on obtient 
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3 = (27) [27 (y)'u2 ES(y)'u] . (3.28) 
Comme dans un dioïde commutatif, (a @ b)* — a*b* (cf. Chapitre 2), 


27) (y) =(2@1)7) = 27) 


d’où 
23 = 5(2y)'ui 27 (2y)'u . 
Ainsi, 
y = 3x3 @ yx2 = 8(2y)'u1 57 (27)'u2 y (y) (yu2 P du) . (3.29) 
On a 


YAy)' 17 <5@y)'y, 47(y) < 8027) 
Par exemple, la première inégalité résulte de 17? < 5y et de 


Qyy = Bey < Gr) = MC x by! 


keN keN 


Argument analogue pour la seconde inégalité. Finalement, 
* ES * * u 
y = 8(27)'ui 527)" yu2 =[ 827) 5C27)*y 1| : | 


H = [ 8(2y)* 5(27)'y | mérite le nom de matrice de transfert du système. Le lecteur se 
convaincra que le graphe d'événements à droite de la Figure 3.18 admet le même transfert 
H que le graphe d'événements à gauche de la Figure (par exemple, le circuit avec un jeton 
et deux bâtonnets est responsable de l’étoile (27 )*, et plus généralement, chaque chemin du 
graphe à droite correspond à un terme du transfert). On constate que certains arcs du réseaux 
de Petri précédent étaient en fait redondants et ont été supprimés. 


1 2 
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. 
(1 © 1 2, 
L x2 
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FIG. 3.18: Graphe d'événements simplifié 
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Notes et références bibliographiques 


Nous avons laissé de coté nombre de questions intéressantes relatives aux réseaux de Pe- 
tri généraux. Disons que l’un des problèmes essentiels liés à ces réseaux consiste en l’étude de 
l’ensemble des marquages accessibles (étant donné un marquage M, existe t-il une séquence 
de franchissements qui y conduise ?). Une bonne référence générale est Brams [14]. L’exposé 
de Reutenauer [61], plus court et d’inspiration plus mathématique nous paraît hautement re- 
commandable. Citons également Reisig [60] et Peterson [54]. Le lecteur pourra comparer 
avec profit le formalisme des réseaux de Petri avec celui du Grafcet décrit par Blanchard 
[11]. Pour les graphes d’événements temporisés, on pourra se reporter à la thèse de Chre- 
tienne [22], à Carlier et Chretienne [18], au livre de Baccelli, Cohen, Olsder et Quadrat [4] et 
aux articles [24, 25, 26, 27]. 


Chapitre 4 


Propriétés spectrales et 
asymptotiques des matrices (max,+) 


Introduction 


On a vu au chapitre 3 que certains systèmes à événements discrets en régime autonome 
et en particulier les graphes d'événements temporisés lorsque les entrées ne sont pas contrai- 
gnantes, donnaient lieu à des équations de type x, — Ax,_; dans les algèbres (max, +) 
(équation aux dateurs) ou (min, +) (équation aux compteurs). Cela motive l’étude des pro- 
priétés de la suite {A"},en. Le résultat essentiel de ce chapitre est que, sous une condition 
d’irréductibilité, on a une propriété de cyclicité de la forme A"*° = )°A" pour n assez 
grand. Dans l’algèbre (min, +), À s’interprète comme le taux de production du système. En 
fait, À est caractérisé comme l’unique valeur propre de la matrice A. Ce chapitre présente 
d’abord les propriétés spectrales des matrices (max, +) ou (min, +) (valeurs propres et vec- 
teurs propres). On montre ensuite le théorème de cyclicité. On s’est efforcé de mettre en 
évidence le parallèle avec l’algèbre usuelle, qui rend tous ces résultats très naturels. 


4.1 Valeurs propres, vecteurs propres 


On pourrait présenter la théorie spectrale des matrices (max, +) de manière autonome. 
Cependant, nous pensons qu’elle se comprend mieux par analogie avec un théorème classique 
relatif aux matrices à coefficients positifs ou nuls : le théorème de Perron-Frobenius. Ce 
théorème s’énonce à l’aide de la notion clé de matrice irréductible, que nous présentons 
maintenant. Cette notion qui remonte au moins à Markov et à Frobenius, n’est autre qu’une 


version matriciellle de la forte connexité des graphes. 


4.1.1 Matrices irréductibles 


Dans toute cette sous-section, (S, @, @) désignera un semianneau (non nécessairement 
idempotent), dont le zéro sera noté €, vérifiant les deux propriétés suivantes 


aBbhb=s—=a—set b=e, (4.1) 
a@b=s—a—s ou b=e. (4.2) 
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Nous appelerons positifs! ces semianneaux. 


On rappelle qu’à une n x n-matrice À à coefficients dans un semianneau &, est associé 
le graphe (orienté) à n sommets 1,... ,n et avec un arc du sommet j au sommet ; dès que 
Ai; Æ € 

Dans un graphe orienté, on définit la relation de forte connexité : iRj s’il existe un 
chemin de i à j et aussi un chemin de j à i. Par convention, il existe toujours un chemin (de 
longueur nulle) de à à i. Les classes d’équivalence sont appelées les composantes fortement 
connexes. Un graphe est dit fortement connexe s’il y a une seule classe d'équivalence pour la 
relation R. 


PROPOSITION 4.1.1.1. Soit S un semianneau positif. Soit À € (S*)"*". Les propositions 
suivantes sont équivalentes. 


1. Le graphe associé à la matrice A est fortement connexe. 


2. On ne peut pas partitionner la matrice À sous la forme triangulaire par blocs : 


a=r[? PE (4.3) 


où P est une matrice de permufation. 


3. (1 A)" ! a tous ses coefficients strictement positifs. 
La matrice À sera dite irréductible lorsqu'elle vérifie les conditions ci-dessus. 


Démonstration. Étant donné un naturel non nul n eta e S , nous notons? n.a =a®:.:@a 
(n fois). Lorsque a,b € S commutent, on a la formule du binôme (qui se prouve comme 
d'habitude) : 


@a@b) = ŒCrab"*. (4.4) 
k=—0 


Moyennant cette identité, (iii) (i) résulte de ce que (71 & AY = &, (OR LE En effet, 
cette somme est non nulle ssi l’un des A; est non nul, i.e. compte tenu de l’interprétation 
combinatoire du produit de matrice produite au Chapitre 2,$2.2.4, si et seulement s’il existe 
un chemin de j à à dans le graphe associé à A. 

(i)=() : Après renumérotation des vecteurs de base, on peut supposer P = Id. On constate 
alors qu’il n’y a pas de chemin allant* “du bloc B au bloc D”, ce qui contredit la forte 
connexité du graphe associé à A. La réciproque est aussi facile, et nous laissons le lecteur 
s’en convaincre. En fait, (ii) et (i) ne sont que la même propriété énoncée dans deux langages : 
l’un matriciel et l’autre combinatoire. Ü 


ÎLes semianneaux vérifiant la propriété (4.1) sont appelés zerosumfree par Golan [41]. Les seuls semianneaux 
positifs que nous aurons à traiter ici sont le semianneau des réels positifs ou nuls (RŸ, +, x), le semianneau de 
Boole (B, @, ®) (cf. Chap. 1,1.2.4.2), le semianneau (max,+) (1.e. Rimax) et ses dérivés. Le lecteur pourra s’il le 
souhaite remplacer S par Rmax Ou R+ dans toute la suite. 

2Par exemple, n.a = n x a dans (R*,+, x), n.a = a dans Rmax. Il faudrait dire qu’il y a un unique 
morphisme de monoïde o : (N, +) — (S,@), tel que w(1) = a, que l’on note n.a = (a), et que l’opération 
(n,a) — n.a équipe S d’une structure de (N, +, x) semimodule, qui devient triviale dans le cas des semianneaux 
idempotents, mais ce serait très pédant. 

3a priori de longueur au plus n — 1, mais dans un graphe à n sommets, il existe un chemin p’ de j à i ssi il 
existe un chemin p de longueur au plus n — 1 de j à i. Un tel chemin p s’obtient par exemple en éliminant les 
circuits éventuels du chemin p’. 

4de manière précise, si B est de taille p et D de taille q (et donc p + q = n), pour 1 <i < petp+1< j< 
p + q,iln’y a pas de chemin allant du sommet i au sommet j. 
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L'intérêt de l’irréductibilité tient à ce qu’un vecteur propre d’une matrice irréductible 
a toutes ses composantes non nulles. L’on peut qualifier cette propriété de syndicale (un 
implique tous). 
LEMME 4.1.1.2 (SOLIDARITÉ DES MATRICES IRRÉDUCIBLES). Soit S un semianneau 
positif, À € S"*" irréductible, et u € (R*})" \ {e} tels que Au = ru. Alors u; Æ € pour 
touti e{l,...,n}. 


Démonstration. Soit j un indice tel que u; Æ &. De ru; = (Au); = @, Axux, on tire que 
tout à accessible en un coup depuis j (1.e. tel que À;; €) vérifie u; 5 €. En raisonnant 
de proche en proche, tout à étant accessible en un nombre fini de coups depuis j, on a que 
ui Æ €. [] 


EXERCICE 4.1.1.3. Montrer qu’une matrice carrée À à coefficients dans un semianneau 
positif est nilpotente ssi son graphe n’a pas de circuits. Caractériser l’indice de nilpotence 
(.e. le plus petit & tel que AÏ = £) en termes de graphe. Est il vrai que le graphe d’une 
matrice nilpotente soit un arbre ? une forêt ? 


EXERCICE 4.1.1.4. Montrer qu’un semianneau S est positif ssi l’application x : S$ — 
(où B = {6, e} désigne le semianneau des booléiens), définie par x (e) = & et (x) = e pour 
x ZÆ &, est un morphisme de semianneaux. Application : soient A;,..., 4, des matrices 
carrées à coefficients dans $. Donner un algorithme pour décider s’il existe un entier r et 
un choix d’indices ÿ,,... ,i, tel que le produit À; ...A;, soit nul. Donner une borne (ne 
dépendant que de la dimension ambiante n et du nombre de générateurs k) pour la plus petite 
valeur de r. 


4.1.2 Un peu de culture : théorie classique de Perron-Frobenius 


Cette section a d’abord une fonction culturelle et peut être sautée par le lecteur qui ne 
s'intéresse qu'aux systèmes à événements discrets. 

La théorie spectrale des matrices max-plus, que l’on va présenter, est motivée par la 
nostalgie du résulat suivant, prouvé par Perron en 1907 pour les matrices à coefficients tous 
positifs et par Frobenius en 1912, pour les matrices irréductibles à coefficients positifs ou 
nuls. 


THÉORÈME 4.1.2.1 (PERRON-FROBENIUS). Soit À € (R'})'*" irréductible. Le rayon 
spectra® de la matrice À est valeur propre de A. L'espace propre associé est une droite 
engendrée par un vecteur u à composantes toutes strictement positives. 


Avant de montrer ce théorème, on va en proposer une justification heuristique partielle. 
Considérons la matrice 
1 2 
À = 
2 1 


et introduisons le cône positif C = (R*}°. Soit f l’endomorphisme associé à la matrice A. 
La positivité des coefficients de À entraîne que C est envoyé sur un nouveau cône 


a=ro=r|,|+r|i]ec 


De même, C; est envoyé par f sur un cône C2 € C1, et ainsi de suite. En peut donc définir 
Co = (,1 En. On a représenté ces cônes sur la Figure 4.1. On constate que C, se rétrécit et 


Sie. le plus grand module des valeurs propres de À, soit o(A) = sup{|A| | À vap de A}. 
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converge vers la demi-droite D formée des multiples (positifs) du vecteur [1, 1], demi-droite 
qui est forcément propre (car f(D) € D par passage à la limite dans f(C,) € C,), ce qui 
montre que À admet au moins une valeur propre positive associée à un vecteur propre positif. 
Une démonstration de ce résultat consiste en suivant cette idée à se placer dans l’espace 


(ee) 


CG 
C 


X1 


FIG. 4.1: Interprétation géométrique du théorème de Perron-Frobenius 


projectif® et à montrer que l’application induite par f admet un point fixe (c’est à dire une 
droite invariante, donc un vecteur propre à composantes positives), ce qui peut se faire en 
appliquant le théorème du point fixe de Brouwer (une application continue d’une boule de R' 
dans elle même admet un point fixe), on montre alors grace à l’irréductibilité que l’application 
linéaire projective associée à f envoie dans lui même un compact convenable homéomorphe 
à une boule. On peut aussi avec Birkhoff montrer que f est contractante pour une métrique 
projective [10]. Maintenant que le lecteur devrait s’être convaincu du résultat, on s’autorise 
à parachuter la preuve suivante (due à Wielandt) a priori moins naturelle mais beaucoup plus 
précise. 


Démonstration. On introduit la “sous-valeur propre maximale” : 
r = sup{t e R* | Au e (R*}'\ {0}, Au > tu} . 


Si u est un vecteur propre associé à la valeur propre À, en notant |#| le vecteur dont les 
coordonnées sont les valeurs absolues des coordonnées de 4, on a par l’inégalité triangulaire : 


Qu = Allu| = |Aul < Alul 


d’où il résulte que r > |À|, et en prenant le sup sur les valeurs propres, r > p(A). Soit une 
suite (ux,r£) telle que uy € (R*})", Aux > ru, rx — r. Quitte à normaliser wy, on peut 
supposer que u4 appartient à la sphère unité S. Par compacité de S N(R* }”, la suite uy admet 
un point d’accumulation u € S N (R* }" qui vérifie 


Au > ru . 


Supposons que (Au); > ru;. Comme (7 + A)" a toutes ses composantes positives, en posant 
v = (1 + A)'u, on obtient : 


Av = (1 + A)" Au > r(I + A)'u =rv 


avec maintenant inégalité stricte pour chaque composante de v, ce qui contredit la maximalité 
de r. On a donc Au = ru etr — p(A). Nous montrons maintenant que l’espace propre 
associé est de dimension 1. Soit w vecteur propre (non nécessairement réel) associé à o (A). 
De Aw = rw, on tire 


6j. le quotient de (R*}"* \ {0} par la relation d’équivalence : x < y ssi 34 scalaire, x = y 
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Alw| 2 rlwl, 


et en raisonnant comme ci-dessus, A[w| = r|w|. Il résulte du Lemme 4.1.1.2 et de la re- 
marque qui précède qu’un vecteur propre pour r a toutes ses composantes non nulles. Soient 
x et y deux vecteurs propres associés à la vap r. Le vecteur z = x, — y,x à sa première co- 
ordonnée nulle et vérifie Az = rz, donc z = 0, ce qui montre que x et y sont proportionnels. 
Donc l’espace propre associé à r est la droite vectorielle de direction u. Ü 


Nous allons voir maintenant de quelle manière ces résultats s’étendent au cas (max, +). 


4.1.3 Valeurs propres et vecteurs propres (max, +) 


On étudie les valeurs propres et vecteurs propres de la matrice À, 1.e. les vecteurs u € 
Rx \{e} et scalaires À € RAA tels que 


Au = Àu . (4.5) 


4.1.3.1 (Interprétation physique des vecteurs propres). Soit un système de production en 
régime autonome vérifiant l'équation aux dateurs 


x(n +1) = Ax(n) . 


Si u est vecteur propre de À, et si l’on prend la condition initiale x, = u, on obtient x(n+1) — 
À @ x(n) — À + x(n). Autrement dit, une pièce est produite toutes les À unités de temps : À 
est donc l’inverse du taux de production. En outre, le vecteur propre # donne une condition 
initiale autorisant un fonctionnement périodique (de période 1) du système. 


De manière analogue au cas classique (cf. Proposition 4.1.1.1) on dira qu’une matrice est 
irréductible si le graphe associé est fortement connexe, ou, ce qui est équivalent, si la matrice 
(1 & A}"-! a tous ses coefficients > &. 


EXEMPLE 4.1.3.2. Voir Figure 4.2. On a représenté le coefficient a;; sur l’arc j + à. La 
matrice est irréductible (il y a un chemin entre n’importe quels sommets). Par contre, la 
matrice 


= ? 
1 [ = + ROUES 
KE) e | 


FIG. 4.2: Une matrice et son graphe 


THÉORÈME 4.1.3.3. Une matrice À € (Ruax)'*" irréductible admet une unique valeur 
propre, notée p(A), donnée par 


n 


p(A) = ee) ANR, (4.6) 


k=1 


64 4. Propriétés spectrales et asymptotiques des matrices (max,+) 


Avant de montrer ce résultat, nous faisons quelques remarques sur la quantité p(A) 
définie en (4.6). On note que dans l’algèbre usuelle, disons pour une matrice À dont les va- 
leurs propres sont réelles positives, le rayon spectral de A est égal à lime o(trAUE, comme 
on le voit aisément en mettant À sous forme triangulaire, de sorte que (4.6) est l’analogue du 
rayon spectral. En explicitant tr A, on peut écrire 


p(A) = en & (Ai te AE 


= in. 
En définissant le poids wA(c) du circuit (£1,... ,i) de longueur [(c) = k par 
Wal) = A AeAG S (4.7) 
on a 
p(A)=@ w(c)S . (4.8) 


p(A) est donc égal à la moyenne géométrique” maximale des poids des circuits du graphe. 
Si un circuit c = cjc2 est le produit de deux circuits élémentaires, on constate que le poids 
moyen de c est dominé par le plus grand des poids moyens de c, et de c2. On pourra donc 
se limiter aux circuits élémentaires pour calculer le max (4.6). En outre, comme les circuits 
de longueur supérieure à n sont composés, on pourra écrire, quitte à rajouter des termes 
redondants à (4.6) : 


OO 


p(A) = tu AN. 


k=1 


EXEMPLE 4.1.3.4. Pour la matrice À de la Figure 4.2, le graphe associé n’admet que les 
trois circuits élémentaires (1, 1), (3, 3) et (1,2, 3, 1), d’où 


A3 + An + Ai 
p(A) = max(A:11, 433, DE PRE. =3., 


EXERCICE 4.1.3.5. Soient P une matrice de permutation, D une matrice diagonale à coef- 
ficients diagonaux non nuls et «& un scalaire. Montrer les propriétés suivantes : 
1. o(œA) — ap(À). 
p(PAPT') = p(A) 
p(DAD”") = p(4) 
p(A) = ®; À; 
si À symétrique, alors p(A) = @;; À;; 
VkeN, p(A*) = (p(4)) 


AN OR D 


Nous attaquons maintenant la preuve de 4.1.3.3. D'abord, une remarque structurelle 
élémentaire. 


LEMME 4.1.3.6. Si À est irréductible, un vecteur propre u a toutes ses composantes 
différentes de &. 


2 


Démonstration. C’est un cas spécial du Lemme “de solidarité” 4.1.1.2, car Ra, est un se- 
mianneau positif. [] 


7 géométrique dans le dioïde, i.e. arithmétique dans l’algèbre usuelle 
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On rappelle que l’étoile de la matrice À est définie (dans un dioïde complet) par 


Dans un dioïde non complet, il faut que les matrices A admettent une borne sup pour que 
cette écriture ait un sens. 


LEMME 4.1.3.7. Dans le dioïde R!*", la matrice A* est bien définie Ÿ si et seulement si 
p(A) <e. 
Ce n’est que la version duale de la Proposition 2.2.5.2 du Chapitre 2. 


On découpe maintenant le théorème en deux parties. 
LEMME 4.1.3.8 (“SUR-VALEUR PROPRE”). Soit À € (Rnax )*”" irréductible, À € Rnax ef 
D(ÀA) défini comme en (4.6). Les deux assertions suivantes sont équivalentes : 
1. Ilexiste u € KR, \ {e} tel que Au =< }u, 
2. p(A) < À. 
Le vecteur u associé sera qualifié de “sous-vecteur propre” et À de sur-valeur propre. Le 


lemme affirme que lorsque À est irréductible, la moyenne géométrique maximale p (A) est la 
plus petite sur-valeur propre. 


LEMME 4.1.3.9 (“SOUS-VALEUR PROPRE”). Soit À € (Rmax )'*”", À € Rmax ef DP(A) défini 
comme en (4.6). Les deux assertions suivantes sont équivalentes : 


1. Iexiste u ER, \ {e} fel que Au Y Au, 
2. p(A) x À. 


Le vecteur u sera qualifié de sur-vecteur propre et À de sous-valeur propre. Le Lemme 
affirme que p(A) est sous-valeur propre maximale (même si À n’est pas irréductible). 


Avant de montrer ces lemmes, donnons un procédé simple pour construire des sous-vecteurs 
propres. 


NOTATION 4.1.3.10 (MATRICES À;, À). Soit À un scalaire non nul. On définit À, — 
À l'A. Lorsque p (A) n’est pas nul, on pose À — Ap(A) = (p(A))" A. 


LEMME 4.1.3.11. Si À > o(A) £ 6, alors toute colonne non nulle de la matrice (A;)* est 
sous-vecteur propre pour À. En particulier, toute colonne non nulle de A* est sous-vecteur 


propre pour p(A). 
Démonstration. On a par homogénéité p(A;) = À! p(A). La convergence de la matrice 
(A;)* résulte du Lemme 4.1.3.7. L’assertion découle de l’observation suivante” : 

A GAY = (4) = GA)" < (A). (4.9) 


n>1 
Ü 
Preuve du Lemme 4.1.3.8. (1) (2). On peut écrire : 
Vk, Alu < Mu, 
d’où 


8de manière équivalente, en se plongeant dans Rinax = Rimax U{+ 00}, tous les coefficients de A* sont < +00. 
%On rappelle que la matrice AT est définie par AT = A @ A7@... 
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Vi,k Afu; < Ab 


et comme tous les u; sont non nuls (via 4.1.3.6), on a en simplifiant 


Vi,k (AË)E < à, (4.10) 
soit en sommant : 
PA) = irA9f = BA < 2 ID) 
k=1 = 


(2)=(1) Réciproquement, le Lemme 4.1.3.11 garantit l’existence d’un sous vecteur propre 
lorsque À £ €. Si À = €, on a p(A) — €, la matrice À est sans circuits, donc nilpotente. Soit 
p le plus grand p tel que A? Æ &. On aura Au < eu en prenant u égal à une colonne non 
nulle de (4°). Cela montre le Lemme 4.1.3.8 (| 


Preuve du Lemme 4.1.3.9. (1) = (2). Choisissant u; £ 0: 
Au; < (Au); — Ai pli (4.12) 


pour un certain i réalisant le max dans (Au);. En appliquant le même argument à ÿ>, et en 
continuant de même, on obtient une suite {i,},-1 à valeurs dans {1,... ,n}, et donc on a des 
indices g et k < n tels que i,,4 = i,. On peut supposer q = i. On a alors : 


Au; < (Au); — Aiili,: Au, < À;,. nl, ae ui, < Ai ui, (4.13) 

ce qui, après des substitutions répétées, donne : 
RG < Ai Ai i gas Ai, jui A (4.14) 
Comme u; Æ 0, la moyenne géométrique du circuit (é, i2,... ,ix) est minorée par À, d’où 


À < p(A). 
(2)=(1). Supposons tout d’abord p(A) À &, et introduisons la notation suivante : 


NOTATION 4.1.3.12 (GRAPHE CRITIQUE). Les circuits qui réalisent le max dans (4.8) sont 
qualifiés de critiques. On appelle graphe critique, noté GC(A), le sous graphe de G(A) 
formé des sommets et arêtes appartenant à un circuit critique. On notera C° = (C1,...,C;) 
l’ensemble des g composantes fortement connexes du graphe critique. 


EXEMPLE 4.1.3.13. On a représenté sur la Figure 4.3 une matrice et son graphe. À est de 
rayon spectral 2. Le graphe critique est en traits gras. On a les deux composantes fortement 
connexes C} = {1} et C5 = {2,3}. 


KA ii 


FIG. 4.3: Graphe critique (en gras) 
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Soit À un indice appartenant à un circuit critique de À. Ce circuit a un poids e pour la 
matrice A. On a alors (A) — e, et donc les colonnes 


Ai = BA: À = BA: 


n>0 n>1 


qui ne diffèrent éventuellement que sur la diagonale en raison du premier terme de la somme 
(A9 = Id) sont égales. Il résulte de 


AG). = (À): = (À*).; 


n 


que le vecteur 4 — (À*).à est vecteur propre de À. Le cas p(A) = € se traite de manière 
analogue à 4.1.3.8,(2)—(1). Cela achève la preuve du Lemme 4.1.3.9. Le Théorème 4.1.3.3 
résulte immédiatement des Lemmes 4.1.3.8 et 4.1.3.9. (I 


Nous caractérisons maintenant le moduloïde!° des vecteurs propres de À. 
THÉORÈME 4.1.3.14. Soit À € (Rnax)”*" irréductible et C,..., C, les composantes for- 
tement connexes du graphe critique. Choisissons i1 € CT,...,i, € C,. La famille formée 
des colonnes suivantes de A* : 


(A sde 


He "lg 


} (4.15) 


est une famille génératrice minimale du moduloïde propre de A pour la valeur propre p(A). 
Le lemme suivant montre que la famille (4.15) ne dépend pas!! du choix des indices 

liens ours 

LEMME 4.1.3.15. Si i, et i, appartiennent à la même composante fortement connexe du 


De def + def + 4 
graphe critique, alors les vecteurs v;, = A? et vi — A*. sont proportionnels. 


Démonstration. Montrons que les deux colonnes A}. et A}, sont proportionnelles. Comme 


i, et i, appartiennent à la même composante fortement connexe du graphe critique, on a 
A*. Ai; = e, et donc, pour tout r : 


itis 


À}; = ÀŸ, À}, 
= AŸ: d'A 
= À}, 
d’où l’on conclut que 
Anse AT AS 


[] 


Il est donc sensé de choisir arbitrairement un indice par composante fortement connexe, 
comme on l’a fait. Nous attaquons maintenant la preuve du théorème. 


1/ La famille (4.15) engendre le moduloïde propre Posons v; — A: On prétend qu’un 
vecteur propre 4 pour la valeur propre p(A) vérifie 


10 
11 


1.e. module sur un dioïde, défini de manière analogue aux modules sur un anneau 
à une permutation et à la multiplication par des scalaires inversibles près 
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u = Àtu= Œoiu, (4.16) 
S 


ce qui montrera que la famille est génératrice. On a, comme les v;, correspondent à des 
colonnes de A* : 


Au = QE Au x Qu, | (4.17) 
k S 


On a d’autre part en choisissant des indices réalisant le max dans u — Au : 


Ui = Ai li 
Ui, = Àj,ili, 

En poursuivant de la sorte, on finit par obtenir un indice i, appartenant à un circuit c tel que 
wz(c) = e,1.e. un circuit critique. On a alors : 


A+ À .À À + 
A; ;ui = Ai Ai, ... À Uiy © A 


rik Ui, 


ik ik 
ce qui montre l’inégalité opposée à (4.17). 


2/ La famille (4.15) est génératrice minimale Ce résultat est un peu plus délicat et peut 
être sauté en première lecture. Pour montrer ce résultat de type indépendance linéaire, on 
fait appel à une structure algébrique symétrisée. On introduit le dioïde des couples, soit R2,, 
muni des lois 


(a,a')@(b,b)={(aga,b@æb'), (a,a')@(b,b') = (abæa’b',ab'@ ab). 


Sur R?.. est défini un signe moins : @(a, a’) = (a', a). On a @(x @ y) = (@x) & (6y), 
O(xy) = (@x)y = x(@y) et 6 @ x = x comme le lecteur le vérifiera aisément. On peut 
identifier l’élément a € Rx au couple (a,£s) € R2.- On plonge de la sorte Rx dans 
Rx: Moyennant cette identification, on peut écrire (a, a’) = a @ a’. Autrement dit, R?,, 
est obtenu à partir de Rx par doublement de structure, ce qui revient à manipuler un signe 
moins formel. On peut alors définir le déterminant d’une matrice B e RX* (à valeur dans 


R2,,) par la formule usuelle : 


det B = Æ sen(o) Q Bio) 
o i=1 


où sgn(o) vaut e si o est paire et @e sinon. D’ordinaire, si les vecteurs colonnes d’une 
matrice sont linéairement dépendants, alors le déterminant de la matrice est nul. Ici, on a de 
manière analogue : 


LEMME 4.1.3.16. Soit B une matrice carrée telle que B.x — Prix AB.1 Alors, le 
déterminant de B est de la forme a & a. 


Démonstration. Le déterminant se développant par multilinéarité comme d’habitude, il suffit 
de voir que le déterminant d’une matrice ayant deux colonnes identiques est de la forme a6a. 
Or un couple x est de la forme a 6 a ssi x — @x (car alors x = x @ x), donc il suffit de 
montrer que det B = € det B. Supposons B.; = B. ;, et soit 7 la transposition d’indices à, j. 
On a 
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detB = Æsen(o) GB 
— ee sgn(o ot) © Bicor(i) 
— ee sgn(o ot) © Bic (i) 


— @detB 
d’où le Lemme 4.1.3.16. (E] 


On peut maintenant conclure la preuve du point 2/. Si la famille n’est pas génératrice mini- 
male, on a une combinaison linéaire de la forme 


A = PAS, 


1£k 


Soit À la matrice extraite de À* formée des lignes et colonnes d’indices i1,... ,1,. Il suffit 
de voir que det À n’est pas de la forme a @ a. Tout les circuits de la matrice À étant de poids 
inférieur ou égal à e, il s’ensuit que pour toute permutation ©, @, Ac(s) < e (décomposer 
o en produit de cycles). En outre, l’égalité est atteinte seulement pour la permutation iden- 
tité (sinon, le poids d’un cycle non trivial serait égal à e, contredisant le fait que à, et à, 
appartiennent à des composantes fortement connexes différentes du graphe critique). Aïnsi, 
det À — e — (e, b) avec b < e et n’est donc pas de la forme a € a. [] 


EXERCICE 4.1.3.17. Valeurs propres et vecteurs propres de la matrice suivante 


£ —2 —3 
À = 1 £ 2 
—5 —4 2 


EXERCICE 4.1.3.18 (4). Montrer que 
p(A) — min{||DAD- || | D diagonale inversible} , 


où la “norme” d’une matrice M est définie par |M|| = &@, ; Mij. Indication : on pourra faire 
usage de certaines des inégalités de l’exercice 4.1.3.5. 


4.1.4 Extension aux matrices réductibles 


Cette section peut être sautée en première lecture. Nous traitons seulement le cas d’une matrice 
ayant deux composantes connexes. Quitte à remplacer À par P!A P (P matrice de permutation), on 
peut supposer que 


B C 
I (4.18) 


avec C À & (sinon, la matrice est bloc diagonale et on est ramené à la théorie précédente) et en outre 
que les matrices B et D sont irréductibles!?, soit le schéma ci dessous : 

D 

NT 
Op 
pe \ ) 


l2Le cas à plus de deux composantes connexes n’est pas essentiellement plus difficile 
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Au vu du dessin, le bloc D peut être qualifié de bloc amont (les arcs de C vont de D vers B). 


THÉORÈME 4.1.4.1. La matrice (4.18) admet les deux valeurs propres distinctes p(D) et p(B) ssi 
p(D) > p(B). Sinon, elle admet l'unique valeur propre p(B). 


Dans tous les cas, on retiendra que o(A) — max(p(B), p(D)) est valeur propre (le rayon spectral 
(max, +) est toujours valeur propre, même dans le cas réductible). 


Démonstration. En partionnant un vecteur propre u sous la forme [v, w], il vient 


Av = Bv@Cw 
Àw = Dw . 
1/ Cas w = &. On se ramène au problème spectral pour B. p(B) est donc toujours valeur propre de 
À. 2/ Cas w £ &. La seconde équation entraîne b(D) = X et la première entraîne d’après le Lemme 
relatif aux sur valeur propres que À > p(B). Il est donc nécessaire que p(B) < p(D). Alors, on 
introduit comme précédemment la matrice À = p(D)-! A, dont l'étoile converge. En raisonnant de 
même, on voit que toute colonne de À telle que À = — e est vecteur propre de À pour la valeur 
propre p (D). Le théorème en résulte. (n] 


EXEMPLE 4.1.4.2. On a représenté une matrice et sa transposée sur la Figure 4.4. On a pour À deux 
composantes fortement connexes C1 = {1}, C2 = {2}, C2 étant en aval de C1. D’après le théorème, 
on a donc l’unique valeur propre p (C2) = 3 associée à un vecteur propre de support {2}. Par contre, 
cette condition est vérifiée pour AT, et donc le spectre de AT est égal à {p(C1), p(C2)} = {1,3}. 
Contrairement au cas usuel, le spectre d’une matrice diffère du spectre de sa transposée. 


Ci :1l 


C2 aval, p(C2) 4 p(Ci) | C1 aval et p(C1) < p(C2) 
une seule valeur propre À = p(C2) = 3 deux valeurs propres p(C1) = 1, o(C2) = 3 


.. | € 
vecteur propre associé | | 
e e 


; e —3 
vecteurs propres respectifs | ë | | | 


FIG. 4.4: Condition de placement des valeurs propres 


4.2 Cyclicité 


4.2.1 Analogie avec l’algèbre usuelle 


On étudie maintenant le comportement asymptotique de la suite 4”. Ici encore, nous pen- 
sons utile de rappeler certains théorèmes asymptotiques relatifs aux matrices à coefficients 
positifs ou nuls dans l’algèbre usuelle. C’est donc encore un paragraphe culturel. Soit tout 
d’abord À € (R* \{0})"*". Sur le cas particulier de la Figure 4.1, on constate que l’image par 
Af du cône positif tend vers une droite vectorielle. En normalisant, on pose B = p(A) 'A, 
on peut montrer que B! converge vers une matrice fixe de rang 1. Lorsque l’on considère des 
matrices dont certains coefficients peuvent être nuls, la situation est plus compliquée : par 
exemple, pour À matrice de permutation d’ordre p, la suite {A*}£en cycle sur p valeurs. La 
notion pertinente qui permet d’éliminer cette difficulté est la suivante : 
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DÉFINITION 4.2.1.1. On appelle cyclicité de la matrice irréductible À, notée c(A), le pgcd 
des longueurs des circuits du graphe associé à A. Lorsque c(A) = 1, la matrice À est dite 
primitive. 

EXEMPLE 4.2.1.2. La matrice À de la Figure 4.5 à pour cyclicité 2. La matrice B est pri- 
mitive. La matrice de permutation C est cyclique d’ordre 3 (la cyclicité est égale à l’ordre de 
la permutation). Le lecteur notera enfin qu’une matrice irréductible dont l’un des coefficients 
diagonaux est non nul est primitive. 
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FIG. 4.5: Matrices et graphes associées 


La suite des puissances d’une matrice primitive est ultimement positive : 
LEMME 4.2.1.3. Soit À € (R*)"*". Les propriétés suivantes sont équivalentes : 


1. À est primitive 
2. existe NENKZ=N=Vi,jef(l,...,n}, A > 0. 


Démonstration. (2)=(1) On a en effet AY > 0, AN #1 > 0,..., soient des circuits de lon- 
gueur N,N +1,N +2,...,donc c(A) = pgcd(N,N +I,N +2,...)=1. 

(1)=(2) Soient c1,...,c, des circuits!? tels que pgcd(|c1l, ..., lc) = 1. Choisissons des 
points i1,...,i,, respectivement dans c1,...,c,. Le graphe associé à A étant fortement 
connexe, 1l existe un chemin p, de j à ÿ,, un chemin p; de ij à à, ... , un chemin p,y1 
de i, à i. Pour tout « € N?, on a donc le chemin de j à i : 


Pa = PiCT Pac Pa 
de longueur 
Pal = Ipil + lat +... +alc| (4.19) 


Comme les |c;| sont premiers entre eux, il résulte facilement du théorème de Bezout que tout 
k assez grand est égal à |P,| pour un certain qg-uplet d’entiers naturels &, et donc que A; > 0 
pour k assez grand. Ü 


Nous donnons pour mémoire le résultat suivant [35, 64], que nous ne prouvons pas : 


THÉORÈME 4.2.1.4. Soit À une matrice primitive et P le projecteur spectral\* associé à la 
valeur propre p(A). Alors 


k k k 
A° = (p(A)) P +o((p(A))) . (4.20) 
BOn note [c;| la longueur d’un chemin ou circuit c; 
en général, si Àj,...,AQ sont les valeurs propres de À de multiplicités respectives n1,...,nQ, on a la 


décomposition en somme directe des espaces caractéristiques : R° = ker(ài — A)"! +...+ker(ky — A)"4. Le 
projecteur x; d’image ker(1; — A)"; et de noyau la somme des autres espaces caractéristiques est le projecteur 
spectral associé à À;. Ici, la situation est très simple, puisque le rayon spectral d’une matrice primitive est racine 
simple du polynôme caractéristique [35, 64]. Soient w et v des vecteurs propres de À respectivement à gauche et 
à droite (wA = p(A)w et Av = p(A)v), alors on montre que P est donné par P = (wv)l Uw. 
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En mettant À sous forme de Jordan, ce théorème revient à montrer que pb (A) est racine 
simple du polynôme caractéristique et que toutes les autres racines sont de module stricte- 
ment inférieur. 

Le proposition suivante ramène l’étude asymptotique des matrices irréductibles non pri- 
mitives à celle des matrices primitives. 


PROPOSITION 4.2.1.5. Soit À irréductible. La matrice A est bloc diagonale, chaque 
bloc étant primitif. 

LEMME 4.2.1.6. Les longueurs de deux chemins de mêmes extrémités sont congrues modulo 
c(À). 

Démonstration. Soient p, et p2 deux chemins de j à i. À étant irréductible, on a un chemin 
r de i à j, donc 


rpil=lrl+ipil=0 (mod c(A)) 
Irp2l = rl +1p2l = 0 (mod c(A)) 


et en faisant la différence |p1| = |p2| (mod c(A)). (| 


Introduisons la relation d’équivalence suivante sur {1,... ,n}:“i R j ssi la longueur 
de tout chemin de j à i est multiple de c(A)”. 

On conclut maintenant la preuve de la proposition. Singularisons un point i. On a c(A) 
classes d’équivalences pour la relation R, Co, ... , Ceça)_1, caractérisées comme suit : 


j € C & la longueur de tout chemin de j à à est congrue à / modulo c(A) 


On constate en outre qu’un circuit parcoure nécessairement toutes les classes successivement 
(quand on rajoute un arc, on passe dans la classe suivante). Ainsi, d’une part, tout chemin de 
longueur multiple de c(A) aboutit nécessairement dans la classe d’où il part, ce qui montre 
que A‘ est une matrice bloc diagonale, et d’autre part, l’ensemble des circuits passant 
par une classe donnée est égal à l’ensemble des circuits de À, ce qui montre que le pgcd 
des longueurs de ces circuits est égal à c(A), donc que la cyclicité d’un bloc irréductible 
quelconque de A4) est égale à 1, ce qui montre la proposition. (m 


4.2.2 Cryclicité des matrices (max,+) 


4.2.2.1 (Cyclicité dans le cas (max,+)). On définit comme suit la cyclicité d’une matrice 
A = on xXn 
ax * 
1. La cyclicité d’une composante connexe du graphe critique est égale au pgcd des lon- 
gueurs des circuits de cette composante. 


2. La cyclicité de la matrice À est égale au ppem des cyclicités des composantes connexes 
du graphe critique. On la notera c(A). 


Le théorème suivant étend la formule asymptotique 4.2.1.4 au cas (max, +). La situation 
(max, +) est plus précise et plus simple, puisque le régime asymptotique est atteint en un 
temps fini au lieu d’être atteint asymptotiquement. 


THÉORÈME 4.2.2.2 (FONDAMENTAL). Soit À € R'*" irréductible. Il existe un entier N tel 
que 


k> N = AË+@ 2 (p(A)) 4) A . (4.21) 


LEMME 4.2.2.3. Tout circuit du graphe critique est critique. 
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Preuve du Lemme. Soit © = (ii,... ,ix) circuit du graphe critique. On a donc un chemin 
P1 tel que p1 U (iii) soit un circuit critique, p> tel que p2 U (i2i3) soit un circuit critique, 
et ainsi de suite. Le circuit composé p1 ...p,c est donc critique. Si le poids moyen de c 
était inférieur strictement à p(A), alors le poids moyen du circuit p1 ... px serait supérieur 
strictement à p (A) : absurde. [] 


On traite d’abord le cas de matrices de cyclicité 1. Quitte à normaliser À, on peut sup- 
poser p(A) = e. 
LEMME 4.2.2.4. Soit A telle que p(A) = e et c(A) = 1. On a At! = AË pour k assez 
grand. 


Démonstration. 1/ Pour i appartenant au graphe critique, AË = A — e à partir d’un 
certain rang. Soient en effet c1,...,c4 des circuits de la composante connexe du graphe 
critique contenant à tels que 


pecd(lcil, ... , cel) = 1 . 
Il existe un circuit & passant par i et rencontrant chacun de c. On en déduit pour tout 
(n1,... ny) e N° l’existence d’un circuit passant par i et d’image commutative 
Benne à 


Via le Théorème de Bezout, pour tout n assez grand, il existe un k-uplet d’entiers 
(1,...,nx) tel que n = [al + nilcil +... + nylcxl. Ainsi, pour n grand, on à un cir- 
cuit de longueur n passant par i. D’après le Lemme 4.2.2.3, un tel circuit est critique, d’où 
Are: 

2/ Pour j appartenant au graphe critique et i quelconque, À, — A; à partir d’un certain 


rang. Soit en effet ko tel que Aï = A. Ona 


pour k assez grand d’après le point 1/. Comme l’autre inégalité est triviale, le point 2/ est 
acquis. 


3/ Pour j n'appartenant pas au graphe critique, on va montrer qu’à partir d’un certain rang : 
k k A+ 
af. = ee AA (4.22) 
i critique 


Soit g tel que AŸ, — Aï. On a pour k assez grand 


tJ 


i critique 


et comme d’après 2/ pour k assez grand, A2 — A, on obtient pour & grand : 


i critique 


On montre maintenant que pour & grand, les chemins réalisant le max dans AF. passent né- 
cessairement par le graphe critique. On note n le poids moyen maximum des circuits non cri- 
tiques (le ‘deuxième poids moyen”). Un tel chemin p se factorise en un chemin élémentaire 
(de longueur au plus n — 1) et des circuits de poids au plus 7, d’où la majoration 
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w(p) < "TA | 


On en déduit que pour k grand (ne dépendant pas de p), w(p) est dominé par 


D As; 


lcritique 
ce qui montre l’égalité (4.22). Ü 


Cas d’une matrice de cyclicité quelconque On pose B — A‘). D’après la Proposition 
4.2.1.5, B est bloc diagonale, les blocs étant irréductibles de cyclicité 1. On peut donc écrire 


k > N = pri — p(B)B* — Afc(A)+c (A) _ p (A) AË‘(A) 
d’où 
k > c(A)N = A#+c(A) _ p(A)"® AË 


ce qui achève la preuve du théorème. Ü 
EXEMPLE 4.2.2.5. Soit 


On a 
-n 9 
: sin" < 3 
; —1 e 
de 2 
Gas sin" > 3 
—1 _e 


Ainsi, le régime périodique peut être atteint au bout d’un temps arbitrairement grand. 


4.2.3 Application : régime périodique et taux de production 


Soit un graphe d'événements temporisé fonctionnant en régime autonome. Un tel graphe 
se représente par l’équation récurrente x(k) — Ax(k — 1). Lorsque A est irréductible, le 
théorème de cyclicité entraîne que le comportement est périodique, 1.e. 


x (k+c) = p(A)'x(k) 
pour k assez grand. En revenant à l’algèbre usuelle, cela s’écrit 
x(k+c)=cx p(A)+x(k) . 


Autrement dit, c événements se produisent toutes les co(A) unités de temps, et p(A) est le 
temps moyen de production d’une pièce (l’inverse du taux de production). Pour conclure, on 
peut donc énoncer le résultat suivant, fondemental pour la compréhension des systèmes de 
production : 
un graphe d'événements temporisé fortement connexe fonctionnant en régime autonome 
atteint après un régime transitoire fini un régime périodique, le taux de production étant 
l'inverse de la valeur propre de la matrice de la dynamique (max, +) du système. 
ISL’extension au cas réductible ne présente pas une grande difficulté mais est fastidieuse (on décompose À en 


composante irréductibles, à chaque composante est associée un régime périodique d’après ce qui précéde, etc.). 
Dans la section suivante, nous établissons un résultat général plus faible 
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4.2.4 Temps de Cycle 


On vient de voir que le théorème de cyclicité répond complètement au problème de 
déterminer le comportement asymptotique de la suite x(k) — A*x(0), quand la matrice À est 
irréductible. Une conséquence bien remarquable de l’hypothèse d’irréductibilité est que si 
x(0) a toutes ses composantes finies, limz x(k);/k est égal à la valeur propre de À, pour tout 
i. Autrement dit, toutes les composantes du système fonctionnent au même rythme. Une telle 
condition n’est pas toujours réalisée en pratique. Par exemple, dans une ligne de production 
gérée naïvement, il se peut très bien que les cellules en amont produisent plus rapidement que 
les cellules en aval, ce qui fait d’ailleurs exploser le stock intermédiaire (d’ou le qualificatif 
de naïf pour la gestion). On voudrait bien sûr détecter ce genre de phénomènes, quand ils ont 
lieu. 

Il est donc naturel de calculer l’asymptotique de Afx(0) lorsque À est réductible. Dans 
le semi-anneau Ra (mais aussi, mutatis mutandis, dans le semi-anneau (R*,+, x))ona 
la pathologie suivante. Considérons la matrice À de la Figure 4.6. . On trouve les valeurs 


“2 | £ € € € € 
Gr ET . £ & | € € 

7 
TS ES bons 
us eme 0e £ 8 e 2e 
1 | É € EE E € 


FIG. 4.6: Le comportement de A dépend de la parité de k 


suivantes de la suite {A }k>0 : 8,8,0,1,—4,3,—8,5,—12,7,—16, ... On peut le voir graphi- 
quement : Af. est égal au chemin de poids maximal de 5 à 1 dans le graphe associé à À, 
représenté à gauche de la figure. Il est clair sur le dessin que si k est pair, un tel chemin passe 
nécessairement par la composante connexe {4} de rayon spectral égal à —2 et pas par la com- 
posante connexe {2, 3} de rayons spectral égal à 1. Au contraire, les chemins de longueur 
impaire passent par cette composante connexe de rayon spectral 1, ce qui explique les taux 
de croissances différents de Af, suivant la parité de k. 

Une manière d’exclure cette pathologie est de recourir à une opération de moyennisa- 
tion, bien naturelle dans les applications aux systèmes à événements discrets. Le vecteur 
initial x(0) n’a pas en général des coordonnées arbitraires dans R;X : on peut la plupart du 
temps supposer que les coordonnées de x(0) sont finies (des coordonnées infinies ne sont pas 
dénuées de sens physique et d’intérêt, prendre x(0); — —co, c’est faire démarrer le nœud 
i avant tout le monde, i.e. lui donner infiniment d’avance, mais il s’agit d’une autre histoire 
que celle qui nous occupe ici). 

Sur l’exemple de la matrice À ci-dessus, le lecteur vérifiera que lim,(A*x(0));/k = 1, 
pour tout À, dès que tous les coefficiens de x(0) sont finis : on a ainsi moyenné les compor- 
tements cycliques des différents coefficients de la i-ième ligne de A* en multipliant à droite 
par x(0). 

Cette observation nous conduit à définir la quantité suivante, baptisée temps de cycle par 
Gunawardena [47]. 


DÉFINITION 4.2.4.1 (TEMPS DE CYCLE). Supposons que la matrice À € (Ryax )"*" a au 
moins un coefficient fini par ligne. On appelle temps de cycle de À, la limite 


1 
x(A) = lim Tr * (AËx) , (4.23) 
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où x € R" est un vecteur quelconque. 


La notation è x (A x) est une hybridation un peu monstrueuse des notations classiques et 


des notations max-plus (la division est dans l’algèbre usuelle, la puissance de matrice en max- 
plus). Les puristes ou les esthètes préférons écrire, en max-plus : y; (A) — limy_oo(AËx)}/ d 


A 
Comme le remarque Gunawardena, la notion de temps de cycle est intéressante plus 

généralement pour toute application f : IR° — R" qui vérifie les propriétés d’homogénéité 
et de monotonie déjà évoquées au Chapitre 1 : 

(H)  fA+x)=1+/f0), 

(M) x<y= fa) <fO), 
pour tous À € R, x, y e R”, avec, par définition, À + x = (À +x1,... ,1+Xx,),etx < ysi 
x; < y, pour i = 1,...,n. Cela résulte en effet de l’observation élémentaire suivante due à 
Crandall et Tartar. 


PROPOSITION 4.2.4.2. Si f : R' — KR’ vérifie (H) et (M), alors elle vérifie : 
ON) IfG)—fOI<Ix- 71, 


où || : || désigne la norme sup. Réciproquement, si f vérifie (H) et (N), alors elle vérifie (M). 
Autrement dit, sous l’hypothèse d’homogénéité, les propriétés de monotonie et de non- 
dilatation pour la norme sup sont équivalentes. Avant de montrer cette proposition, nous en 


déduisons sinon l’existence du temps de cycle y, du moins le fait que x ne dépend pas de la 
condition initiale x € R". 


COROLLAIRE 4.2.4.3 (UNIVERSALITÉ DU TEMPS DE CYCLE). Si f : IR' — KR" vérifie 
(H) et (M), la limite x(f) = lime KT! FA (x), si elle existe pour un x € KR' particulier, 
existe pour tout y € R”, et a la même valeur. 

En effet, d’après la proposition qui précède, f vérifie (N), et donc || f(x) — f(>)|| < 
x — y]. Le corollaire en résulte aussitôt. Passons maintenant à la preuve de cette proposition. 


Démonstration. On peut écrire : 
—Ix—yl+7<x<y+lx- y. 
Sous les hypothèses (H) et (M), on en tire : 
—|lx — ll + FO) < fG) < FO) +]Ix —yl , 


donc [| f(x) — f(y)Il < [x — y]. On a montré la première assertion de la proposition : 
((H) et (M)) = (N). 

Passons à la seconde assertion. Si x < y,ona0 < y—x < ||y—x||, et donc —|y —x| < 
y—x—{||y —x| < 0, ce que nous réecrivons : 


x£<y=|y-@+1y—xlh < Ir xl . (4.24) 
En supposant en outre que l’hypothèse (N) est vérifiée, on tire de (4.24) : 
x £y > |fO)—- f@+Iy—xIDI <Iy — xl : 
En utilisant (H), il vient 
x £y > 1|fO) — (Cf) +17 —xII < [y —xl . 


On peut réécrire la dernière igégalité sous la forme [[U — W|] < W, avec W = ||y — x] et 
U = f(y) — f(x). On a alors : U > W — [U — W| > 0, ie. f(y) > f(x). On a montré 
((H) et (N)) = (M). Ü 
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En général, les hypothèses d’homogénéité de de monotonie (H) et (M) ne suffisent pas à 
garantir l’existence du temps de cycle x(f). Dans le cas des applications max-plus linéaires, 
le temps de cycle existe bien et se calcule comme suit. On note (A); le max des poids 
moyens des circuits qui sont en amont de i dans le graphe associé à A (rappelons que le 
poids moyen du circuit (é1,... ,ig,t1) Vaut k (A, ++ AÀ;;)). 

PROPOSITION 4.2.4.4. Si À € (Rnx)'*”" a au moins un coefficient fini par ligne, on a 
x (A) = (A). 

En attendant la version suivante du polycopié (c’est un trou), la preuve —-pas très longue— 

est laissée en exercice au lecteur (voir [47]). 


4.2.5 Algorithme de Karp pour le calcul du temps de cycle et du rayon spectral 


Le rayon spectral de la matrice À peut se calculer simplement. Un algorithme naïf 
consiste à évaluer!6 la formule (4.6) ce qui suppose la tabulation des puissances succes- 
sives de la matrice À jusqu’à l’ordre n (= taille de À). Pour une matrice pleine, on a ainsi un 
algorithme de temps d’exécution O(n*). 

Mais on peut faire beaucoup mieux. Une méthode plus efficace (en temps O(n*)) per- 
mettant de calculer la valeur propre d’une matrice consiste à appliquer une formule due 
à Karp [49], que nous présentons ci-après. En outre, une autre classe de méthodes, dites 
d’itérations sur les politiques, battent l’algorithme de Karp d’un ordre de grandeur avec un 
temps moyen d’exécution expérimentalement pratiquement linéaire en le nombre de coef- 
ficients de la matrice (=O(n?)), mais elles sont plus longues à décrire, et on ne sait pas 
prouver, pour l’instant, que le temps d’exécution dans le pire des cas est polynomial, voire 
O(n*), comme les expériences le laissent conjecturer. Nous renvoyons provisoirement le lec- 
teur intéressé à la section relative à l’algorithme de Howard sur la page web de l’auteur 
(http ://amadeus.inria.fr/gaubert). 

Karp a énoncé son algorithme dans le cas des matrices irréductibles. Dans le cas des ma- 
trices réductibles, nous allons montrer que le même algorithme calcule en fait naturellement 
le temps de cycle. 

Le cœur algébrique de la formule de Karp est le résultat suivant “à la Wielandt”, analogue 
à la caractérisation du rayon specttal des matrices irréductibles comme sous-valeur-propre 
maximale ou sur-valeur-propre minimale donnée plus haut. 


THÉORÈME 4.2.5.1 (SUPER-CONVOLUTION). Soit y € (Rnax )'*! \ {€}, ef À € (Rax )*" 
une matrice ayant au moins un coefficient fini par ligne. Posons 


A={XER] 74" < @ xy4"#. 


1<k<n 
On a 
max Xi (A) = min À . (4.25) 
VIRE 
Démonstration. Notons x la quantité à gauche de (4.25). Pour tout j € {1,... ,n}, on peut 
choisir i € {1,... , n} tel que (yA”); = y;4;;, et un chemin p de longueur n de j à i tel que 


A; = w(p). Un tel chemin passe au moins deux fois par un même sommet, que nous notons 
l. Aïnsi, on peut écrire p comme la concatenation des chemins p2cp;1, où c est un circuit 


1 
160n se souviendra que ak en (max,+) vaut 1/k x a avec les notations usuelles, de sorte que (4.6) n’est autre 
que la moyenne arithmétique maximale des circuits de la matrice À. 
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de longueur & > 1 passant par /, p, un chemin de j à / (peut être vide), et p, un chemin 
de / à i (peut être vide). Puisque À est en amont de /, on a, en utilisant la caractérisation 
du temps de cycle donnée dans la Prop. 4.2.4.4, w(c) < n°. Donc, w(p) < HA", et 
(YA); < A (y A") ;. Une telle inégalité ayant lieu pour tout j (le cas échéant avec un k 
différent), on a montré que x € A. 

Il reste à voir que x est le plus petit élément de A. Prenons un quelconque À € A. On a 


yA" < & AA 
1<k<n 


Par une récurrence immédiate, on prouve : 


Vp > 0, yA"TP < ee) APP MATE , 


1<k<n 
En multipliant à droite par x € KR”, on obtient : 


yA"TPx < Œ& RÉ AT- E 


1<k<n 
D'où on tire facilement : 
dis 
max x:(A) —= max lim (4"*Px)"? 
1<i<n 1<i<n po 
NiRE Vie 
1 
n+p 
UE 
— lim | max (4"*?x); = lim (yA"*Px)mr <X. 
L 
po À'1<i<n P—00 
VIRE 
Ainsi, u < À, pour tout À € À, ce qui achève la preuve. O 
U = 


COROLLAIRE 4.2.5.2. Sous les hypothèses précédentes, on a 


OA"); — GATN); 


max x;(A) = max min (4.26) 
1<i<n 1<j<n 1<k<n k 
VAE OA") jÉE 


Avec les notations max-plus, cela s’écrit : 


(A) = 2) 
sa ou a A (or | (4.27) 


Vie OAÏ);£e 


Démonstration. Les assertions suivantes sont clairement équivalentes 


yA" < ee) AT 


1<k<n 
OA"); Ze kefl,...,n}, (A1); < ME (yA"E) 
: A"); \? 
GA"); Ze kEf{L,...,n}, (LL) 


1 

A); \E 
1<j<n 1<k<n CAE) 
GA) 


Ainsi, on voit que le Corollaire 4.2.5.2 n’est qu’une simple reformulation du 
Théorème 4.2.5.1. (] 
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Par spécialisation du Corollaire 4.2.5.2, nous optenons la formule donnée par Karp pour les 
matrices irréductibles. 


COROLLAIRE 4.2.5.3 (FORMULE DE KARP). Pour toute matrice À € (Ra) *" 
irréductible, et pour tout i € {1,... ,n},ona: 

A SAR 

PCA) = xi(A) = max min 


<j<n 1<k<n 
A’; #E 


Démonstration. Appliquer le Corollaire 4.2.5.2 au vecteur y tel que y; = e et y; = € pour 
j Ai. [] 
Cette formule permet aisément de calculer l’une quelconque des coordonnées du temps de 
cycle et a fortiori le rayon spectral o (A) : il suffit de prendre y comme ci-dessus, de mettre 
les vecteurs y A, k = 1,... , n dans un tableau, et d’appliquer la formule (4.26), ce qui exige 
un espace mémoire O (n°) et un temps O(n*). Dans le cas où la matrice est codée de manière 
creuse, en notant Æ le nombre d’arcs du graphe associé, il faut un temps O(nE). L'espace 
mémoire naïf O (n°?) est rédhibitoire pour bien des applications pratiques (pour fixer les idées, 
n = 10°, E = 3n). Nous laissons le lecteur trouver la petite astuce qui permet d’evaluer la 
formule de Karp avec seulement une place mémoire O(n) (et un temps d’exécution encore 
O(nE), mais cependant doublé par rapport à l’évaluation naïve, on n’a rien sans rien). 


4.2.6 Calcul direct du taux de production à partir du graphe d’événements 


Dans ce qui précède, on a associé à un graphe d’événements donné une matrice À € 

ms et montré que le taux de production était égal au rayon spectral de À. On peut aussi 
caractériser directement le taux de production à partir du graphe sans avoir à construire la 
matrice À, aussi le résultat suivant nous semble-t-il le plus important de ce chapitre en ce qui 


concerne les applications. 
THÉORÈME 4.2.6.1 (FONDAMENTAL). Le taux de production À d’un graphe d'événements 
temporisé fortement connexe est caractérisé par 

N(c) 


c circuit élémentaire T'(c) ’ 


À = 


(4.29) 
où N(c) dénote le nombre total de jetons du circuit et T (c) la somme des temporisations des 


places du circuit. 


Démonstration. On a vu au Chapitre 3,$3.2.5 qu’en l'absence d’entrées, un graphe 
d’événements admettait les équations suivantes dans Rx 


x(k) = ASAix(k — 1)... ASA,;x(k — a) , (4.30) 


où x désigne le vecteur des fonctions dateurs. Cela s’écrit de manière équivalente sous la 
forme d’état X(k) = AX(k — 1) avec 


| x(E) | AjAi AÿA42 .. AjA 
Id E me € 
X( = Az 


a+ | ea UE 
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Il importe ici de se souvenir comment les matrices Ao,..., A, sont obtenues à partir du 
graphe d'événements initial. Supposons par exemple une unique place P sans jetons avec un 
temps de séjour de f, unités de temps, située immédiatement entre les transitions j et i. On a 
l’inéquation x;(k) > x;(k) +f,, autrement dit, le coefficient #, associé à une place à 0 jetons 
apparaît dans la matrice À, : on a A’, = 1). De même, (A;);; correspond à la temporisation 
maximale des places situées immédiatement entre j et ; et contenant un unique jeton, A) 
correspond aux places contenant 2 jetons, et ainsi de suite. À chaque arc du graphe de À, 
nous associons un chemin du réseau de Petri de la manière suivante. D’une part, à un arc 
correspondant à un coefficient de la forme (A$A,);; de À, nous associons un chemin du 
réseau de Petri, partant de la transition j, empruntant une place avec r jetons, poursuivant 
par un nombre quelconque (fini, peut être nul) de places sans jetons, finissant en la transition 
i, et de poids total (A5A,);; (le poids est ici le temps total de séjour). L'existence d’un tel 
chemin résulte clairement de l’interprétation de AS A, en termes de graphe (ce chemin n’est 
en général pas unique). D'autre part, à un arc correspondant à un coefficient de la forme Id;; 
de À, on associe le chemin de longueur et de poids nul dans le réseau de Petri, allant de à à 5. 
Plus généralement, à un chemin p de À, on associera le chemin du réseau de Petri formé en 
mettant bout à bout les chemins ainsi associés aux différents arcs de p. Pour tout circuit c de 
longueur k de À, il est clair que le terme w(c)/|c| est égal à T(c’)/N(c’) pour le circuit c’ du 
réseau de Petri ainsi associé à c. Réciproquement, si c’ est un circuit vérifiant l’égalité dans 
(4.29), c’ est associé à un circuit c de À par la correspondance décrite. La formule en résulte, 
moyennant l’observation que la forte connexité du réseau de Petri entraî ne la forte connexité 
du graphe de À, laquelle est un sous-produit immédiat de cette correspondance. Ü 


EXERCICE 4.2.6.2. Quel est le taux de production du graphe d'événements de la Figure 
3.12 du Chapitre 3 ? 


4.2.7 En résumé 


Nous pensons utile de conclure cet exposé de la théorie spectrale (max, +) en offrant un 
tableau comparant les résultats avec ceux de la théorie de Perron-Frobenius usuelle. 


CT em 


irréductible = 
unique valeur propre 
graphe critique irréductible = 
une seule droite propre 


irréductible = 
unique valeur propre positive 
t.q. vecteur propre associé positif. 


Propriétés 
spectrales 


p(A) = max{r | Au € (RŸ)"\{0}, Au > tu} | p(A) = max{r | 3u € (Rmax)! \{e}, À Qu > tu} 


1/ si une seule composante f.connexe 
du graphe critique, 

c(A) = pgcd longueurs 

des circuits critiques 

2/ en général, 

c(A) = ppem des cyclicités 

des composantes f. connexes 

du graphe critique 


Propriétés cyclicité 
asymptotiques c(A) = pgcd des longueurs 
(cas irréductible) des circuits de A 


A1c(4) u p(A}"<(@) p n>N> An+c(4) ue p(A)"(4) Ar 
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4.2.8 Quelques mots sur l’extension au cas probabiliste : exposant de Lyapu- 
nov des produits de matrices aléatoires 


Il s’agit ici d'indiquer brièvement comment ces résultats peuvent se généraliser au cas 
des matrices aléatoires. Nous renvoyons le lecteur à [4] pour un traitement détaillé. 

L'étude du comportement asymptotique de A" était motivée par le fait que les co- 
efficients de À représentent essentiellement les durées (fixes) de certaines tâches. Plus 
généralement, on peut s’intéresser au cas où ces durées sont aléatoires. Soit x(n) le vecteur 
des dateurs du graphe. On aura pour le système autonome une représentation dans l’algèbre 
(max, +) du type 


x(n) = An — 1)... A()A(O)x(0) , 


où les A(k) sont des matrices aléatoires : À;; (k) représente le temps minimal (aléatoire) de 
transfert du jeton numéroté k entre les transitions j et i, et l’on s’interroge sur l’existence de 
la limite de l’espérance 


def. < 1 
L = lim E[- x x(n)] (4.31) 
n 


n— 00 


qui est l’inverse du taux de production moyen (1.e. Le représente le nombre moyen de tirs 
de la transition i par unité de temps). Dans le cas où lim, E[x(n + 1) — x(n)] existe, on 
constate immédiatement (par une moyenne de Cesaro) que L = lim, E[x(n + 1) — x(n)]: 
L est alors le le temps moyen entre deux événements, ce qui justifie le nom traditionnel de 
temps de cycle. Nous ferons pour simplifier les hypothèses suivantes : 


HYPOTHÈSES 4.2.8.1. 


1. Les matrices aléatoires A(0), A(2)... sont i.i.d. (c’est-à-dire indépendantes identi- 
quement distribuées). 

2. Les variables aléatoires A;; (C9 _ Ai; (k) @ e = max(A;; (k), 0) sont intégrables (i.e. 

ELA; (K)*] < oo). 

3. La condition initiale x(0) est une variable aléatoire telle que Elmin; x(0);] > — et 

E[max; x(0);] < +oo. 


Considérons dualement À;; (k) 7 _ min(A;; (k), 0). On a A;;(k) = Ai;(k)* + A;;(k)., 
donc l’hypothèse (2) nous permet de définir sans ambiguïté E[A;; (&k)] = E[A;; (R)*] + 
E[A;; (k)7], éventuellement égale à —c (l'hypothèse d’intégrabilité a éliminé le cas am- 
bigu +00 + —oo). Cette dissymétrie entre +oo et —oo = € vient du fait que & est une “bonne 
valeur” qui correspond à l’absence de certains arcs dans le graphe, et qu’il est essentiel de 
pouvoir traiter, alors que + correspond à des temps moyens d’attente infinis, qu’il est na- 
turel d’exclure dans une première étude. 


Considérons maintenant 


x C1 xt) 


(autrement dit, on ne regarde que le sup des coefficients de x(n), soit la date la plus tardive 
d’occurrence de l’événement nr à tous les nœuds du graphe). On se contentera du résultat 
élémentaire suivant qui montre l’existence d’un “taux de production” dans le cas stochas- 
tique. 
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THÉORÈME 4.2.8.2. Sous les hypothèses 4.2.8.1, on a l'existence de 


ET LSeUTR 


(4.32) 


Démonstration. Posons plus généralement étant donnée une matrice À non nécessairement 
g def 
carrée : || A|| = ; A 
LEMME 4.2.8.3 (SOUS-ADDITIVITÉ DE LA NORME). On a pour toutes matrices À, B de 
tailles compatibles : 
IA @ BI < IAÏ & BI = AN +181 
Démonstration. || AB|| = @;;(AB);; = x Aix Biÿ < Ou Al Bx; < IA BP By; = 
IAB. Û 
Introduisons la notation suivante pour le produit de matrices 


Zuo © A(n — 1)... A(DA(O) . 


Il vient à partir de 4.2.8.3 : 
IX GI = 1 Zaox (OO) < IZroll @ IIx(O)| , (4.33) 
et d’autre part 


IxGI Z 1Zroll & min x); = Zroll + minx(O) , (4.34) 


d’où moyennant l’hypothèse 4.2.8.1,(3) 


.…. Hx@)  .  [Zrll 
Lim ——— = [jim—— . 


n n n n 
Introduisons plus généralement pour p > q 


def — 
Zpg = A(D—1)...A(@) , 294 = ElZrqll : 


On a pour p > r > 0, 
Zp0 = ZprZro ; 
donc, par sous additivité de la norme 
Zpq € Zpr © Zr0 = Zpr + Zr0 : (4.35) 


Le point essentiel est que Z, , est de même loi que Z,_,, (cela résulte du caractère 1.1.d des 
A(k)), donc 


Zp,r — Lp-r0 : (4.36) 


Posons maintenant 


def s 
Un — no, de sorte que Ellx(n)|| = v, . 


Il résulte de (4.35) et de (4.36) que v, vérifie l’inégalité sous additive 
Vprq € Up + Ug : (4.37) 


De telles suites vérifient la propriété suivante, qui est bien classique. 
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LEMME 4.2.8.4. Soit v, € RU {—co} une suite sous-additive. Alors, lim, v,/n existe dans 
RU {—co}, et en outre 
. LU 
lim— = inf — . 
nn 
La preuve est laissée au lecteur en guise d’exercice (élémentaire). On déduit aussitôt de 
ce lemme l’existence dans R U {—co} de la limite suivante, dite exposant de Lyapunov 


1 
&= lim 2 = lim El A(p — 1)... A(0)|? = lim— x El A(p — D @...@ AC) . 
PP _p P P D 
(4.38) 


[] 


En fait, on peut montrer beaucoup mieux, à savoir que jus converge vers l’exposant de 


Lyapunov £ de manière presque sûre (via le théorème ergodique sous additif de Kingman), 
et même que lim, x(n);/n = £ p.s pour tout i moyennant des conditions d’irréductibilité. 
Les hypothèses de type i.1.d. peuvent aussi être naturellement relaxées en des hypothèses 
d’ergodicité et de stationnarité [4]. Cette théorie constitue une extension au cas max, + de 
la théorie des produits de matrices aléatoires maintenant classique [13, 34], à laquelle le 
terme “exposant de Lyapunov” a été emprunté. Signalons que le problème du calcul effectif 
de l’exposant de Lyapunov £ est une question ouverte. On sait seulement traiter des cas 
particuliers de manière exacte, et donner des bornes dans les cas généraux. 


EXEMPLE 4.2.8.5. Dans le cas où A(0), A(1) ... sont des scalaires, (4.38) est une version 
affaiblie de la loi des grands nombres qui affirme que 


AG D+...+40 
P 


E[A(O)] p.s. 


& = lim(A(p — 1). A(0))? = lim 
P P 


EXEMPLE 4.2.8.6. Voici un cas particulier élémentaire où l’on peut calculer l’exposant 
de Lyapunov grace à une chaîne de Markov finie sous jacente. Supposons que la variable 
aléatoire A(0) ne prenne que les deux valeurs : 


€ € 


On remarque tout d’abord que 
x) = D RG +D-x@]+x(0) = D TAG) — xG)1+x(0) 
i<n—]1 i<n—]1 


D GG), AG) +x(0) (4.39) 


i<n—1 


[-Y 
Ke] 
h 


où f(t, D) = Dit — f est invariante par addition d’un scalaire sur le premier argument, soit 
VLeR, fA@ft,D)=f{A+1,, D)= f(,D). (4.40) 


Introduisons l’opérateur de normalisation 


Nu Rue [5 eu | 


X2 — X] 


et définissons l’état normalisé z(i) = N(x(i)). D’après (4.40), on a 
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«on _ 1 
DIE 


1 
D f GG), AG) + 2x0) (441) 


i<n—1 


ce qui exprime x(n)/n comme une moyenne ergodique qui ne dépend que de z et À, à l’in- 
fluence près de la condition initiale qui est négligeable (terme Lx (0)). Il est donc loisible de 
choisir x(0), et nous prenons l’état initial (déterministe) x(0) = [e e]!. Grace à la norma- 
lisation, z est une chaîne de Markov n’ayant que les deux états Zo = W(x(0)) = [e ef et 
Zi = N(Bz(0)) = [e — 1]7. Le graphe ainsi que la matrice M de cette chaîne de Mar- 
kov sont représentés sur la Figure 4.7. L’unique mesure invariante est x = [1 — p p], donc 


F0 — 4 AP 


FIG. 4.7: Chaîne de Markov après Normalisation 


d’après le théorème ergodique pour les chaînes de Markov et vu la formule (4.41) 


€=imElx(n)/n1 = Ef(o, A(0))xo + Ef(z1, A(O)}ri 
bXx1+gx2)(1—p)+(p x 1+gq x2)p = p +2q 


Notes et références bibliographiques 


Les résultats de valeurs propres et vecteurs propres sont l’un des points centraux de la 
théorie que les différentes écoles (max,+) ont prouvé ou reprouvé avec des motivations et des 
degrés de précision variés. Les premiers résultats semblent remonter à Cuninghame-Green, 
Vorobyev [68], Romanovskÿ[62], Gondran et Minoux [42, 43, 46]. Le théorème de cyclicité 
dans la forme présentée ici est du à Cohen, Dubois Quadrat et Viot [23], repris dans [4]. Un 
théorème de cyclicité valide dans des dioïdes plus généraux (loi multiplicative simplifiable 
et condition de type archimédien) a été obtenu par Dudnikov et Samborski[32, 52]. Le lec- 
teur pourra aussi consulter le livre de Chretienne et Carlier ([18],Chapitre V, “Méthodes de 
Circuits Critiques”) où les questions de périodicité sont traitées du point de vue de la théorie 
des graphes bivalués. Les graphes d'événements stochastiques sont traités en détail dans la 
dernière partie du livre [4]. 


Chapitre 5 


Automates et Systèmes à Événements 
Discrets 


Introduction 


On a considéré dans l’introduction de ce cours deux classes de systèmes. 


(a) Les systèmes déterministes présentant des phénomènes de synchronisation et de sa- 
turation, et dont le comportement temporel se représente par des équations linéaires 
récurrentes sur un dioïde réel (algèbre (max, +) ou (min, +)) Voir par exemple le li- 
mitateur de débit, les processus d’assemblage ..…. 


(b) Des systèmes plus généraux présentant en sus des phénomènes de concurrence (par 
exemple, deux utilisateurs partageant une ressource). Ces systèmes ne peuvent plus être 
considérés comme “déterministes”, dans la mesure où plusieurs événements peuvent 
se produire (plusieurs trajectoires partent d’un même état initial). Cependant, on a 
vu sur des cas particuliers que l’ensemble des trajectoires possibles était déterminé 
comme l’ensemble des chemins d’un graphe, et était solution d’équations linéaires sur 
un dioïde de langages. 


Ce chapitre d’intéresse aux systèmes de type (b). L'outil mathématique privilégié pour traiter 
ces systèmes d’un point de vue “logique” est la théorie des langages formels et des automates. 
Disons grossièrement qu’on se donne un alphabet représentant les différents événements et 
que le comportement d’un système à événements discrets sera représenté par le langage 
reconnu par un automate. On s’attaquera tout particulièrement aux problèmes de contrôle : on 
veut garantir un “bon comportement” de l’automate en interdisant éventuellement certaines 
transitions. Il s’agit là principalement de résultats relatifs au comportement logique de ces 
systèmes. Ce chapitre est conçu comme une initiation à cette théorie, et nous n’avons pas 
cherché à inclure les nombreuses et intéressantes extensions. Nous nous en sommes tenus 
aux résultats de base qui recèlent à notre avis l’essentiel des idées. Voir les notes à la fin pour 
un bref survol des prolongements. 


5.1 Langages et Systèmes à Événements Discrets 


5.1.1 Spécification d’un Système à Événements Discrets par un langage 


Soit Z un alphabet fini. Chaque lettre de Y représente un événement. Une suite 
d'événements correspond à un mot w = a1...4, € X*, ce qui se lit : d’abord l’événement 
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ai puis l'événement a, .... Le mot vide (unité pour le produit de concaténation) sera noté 
e. On rappelle que u est un préfixe de w s’il existe v tel que uv = w, ce que l’on note u < w. 
L'ordre < est qualifié de préfixiel. Une partie X € X* est dite préfixielle si 


weXetu<w=Ruex, 
(X est stable par préfixe). Etant donné une partie X € Y*, on note 
X={ueXZ*| 1wEX, u <w} 


la clôture préfixielle de X. On a évidemment X = X, X > X, X = X & X est une partie 
préfixielle, ce qui justifie le terme “clôture”!. On spécifie le comportement d’un SED par 
une partie L € X* (on dit aussi un langage). Si w = a,...a, € L, cela signifie que la suite 
d'événements a, ...a, est compatible avec le fonctionnement du système. Il est clair que tout 
préfixe de w est également admissible. On spécifiera donc un système à événement discrets 
par une partie préfixielle L de Y*. 
EXEMPLE 5.1.1.1. Soit un poinçonneur automatique de tickets de métro. Soit l’alphabet 
Z = {a, b, c}. On associe aux lettres a, b, c les événements suivants : 

a: ticket introduit 

b: ticket accepté 

c: ticket refusé 


On se convainc que l’ensemble des trajectoires possible est représenté par le langage suivant : 


L = (ab@ac}* = ab @ac @ (ab)? @& (ac)? & abac ® acab & 
— a@ab@ac@aba & aca & abab & ..…. 
— (e&a)(ab& ac) (5.1) 


Selon l’usage, on a identifié dans ce qui précède la lettre x et le singleton {x} € X*,et l’on a 
noté l’union &. Ainsi, a @ b désigne le langage {a} U{b} — {a, b} € Y*. La notation additive 
est justifiée par le fait que (P(Z*), U, .) constitue un dioïde (c’est le dioïde des langages sur 
ZX). On peut déjà noter au vu de (5.1) que le langage décrivant le système n’est pas quel- 
conque, mais constitue une partie rationnelle de X*, 1.e. que L s’écrit comme une expression 
finie faisant intervenir des parties finies et les opération @, @, x (cf. Chap. 2,2.2.4.7). 


Notre problème est ici de manier effectivement (1.e. d’un point de vue calculatoire) des 
langages a priori infinis. On peut dans certains cas comme ci-dessus décrire un langage infini 
à l’aide d’une expression rationnelle. On peut aussi le décrire comme l’ensemble des mots 
reconnus par un automate. Nous précisons ce point de vue dans le numéro suivant. 


5.1.2 Rappels sur les automates 


DÉFINITION 5.1.2.1. On appelle automate (déterministe) sur l’alphabet Y un quadruplet 


A = (Q, ô, go: Q:) , 


où © est un ensemble (ensemble des états), gp € © un état initial, ô : 5 x Q — Q une 
fonction en générale partielle (1.e. non définie pour certaines valeurs) qualifiée de fonction de 
transition, et Q, € Q un ensemble d’états dits finaux ou encore appelés éfats d’acceptation. 


lEn théorie de la résiduation et des correspondances de Galois [12], une cloture est une application croissante 
f d’un ensemble ordonné dans lui même telle que fof = f, f > Id. 
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Dans la suite, les automates ne seront pas supposés finis, sauf spécialement lorsqu'il 
s’agira de calcul effectif. 

La manière la plus simple de spécifier un automate est graphique : l’ensemble Q des 
états est représenté par l’ensemble des sommets d’un graphe orienté. On dessine une flèche 
qg — gq' valuée par la lettre a € Y ssi ô(a,q) = q’. La fonction de transition 6 peut être 
donnée sous forme de tableau. L’état initial est repéré par une flèche entrante. Les états finaux 
par des flèches sortantes. Par exemple, prenons Z = {a, b}, O = {1,2}, go = 1 et O4, = {2}. 
Soit la fonction de transition & donnée par le tableau à droite de la Figure 5.1 : par exemple, 


PEN 1 2 
L ) D 1 a 
Ÿ 2116 


fonction de transition 


FIG. 5.1: Représentation graphique d’un automate et représentation matricielle 


la première ligne signifie que ô(a, 1) = 2 et que 6(b, 1) n’est pas défini. On obtient alors le 
graphe à gauche de la Figure. 
On étend la fonction de transition en une fonction partielle sur Z* x Q en posant : 


ô(e,qg)= 4, ô(wo,q) = (0, ô(w, q)) 


(rappelons que e désigne le mot vide). Par exemple, pour l’automate 5.1, 8(ba,2) — 
ô(a, ô(b,2)) — &6(a,1) = 2. Dans la suite, on écrira ô(w, q)! pour “ô(w, q) est définie”. 
Graphiquement, si w — a;...a,, Ô(w, g) est défini ssi il existe un chemin dans le graphe 
partant de g et portant les labels successifs aj...a,. Si ce chemin conduit au sommet 
g', évidemment, ô(w,g) — q’. On dira que le mot w est reconnu par l’automate À si 
ô(w, go) € Q,. Nous noterons 


L(A) = {w € E* | 8(w, go) € Qu} (5.2) 


le langage reconnu par l’automate (1.e. l’ensemble des mots envoyant l’état initial sur un état 
final). 

Du point de vue des systèmes à événements discrets, un automate représente un système 
dynamique, et l’on peut s’interesser à toutes ses trajectoires, y compris celles qui n’arrivent 
pas dans un état final. Soit donc A = {Q, 6, go, Q} l’automate obtenu à partir de À en rendant 
tous les états finaux. On s’intéressera particulièrement au langage suivant 

def 


L(A) = L(A) = {w € E* | 8(w, d)!} (5.3) 


que nous qualifierons de comportement de l’automate? A. On a clairement 


LA) c LA), LA c L(A) = L(A) . (5.4) 


EXEMPLE 5.1.2.2. Pour l’automate de la Figure 5.1, un mot w est reconnu ssi il conduit de 
l’état 1 à l’état 2, ï.e. ssi il commence par un nombre quelconque de circuits ab et s’achève 
par une lettre a, soit 


2Cette définition est locale à la communauté des systèmes à événements discrets. Toutes les autres définitions 
sont standard en théorie des langages. 
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L(A) = a & aba & ababa & ...— (ab)*a . 


De même, £(A) désigne l’ensemble des chemins partant de 1 et arrivant dans un état quel- 
conque, soit 


L(A) = e @ a ® ab  aba ® abab  ... = (ab)*(e E a) 
On définit les ensembles des états respectivement accessibles et co-accessibles par 


Qx = {(q € Q | Aw € E*, 6(w, qo) = 4}, 
Qco = {g E QO | IwE EX”, ô(w,q) € Qu} . 


L’automate est accessible ssi tous ses états le sont. Définition analogue pour la co- 
accessibilité. Par exemple, l’automate de la Figure 5.1 est accessible et co-accessible. L’au- 
tomate est dit émondé si O = Qx = Ov. Soit 


OQëm = Qc N Qc; Êem = Ê|S x Qem . 


On définit l’automate 


Em(4) = (Qëm; Oem: 40: Qm N Oem) 


qui est clairement émondé et reconnaît le même langage que À. On a trivialement : 
FAIT 5.1.2.3. Si À est co-accessible, alors L(A) = L(A). 


Autrement dit, sous l’hypothèse de co-accessibilité, £(4) n’est autre que l’ensemble des 
mots qui se prolongent en des mots reconnus. 


EXEMPLE 5.1.2.4. L’automate de la Figure 5.1 est émondé. Pour l’automate À de la Figure 
5.2,ona 


Où = {1, 2, 2}; Oo = tl 2; 4} 
On a représenté à droite de la Figure l’automate émondé Em{4). On voit facilement que 


L(A) = ba* = LŒEm(4)), £L(A) = e ® ba* & aa* £ LA = L(Em(A)) = e P ba* . 


CRE N a 
: œ 3 ) / A 
A ae OS 
DS 
A 22 Fr |Em(4) 2} 
Fr A | A ) 


FIG. 5.2: Un automate et l’automate émondé associé 


EXEMPLE 5.1.2.5. On a représenté sur la Figure 5.3 l’automate reconnaissant le langage du 
poinçonneur de métro. 


Le lecteur pourrait se demander quel est le lien entre les langages rationnels (cf. Chap 
2,2.2.4.7) et les langages reconnus par des automates. On dira qu’un langage est régulier s’il 
est reconnu par un automate fini (1.e. avec un nombre fini d’états). Nous citons pour mémoire 
le plus célèbre des théorèmes de la théorie des langages formels : 
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S GE) L(A) = (ab & c))*a 


c L(A) = LA) = (a(b @ c))*(e & a) 
FIG. 5.3: J'suis l’poinçonneur des Lilas (cf. 5.1.1.1) 


THÉORÈME 5.1.2.6 (KLEENE). Soit Z un alphabet fini. Un langage est régulier ssi il est 
rationnel. 


Nous ne prouverons pas complètement ce théorème. Notons simplement que l’implica- 
tion régulier — rationnel a déjà été prouvé au chapitre 2 lorsque l’on a montré que l’en- 
semble des chemins entre deux sommets d’un graphe (fini) est une partie rationnelle. Le 
langage reconnu par l’automate n’est autre que l’ensemble des chemins du sommet initial 
aux sommets finaux. Pour la réciproque, il suffit de vérifier que l’ensemble des langages 
réguliers est stable par union, produit et étoile, ce qui revient à fabriquer des automates finis 


reconnaissant À @ B, AB et A* à partir d’automates reconnaissant A et B. Nous renvoyons 
le lecteur à [33] pour la preuve. 


5.1.3 Supervision d’un automate 


On considère une partition de l’ensemble des événements : 
Z=ZnUX» . 


Xin désigne l’ensemble des événements ‘“incontrôlables” (une machine tombe en panne, un 
client arrive, ... ) ZX désigne l’ensemble des événements que l’on peut interdire ou bien 
autoriser. Par exemple, si l’événement a € Y4 représente l’entrée d’un train sur une portion 
de voie, on peut interdire l’événement a au moyen d’un feu rouge. 


DÉFINITION 5.1.3.1 (SUPERVISEUR). On appelle superviseur une application 
f : L(A) = PE) 


telle que Vw € L(A), f(w) > Ein. 


f(w) s’interprète comme l’ensemble des événements autorisés connaissant les 
événements passés donnés par le mot w. On impose f(w) D Yi (les événements in- 
contrôlables sont autorisés de manière permanente). 


DÉFINITION 5.1.3.2. On appelle comportement d’un automate À soumis au superviseur f 
le langage préfixiel £ ;(.4) de X* défini récursivement comme suit : 


e € L(A) ue 2. _ (5 5) 
VweY*,aeY, (wa € L(A) & w € Lr(A), a € f(w)et wa € L(A)) : ‘ 
Le langage contrôlé par f est par définition : 


Lr(A) = LH(ANL(A) , 


1.e. l’ensemble des mots du comportement supervisé qui aboutissent dans un état final. L’au- 
tomate supervisé est dit non bloguante si L (A) = L(A). 
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La condition de non blocage exprime que toute tâche commencée par le système super- 
visé, i.e. tout mot de £; (À) doit pouvoir se completer en une tâche achevée (mot de £(4)). 

Pour rendre effective la définition ci-dessus d’un superviseur, on peut introduire un au- 
tomate “contrôleur” réalisant cette supervision. Soit un automate S — (X,Ë,xo, X4), et 
p : X — P(E) une application telle que Vx € X, (x) > Xi. On dira que S réalise 
le superviseur f si pour tout mot w € L (À) : 


F(w) = pE(w, xo)) . 


Autrement dit, on fait lire au contrôleur $ le même mot que l’automate à contrôler À, et l’en- 
semble f (w) des événements autorisés est un feedback de l’état du contrôleur S. De manière 
plus précise, considérons l’ensemble des états (resp. l’état initial, resp. les états finaux) X x Q 
(resp. (xo, go), resp. (X4, Q,)). Soit la fonction de transition £ x 6 définie par 


(E(G(x)),8(0,q)) sio € b(x) 


X ô(©, (x, —= 
$ (6, ,4)) non défini sinon. 


On notera S /A l’automate représentant “A supervisé par S”, soit 


SJ/A = Ac(X x Q,E x à, (x0, 40), (Xa: Qa)) , 


où Ac dénote la partie accessible d’un automate. Il résulte du fait 5.1.2.3 que l’automate 
supervisé est non bloquant ssi S/A est co-accessible. 


EXEMPLE 5.1.3.3. Soit l’automate de la Figure 5.4,(a). Considérons le superviseur S de la 
Figure 5.4,(b). On a représenté l’automate $/A en (c). On a 


L(A) = (@y*B)", L(A) = L(À) . 
L (A) = L(S/A) = (œB)*(e D ay”) . 


Lr(A = (aB) . 


D(x2) = (a. 8. y) 


Aa ND) CN dei oo hi à 
Be + | ts DAS" NPA) À A) y \D_ | 
; B 
E = {a. 8,7} | Gi) ={0,8) 
Zi = {a,y}, 2e = {8} : P(x3) = {a, BY : 
(a) Automate À (b) Contrôleur S (c) Automate supervisé S/A 


FIG. 5.4: Un automate supervisé 


5.2 Synthèse d’un superviseur 


5.2.1 Contrôlabilité 


Etant donné un automate À et un langage K, on se demande à quelle condition il existe 
un superviseur f tel que £(.A) = K. On introduit pour cela la notion suivante. 
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DÉFINITION 5.2.1.1 (CONTRÔLABILITÉ). Le langage K C X* est dit contrôlable si 
KREALANCE : (5.6) 


Autrement dit, tout préfixe de K suivi d’un événement incontrôlable admissible pour 
le comportement de À reste préfixe de K, où “les événements incontrollables ne sont pas 
illégaux”. La proposition fondamentale qui suit justifie le terme de “‘contrôlable”. 


PROPOSITION 5.2.1.2. Soit À un automate. Pour tout K € £(A) non vide, les deux asser- 
tions suivantes sont équivalentes : 


1. il existe un superviseur f tel que L À = K 

2. K est une partie préfixielle contrôlable. 
Démonstration. (i}={ii) : Résulte évidemment de la définition de £ (A) (cf. (5.5). 
Gi)=(). Soit B — (X,Ë6,x0, X) un automate émondé reconnaissant K (comme K est 


préfixiel, on peut prendre tous les états finaux). Définissons, pour x € X, 


{o | 1€ XK, E(s, xo) = x, so € L(A)etso 6 K]} 
{o | se K, E(s,xo) = x etso € K} 


2, 
ZE, 


De manière intuitive, Z° représente les événements qu’il faut prohiber (car il sortent de la 
cible K) et Y! est l’ensemble des événements que l’on peut autoriser. 

Comme l’automate est émondé, l’état x est accessible depuis xo, donc la condition 
35, Ë(s, xo) — x est toujours remplie. Notons tout d’abord que E° NY! = ÿ. Supposons 
en effet o € 2° N EX! avec 


se K,E(s°, xo) = x, so € K 
s'eK,E(s',xo) =x, so €eK. 


Alors, £(o,x) = E(o, E(s°, xo)) = E(s°o, xo) qui n’est pas défini (car so & K), mais par 
ailleurs, £(o,x) = E(o, E(s!,xo)) = E(s!o, xo) qui est bien défini : contradiction. K étant 
contrôlable, on a 5h € 24. Soit p : À — Y une application vérifiant D(x) N x0 — ÿ, 
p(x) D 5: U Yjn. On définit le superviseur f par f(w) = p(Ë(w, xo)). On va montrer que 
£; (À) = K. Clairement, 


L(A) CK. (5.7) 
Posons pour un langage M : 
MŸ=MNEV 2{weM| lwl=)j} 


(ensemble des mots de longueur j). On montre par récurrence sur j que Fa f (4) = KO), 
Cas j=l:sioekK,onaoe Es donc o € £(A). Supposons par récurrence que 


ZT} (4 = K%. 


Soit so € KŸ+D avec s € KU) = ras) Comme x — Ë(s, x0) est défini, so € K entraîne que 


—_— —(j +1 : o . 2 
o € X}, donc o € f(s), d’où so € LF(A) et L > (A) D KÜ+, L'autre inclusion étant 


claire d’après (5.7), le résultat est acquis. Ü 
Le problème de supervision pour le langage reconnu £ se résoud de manière analogue. 
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PROPOSITION 5.2.1.3. Soit À un automate non-bloquant (i.e. co-accessible). Pour tout 
K C £L{(AÀ) non vide, les deux assertions suivantes sont équivalentes : 


1. il existe un superviseur f tel que l’autamate supervisé soit non bloquant et L;(A) = 
K 


2. K est contrôlable et K N L(A) = K. 


La preuve, similaire, est laissée au lecteur. 


EXEMPLE 5.2.1.4. Considérons l’automate représentant un magasin de stockage d’une ca- 
pacité de trois unités à gauche de la Figure 5.5 (cf. 1,1.1). Supposons que les deux dernières 
places de stock soient indisponibles, et que l’on veuille limiter la capacité du magasin à une 
unité. Le langage K admissible est alors la partie préfixielle reconnue par l’automate B de la 
Figure 5.5, soit 


K = (ab)*(e&a) =K . 
Supposons l’événement a contrôlable et b incontrôlable : 
En = (0), 2 = {a}. 
On a 
K3; NL(A) = (ab)*(e @ a)bN L(A) . 
On a (ab)*ab C K et (ab)*b N L(A) = ÿ, donc 
K3 NL(A)CK 


ce qui montre que le langage K est contrôlable. Prenons comme superviseur l’automate B — 
({0’, 1°}, &, 0, {0’, 1°}), avec 


P(0") = {a,b}, (1) = {b} 
et posons f(w) — p(E(w, xo)) On a alors £L(A) — K . Supposons maintenant Xi, — 


a a a ! a 


OS HS D D PT 
Ne. . PSS RE, ” LNCA » LR D AN | Fr. 
. ee A7 | RAS 
À B K = L(B) = (ab)*(e @ a) 


Ein = {b}, ce = {b} 


FIG. 5.5: Le langage K = £L(B) est contrôlable 
{a}, Xe = {b}. On a alors 
K Ein N L(A) = (ab)*(e & a)a N L(A) > {4°} 


et a & K, ce qui montre que le système n’est plus contrôlable. Tout cela est intuitivement 
très clair : l’événement a revient à stocker une pièce supplémentaire. Si l’on ne contrôle pas 
a, on ne peut pas limiter la capacité du stock. 
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5.2.2 Sous-langage contrôlable maximal 


Soit K € £(A). Même si K n’est pas contrôlable, il existe une “approximation 
contrôlable” naturelle de K. Il s’agit du sous langage contrôlable maximal, dont nous mon- 
trons ci après l’existence. Soit C(K) l’ensemble des parties de K contrôlables, 1.e. 


C(K)}=UCK| JEnNL(A) CI}. (5.8) 


Comme Ÿ est contrôlable, C(K) est non vide. Au vu de la condition (5.6), il est immédiat 
que cet ensemble est stable par union (éventuellement infinie). Donc il admet un plus grand 
élément. On notera 


KT = sup C(K) = max C(K) 


ce plus grand élément (sous-langage contrôlable maximal). II faut bien voir que la définition 
de KT comme borne sup ne permet pas de le calculer. On va donner ci-après une ca- 
ractérisation alternative de K° comme point fixe d’un opérateur, ce qui permettra de le 
calculer par une technique élémentaire d’itération de point fixe. Définissons l’application 
Q : P(E*) — P(E*) par 


Q(J) = sup{T CK| TEnNL(A) CI}. 
On obtient facilement les expressions équivalentes de <2 : 
Q(J)={HtERK | 1En NL(A) CJ}= K Nsup{T CE*| TEnNL(A) CI}. (59) 
PROPOSITION 5.2.2.1. On a KT = Q(K°). En outre, 
J=8(J) = Kt35J. 
En d’autres termes, KT est le plus grand point fixe de l’opérateur €. 


Démonstration. . En fait, il s’agit d’un lemme tout à fait général sur les structures ordonnées. 
Soit € un treillis complet, et soient f, g : £ — € deux applications croissantes, f préservant 
les sup quelconques : 


f(V/ u) = \ f(u) pourtoutuU CE (5.10) 


ueU ueU 
(4 éventuellement infinie). Définissons de manière plus générale : 
C(K)=UCK| FU) <g(J)} KT supC(K) , 
ainsi que 
C(K,J)={MCK| FM) <g(J)} ,œ(K, J) © supC(K, J) . 


Il suffit clairement de montrer que 


J=o(K,J) = J<Kt (5.11) 
KT = «(K,K') (5.12) 


Il résulte de (5.10) que K° e C(K, KT), donc 
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K'<o(K,K'). (5.13) 
D'autre part, en spécialisant w(K, J) e C(K, J) à J = K, il vient 
F@UK, KT) < g(K°) < g(œ(K, K°)) 


(par (5.13) et par croissance de g). Ainsi, o(K, KT) e C(K), d’où l’on tire o(K, K?) < KT. 
On a montré (5.12). Enfin, J = w(K, J) entraîne J e C(K), d’où J < K1. [] 

A partir de la Proposition 5.2.2.1, il est naturel d’appliquer un argument de point fixe et 
de définir la suite : 


Ko = K, Kj41 = Q(K;). (5.14) 
En vertu de la décroissance de la suite {K,}, on a l’existence de 


Ks= mr: )k 
n— 00 n>0 
et l’on se demande bien sûr si K,, est le point fixe requis. On a K° = (KT) € K = Ko. 
Donc par récurrence KT = Q(K7) € K; et donc 


AN D ES (415) 


jz0 


En général, K4 n’a aucune raison d’être un point fixe et cette inclusion est stricte. Nous 
donnons une condition suffisante dans la section suivante. 


5.2.3 Synthèse d’un automate reconnaissant K ? 


Nous montrons que la convergence de la suite K, vers le sous-langage contrôlable maxi- 
mal est acquise sous les hypothèses suivantes 


5.2.3.1 (Hypothèses). 
1. À est fini et K est un langage régulier (i.e. reconnu par un automate fini) 
2. KC L(A) 
3. L(A) = L(A) . 


L'hypothèse (1) est la plus simple qui permette de manipuler effectivement des langages. 
L'hypothèse (ii) signifie que le langage “cible” K est une restriction du comportement ad- 
missible du système £(A), cas auquel il est naturel de ce ramener au besoin en remplaçant 
K par K N £L(A). La condition (iii) enfin signifie qu’on se préoccupe de tous les événements 
(et pas seulement des mots acceptés). Cette condition est aussi souvent réalisée. Le lecteur 
pourra se reporter aux exemples de la section suivante pour juger de la classe de problème 
prise en compte. Moyennant ces hypothèses, on va montrer que chacun des langages K'; de 
la suite définie ci-dessus et convergeant vers K+ est régulier. On va construire pour cela à 
chaque étape un automate C; reconnaissant K';, et donner la règle produisant C;,, a partir de 
C;. On montrera que la suite C; converge en un nombre fini d’itérations vers un automate C 
reconnaissant K = K'. 

Dans la suite, on notera 


Z4(g)={ae >| &8(a,q)!} 


(c’est l’ensemble des événements possibles dans l’état g). Soient 
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A = (Q,5, qo, Qa) et B = (X,8, x0, Xa) 
deux automates émondés tels que L(B) € £(A). On dit que B est un raffinement de A si 
Vs,1 € L(B), Es, xo) = (6, xo) = 6(s. go) = (E, go) : 
On a alors une unique application h : X + Q telle que 
Vs € L(B), h(Ë(s, xo)) = ô(s. do) : 


Par exemple, l’automate de la Figure 5.6,(b) raffine celui de la Figure 5.6,(a), la fonction h 
étant tabulée sur la Figure. L'intérêt des l’hypothèses 5.2.3.1 est que l’on peut fabriquer un 


G 
à Fi oc h(x1) = qi 
mn Re ue h(x2) = g2 
HS. y T Mende = = 
“gi MH, CEE ne PR = 
Ts | . Te dl 4) — 43 
{ya 
NEZ 


(a) (b) 


FIG. 5.6: Un automate et un raffinement 


automate raffinant À et reconnaissant K 


LEMME 5.2.3.2. Moyennant l'hypothèse 5.2.3.1, il existe un automate fini Co raffinant À tel 
que L{Co) = K. 


Démonstration. K' étant régulier, il existe un automate fini B — (Z,@,z0, Za) tel que 
£(B) = K. On définit le produit À x B des automates À et B par 


A x B=(Q XZ,8 x @, (go, zo), Qa X Za) 


(Ô x æ)(6,q,1) = (ô(o,q), a (o,2)) . 


Il est clair que À x B reconnaît l’intersection des langages reconnus par À et par B. Comme 
K C£L(A),ona 


L(A x B)=K . 
On ne change pas le langage reconnu par un automate en l’émondant. En posant 
Co = Em(A x B) , 
on a un automate reconnaissant K et raffinant 4. En effet, trivialement 
ô X œ(s, go, zo) = à X @(f, qo, zo) — Ô(s, go) = Ô(F, qo) 


[] 
On considère à nouveau la suite de langages {K';} définie plus haut. 
LEMME 5.2.3.3. Soit À un automate fini émondé et C; un automate fini émondé raffinant À 
tel que K; = L(C;). Alors w € Q(K;) ssi w € K'; et 


Vu ew, x = E(u, xo) = 2 (h(x)) N Ein € 2G (x) . 
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Démonstration. Les propriétés suivantes sont équivalentes 


w € Q(K;) 

w € K; et WEin N L(A) CK; 

we K; et Vu < w, uXi N L(A) C K; 

w € K'jet Vu < w, Vo € Ein, uo € L(A) = uo € K; 
w € KjetVu <w, Vo € ZA(8(u, qo)) N Ein, uO € K; 
w € K;et Vu <w, Z4(8(u, qo)) N Ein € 24 (E(u, Xo)) 


d’où la conclusion. [] 
Ainsi, on peut construire l’automate C;,, reconnaissant K;,1 = &(K';) en supprimant 
les états de €; qui violent la contrainte suivante : 


2x) N En € ZA (x) . (5.16) 
Soient précisément 


X'=frex| ZA) NE CE (x)} , 


X : =X,NX, 
et soit la nouvelle fonction de transition 


E(o, x) si £(o, x) € X’ 


non défini sinon. 


E(o, x) = 


On définit alors 


Em(X',8’,x0, X!) sixo € X’ 
cs _ m(X”, &”, xo, X,,) | X0 € (5.17) 
(7 sinon. 
THÉORÈME 5.2.3.4. Soit C;,1 comme en (5.17). Alors 
LCD=rk; = 0). (5.18) 


En outre, il existe un indice i tel que C; = C;:1 et l’on a 
L(C;) = K* . 


Autrement dit, on obtient en un nombre fini d’itérations un automate reconnaissant le 
sous-langage contrôlable maximal. 


Démonstration. du Théorème. L’assertion (5.18) résulte immédiatement du Lemme 5.2.3.3. 
Montrons la convergence. Comme C;,, s’obtient en supprimant des états de C;, la suite C; 
converge en un nombre fini d’itérations (éventuellement vers l’automate vide), soit pour un 
certain i, C; = C;:1 et donc 


Ke LCG)= LE CG=Rks = Ks). 


K est donc point fixe de Q, et comme K? est le point fixe maximal (Proposition 5.2.2.1), 
K € KT. D'autre part, on a observé dans la section précédente que l’on a toujours KT € 
Koo- Cl 
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5.3 Exemples 
Nous donnons quelques exemples d’application de l’algorithme décrit ci-dessus. 


5.3.1 Exemple 1 (cf. [691) 


Soit 2 = {o,@, B}, Ein = {8}. M = (187 B)B*, K = if? ® æf*.Ona 
représenté sur la Figure 5.7 un automate À tel que M = £(A). Le mot «8? appartient à K, 


À Ê 4 L L. 
DB-O)8 ET 8 d-É-2)8 6) 
OS LA 
K ES \2 ce 
< 
\ à n* 
DENT: 4 7 B 
Re A 
L(A) = 02B* @ «1B°B* L(B) = K 


FIG. 5.7: Automates reconnaissant M et K 


mais «8° g K. Comme B est incontrôlable, cela entraîne que «18? # K1. De la sorte, il est 
intuitivement clair que KT = e @ a2B*. On peut le vérifier plus formellement à partir de la 
construction de point fixe décrite plus haut. On a pour K, la table des transitions suivante : 


| [0 1 2 3 4 | EGG)NEn | 2@) | 
(041 (0%) 


© 


E © D 


Pour À, 2(3) N Ein = {B} n’est pas inclus dans Z#(3), il faut donc éliminer l’état 3 du 
tableau de B. On obtient alors le nouvel automate qui reconnaît K' : 


IST 


0 (041 


Ki = 8 @ f* 


En poursuivant, on obtient 


RE — (RGO) 0 Ein 


K) : 0 (0%) ne 
1 
4 


K: = E 28" 


BEMIEUONE" 


USD 6 
: p p 
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K;3=e@of*. 


La contrainte (5.16) est maintenant satisfaite pour tous les états : on a donc KT = e @ a2B*. 


5.3.2 Chat et souris dans un labyrinthe 


On reprend l’automate de la Figure 1.2 du Chapitre 1 représentant les mouvements pos- 
sibles du chat et de la souris. L’automate B représentant le comportement admissible du 
système s’obtient en supprimant tous les états diagonaux (où le chat et la souris se trouvent 
dans la même pièce), soit l’automate de la Figure 5.8. Considérons tout d’abord le cas fort 


(2 KG) 
/ \ €] 2 
|. © 
Re A 2 il 
{ | | ( | 
@X | | @) :@i 
ml L— ———"f 7 mi 
n ai . ci = 

Gi ) @2 \ @) 
tb) K = L(B) (c) Em(K) 


FIG. 5.8: Comportement admissible du chat et de la souris (c) 


simple où Ein = Ÿ : cela signifie que l’on contrôle toutes les trappes. Il suffit alors d’interdire 
les transitions conduisant aux états diagonaux dès qu’on se trouve dans les états les précédant 
immédiatement. Soit £ la fonction de transition de l’automate (c). Cela revient à définir un 
contrôle de la forme 


f(w) = D(E(W, iojo)), iojo = 12 (état initial), 


où le retour d’état o (p(q) décrit l’ensemble des transitions autorisées au point g = ij) est 
donnée par le tableau suivant : 


12 | g(g) = {c1, m2} 

13 | g(q) = {m2, 2} 

21| d(@)=1{m} 

23 | (g)= {im} 

32| p(a)={a} 
Autrement dit, on interdit à la souris tout mouvement qui la ferait rentrer dans une pièce où le 
chat se trouve déjà et de même pour le chat. Un cas plus intéressant est celui où Xi, = {ci}. 
Le chat peut alors de manière incontrôlable emprunter la trappe c,. Calculons le sous-langage 
contrôlable maximal. Pour l’automate de la Figure 5.8, on a le tableau suivant. 


AD 2128-92 ES mONTSEe 2000 | 
12 m) C1 (a 
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La contrainte Z(h(x)) N Zn € XZ(x) est violée pour x — 13. Lorsque l’on supprime cet état 
de l’automate À, on obtient l’automate € à gauche de la Figure 5.9, pour lequel les états 21 
et 23 ne sont plus accessibles. L’automate émondé correspondant est représenté à droite de la 
Figure. On vérifie en dressant un tableau analogue que la contrainte est maintenant satisfaite 
pour tous les états. Cette automate reconnaît donc le sous-langage commandable maximal de 
K. Ce langage est réduit à c?. Autrement dit, le chat peut circuler librement entre les pièces 1 
et 3, ce qui est normal, l’événement c; étant incontrôlable. En outre, la souris est cantonnée 
dans la pièce 2, comme on aurait pu le voir de manière tout à fait élémentaire. 


@) 
rue * A 
/ \ c1 
J_ \a | 
OR US | | 
PA) LRO | | 
1 ) | 23 
Nr es | _ se | | 
SN pr: | | 
c] \ / a \ . { 
32 \ de ) 
DA < te ci 


FIG. 5.9: Calcul de K? 


Commentaires et références 


La théorie de la supervision des automates est due à l’école de Ramadge et Wonham. Cet 
exposé est largement inspiré des trois papiers [58, 59, 69], que le lecteur pourra utilement lire. 
Signalons quelques problèmes et quelques directions d’approfondissement. 

Dans le cas où les hypothèses du théorème 5.2.3.4 ne sont pas vérifiées, on peut se 
demander si la suite K; converge encore vers K1. La convergence n’a pas lieu en général 
pour des langages non réguliers [69]. Par contre, les conditions £(4) — £(4) et K C £(A) 
du Théorème 5.2.3.4 qui sont commodes pour définir la suite d’automates reconnaissant les 
K; peuvent être relaxées comme suit. Lorsque M est un langage régulier, on note |[W|| le 
nombre minimal d’états d’un automate reconnaissant M. 


THÉORÈME 5.3.2.1. (Ramadge et Wonham [69]) Supposons À fini et K régulier La suite 
{K ;} définie en (5.14) converge en un nombre fini d'itérations vers K T. En outre 


KT < ILCAIIKI + 1 


On s’est limité dans ce chapitre au minimum concernant la supervision d’automates. Il 
y à une méthode de nature différente qui permet de calculer le langage controlable maximal : 
au lieu de calculer une suite d’automates, on soustrait à K° le langage des mots indésirables 
que l’on caractérise comme une sorte de résiduel. Nous avons privilégié ici par simplicité la 
première approche, historiquement la plus ancienne et très naturelle, mais il faut bien voir que 
la complexité de l’approche via la résiduation est meilleure. En effet, la première approche 
émonde de manière progressive un automate raffinant À et reconnaissant K (soit au plus 
IL(CA)IIIIK 11 états), donc le nombre d’itérations est borné par le nombre d’arcs à émonder, 
soit au plus (IILCAIIIK1D?. Le second algorithme exposé dans [16] valable pour une partie 
K préfixielle a une complexité de O(||£(A)|||1K1l). 
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Dans le cadre, de cette introduction, nous avons, pour faire bref, laissé de coté les 
problèmes de contrôle en observation partielle. La théorie présentée ici s’étendre de manière 
élégante, pour peu que certaines conditions naturelles d’observabilité ou de clôture du lan- 
gage légal K par rapport aux observations (normalité) soient vérifiées [701]. 

On a ignoré aussi le principal problème pratique qui est lié à la taille des automates : 
lorsque l’on couple un système à n état avec un système à p états, on obtient en général un 
automate produit (à np états). Pour cette raison, les questions de contrôle décentralisé sont 
essentielles. Il s’agit d’implémenter uniquement des contrôleurs locaux sans avoir à manipu- 
ler des produits d’automates. Certaines conditions de décentralisation de nature algébrique 
ont été mises en évidence. Nous renvoyons le lecteur intéressé à [70, 59]. 


Annexe À 


Travaux dirigés : equations linéaires 
dans les dioïdes 


Les exercices de base sont signalés par le signe S. Les exercices d’une difficulté particulière 
sont signalés par le signe 


A.1 Leçon de calcul 


Que valent les expressions suivantes dans l’algèbre (max, +) : 


2 
283=? 2®@-1=7 2®-00=? 0@3=? 37 2=7 ==? V8? 


A.2 Equation du second degré dans l’algèbre (max,+) 
Résoudre 
ax? ® c = bx 


dans (RU {—co}, max, +). 


A.3 Symétrisabilité d’un dioïde 


Si Un élément x d’un dioïde peut il avoir un symétrique (i.e. un élément x’ t.q. x ® x’ = 
€) ? Un dioïde est-il plongeable dans un anneau ? 


A4 Étoiles de séries 


Exercices 2.2.1.8 et 2.2.1.9 du Chapitre 2. 


A5 Étoile de matrice 


Ÿ Calculer l'étoile de la matrice suivante dans le dioïde Rain - 
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£ 2 3 
A=| —1 € —2 
5 3 2 


Que dire de l’étoile de A dans l’algèbre (max, +) ? 


A.6 Formule du binome dans les dioïdes et applications 


a) Dans un dioïde dont le produit est commutatif, donner une “formule du binome” pour 
(a & b)". Que devient la formule du binome dans le dioïde Rx ? Y a t-il une formule du 
binome dans un doïde non-commutatif ? 


b) ({ sans l'indication à la fin de l’exo.) Si en outre le produit est simplifiable, montrer que 
l’on a 


(ab) =a"@b". 
c) Montrer que sous les mêmes hypothèses 
ab<a @b . 


d) On appelle groupe réticulé un groupe (G, -) muni d’un ordre de treillis tel que Va, b, c : 
(a V bjc = ac V bc, c(a V b) = ca V cb. d.1) Montrer que ces identités sont équivalentes 
à Va, b,c: (a A bjc = ac A bc, c(a À b) = ca A cb. d.2) Quelle est la différence entre un 
dioïde et un groupe réticulé ? 

e) Montrer que dans un groupe réticulé commutatif, l’équation x? = a (p e N, a € G fixés) 
admet au plus une solution (on pourra montrer que x? = e entraîne x — e). 

f) Trouver tous les groupes réticulés finis. 

g) Exhiber un groupe réticulé où l’équation x? — a admet plusieurs solutions. 


[Indication pour b) Comparer (a & b)?(a & b) et (a?  b?)(a &b).] 


A.7 Equation implicite dans un dioïde de langages 


On considère (P({a, b}*), LU, .) le dioïde des langages sur {a, b}. Trouver la plus petite 
solution (pour l’inclusion) de l’équation 


x =a(aæb}x Pb. 
La solution est elle unique ? Indication : pour l’unicité, si on n’y arrive pas à la main, on peut 
recourir à l’argument massue exposé dans l’exercice suivant. 
A.8 Un peu de topologie 
On considère À un alphabet fini, et étant donnés deux langages L, L’, on pose 
œo(L,L")=inffneN]| LANA"Z£L'NA"} 


(avec inf Ÿ — +oco). Soit o : 0 < po < 1. Montrer que l’on définit une distance d sur P(A*) 
en posant 


A.9 Énumération d’un ensemble de chemins 103 


d(L,L')=p%%22), 
On pourra montrer un peu mieux, à savoir que d vérifie l’inégalité ultramétrique : 
d(L, L”) < max(d(L, L'),d(L', L")) . 
On considère une équation implicite de la forme 
x=Lx@L (A.1) 


où L, L' sont des langages, et L ne contient pas le mot vide. Montrer que l’application x + 
Lx & L' est contractante. Que peut on en conclure pour les solutions de l’équation (A.1). 


A.9 EÉnumération d’un ensemble de chemins 


Soit le graphe à gauche de la Figure suivante. 


& LYT 3 EU". # 
1 ) ; —# } { 12) —% ) 
Nu LD a: D DAT 
Le ù A —20 7 
su DPE LA 


circuits(2)=(c @ aby* circuits(4)= ? 


(a) (b) (c) 


L'ensemble des chemins allant du sommet 2 au sommet 2 est égal à 
(ab & cc} = ab@c@ abc @ cab &abab &... 


(on lit les mots de droite à gauche). 


a/ Donner une expression analogue pour l’ensemble des chemins 4 — 4 dans le graphe (b) 
de cette Figure. 


b/ Donner l’ensemble de tous les circuits élémentaires passant par le sommet 4. 


c/ On considère maintenant le graphe de la Figure (c). Quel est le poids maximal d’un circuit 
passant par 4 ? Le poids moyen maximal d’un tel circuit ? Le poids maximal d’un chemin de 
4 à 27 Le poids minimal d’un tel chemin ? 


A.10 Déterminants et formules de Cramer 


Soit D un semi-anneau, À une n x n-matrice à coefficients dans D. On définit le 
déterminant positif : 


det* A — & S Àio (i) 


o permutation paire à 


et det_ À de manière analogue pour les permutations impaires, de sorte que l’on a “formelle- 
ment” det A = det* À — det A. On désigne par A(;, la matrice A privée de la ligne à et de 
la colonne j. On définit la comatrice positive : 
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det* A, sii—+ j pair 
(com* A); = Ho ms 
det Ajj) Sii + j impair 


et la comatrice négative com A en inversant les conditions de parité. 
a/ Montrer que 
(com* A)? A @ (det” A) Id = (com A)7 À @ (det* A) Id 
b/ En déduire le résultat suivant (condition de Cramer). 
Ax =b— (det*A)x ® (com A)7b = (det A)x @ (com*A)/b 


c/ Application. Résoudre dans Rx le système de deux équations linéaires suivant : 


HIHDE 


Annexe B 


Travaux dirigés : Réseaux de Petri et 
Graphes d’Evénements 


B.1 Philosophes Dînant 


Exercice 3.1.1.7 du Chapitre 3. 


B.2 Cellules de production en série 


On considère la ligne de production formée de deux cellules représentée sur la Figure 
B.1. Des pièces brutes arrivent dans un stock de capacité illimitée S;. La cellule C; comprend 
1 machine, qui met 1 unité de temps pour traiter une pièce. A l’issue du traitement, la pièce 
est déposée dans un stock S; situé entre C; et C2. La cellule C est formée de deux machines 
pouvant travailler en parallèle (1.e. traitant des pièces de manière indépendante) et mettant 
chacune deux unités de temps pour traiter les pièces. À l’issue, les pièces sont déposées en 
S3. 


cellule Ci cellule C> 
stock entrée une machine deux machines 
lut 2ut 


FIG. B.1: Ligne de production 
a/ Dessiner un graphe d’événements temporisé représentant ce système. 
b/ Ecrire un système d’équations de la forme 
x(#) = Aox(f) P Aix(t — 1)... Agx(t — k)  Bu(t) 


dans l’algèbre (min, +). Que dire du système analogue d’équations dans l’algèbre (max, +) ? 
Ecrire une représentation sous forme d’état 


xt) = Axt—1)@ But), y) =Cx(t). 
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dans l’algèbre (min, +) en essayant de minimiser la taille de l’état. 
c/ Calculer la série de transfert du système. 


d/ On suppose que le système commence à fonctionner à l’instant 0. A cette date, une très 
grande quantité de pièces brutes arrivent dans le stock S,. Représenter le nombre de pièces 
achevées en fonction du temps { > y(f) (compteur de sortie). Quel est le taux de production ? 
Que peut on dire du stock entre C et C; ? Que se passerait-t-il si les machines en C prenaient 
trois unités de temps pour traiter les pièces ? 


B.3 Simplification de Réseau de Petri 


On considère le graphe d'événements de la Figure B.2, 


a/ Ecrire un système d’équations de la forme 
Xn = AXn_1 ® Bus, n° Cxn (B.1) 


pour les dateurs dans l’algèbre (max, +). 
b/ Calculer la série de transfert du système. 


b’/ Quel est le débit maximal admissible pour chaque transition (nombre de tirs par unité de 
temps) ? 


c/ Peut on remplacer À, B, C par des matrices 4’, B’, C’ de dimension plus petite telles que 
la relation entrée-sortie u > y du système (B.1) soit identique ? Si oui, dessiner un graphe 
d’événements plus simple ayant le même comportement entrée-sortie. 


M en KL , (temporisations en toutes lettres 
S à F2 = dans les places) 
Roue 


FIG. B.2: Quelles sont les équations de ce graphe d'événements ? 


B.4 Calcul de transfert 


On considère une entreprise de vente par correspondance fonctionnant comme suit. En 
u arrivent des commandes par courrier. En x, un employé dépouille le courrier et vérifie les 
commandes, ce qui lui prend 10 minutes par commande. Il passe alors un bon à l’un ou l’autre 
de deux magasiniers (en x,) qui mettent chacun 30 minutes pour servir une commande. Une 
troisième personne (en x;) prépare les colis. Elle peut traiter deux colis en parallèle et peut 
commencer à travailler dès que l’employé en x; lui a transmis la commande, mais doit at- 
tendre que les magasiniers l’aient livrée pour fermer le colis. La préparation d’un colis lui 
prend au minimum 15 minutes (lorsque les magasiniers ne la font pas attendre). On vérifie 
que ce système se représente par le graphe d'événements de la Figure B.3, (a). 
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FIG. B.3: Service de livraison 


a/ Ecrire ce système sous la forme 
x = Ax@Bu, y—=Cx 


où À € (Rmaxl7 1) et B, C sont des matrices à coefficients constants. 


b/ Donner la relation entrée sortie sous la forme y — hu, où h € Raxlly 1]. On suppose 
qu’une infinité de commandes arrivent à l’instant O (et que le système ne fonctionnait pas 
auparavant). Dessiner la sortie y. Quel est le taux de production du système pour cette entrée. 
Que dire du nombre de bons en attente d’être servis par les magasiniers entre x; et x2. 


c/ Pour palier l’inconvénient observé à la fin du b/, on a fait appel à un cabinet de conseil 
en organisation qui à proposé la modification suivante. La personne qui expédie les colis (en 
x3) reçoit le bon des magasiniers et après avoir mis le colis à l’expédition, retourne le bon à 
X1, qui n’accepte maintenant de dépouiller une nouvelle demande que si le nombre de bons 
délivrés n’excède pas de plus de deux unités le nombre de bons retournés par x3. On obtient 
alors le nouveau réseau de Petri de la Figure B.3,(b). Calculer la série de transfert. Calculer la 
sortie y correspondant à la même entrée qu’en (b). Que peut on dire du taux de production ? 
Peut on améliorer la modification proposée par le conseil ? 
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Annexe C 


Travaux dirigés : automates et 
Systèmes à Evénements Discrets 


Modèle simplifié du Triangle de Gagny 


On considère la portion de réseau de chemin de fer de la figure C.1 (croisement en 
triangle). Il s’agit d’éviter la situation de blocage décrite à droite de la Figure. On considère 


FIG. C.1: Croisement en triangle 
l’alphabet 
Z = {ei,51,€2,52,e3,53} . 


On a l’ensemble des événements suivants : 
e, un train rentre sur le tronçon 1 
si un train sort du tronçon 1 
e) un train rentre sur le tronçon 2 


a/ À chaque tronçon i = 1,2, 3, on associe un état g; valant 0 où 1 : g; = 0 signifie aucun 
train sur le tronçon, g; — 1 signifie qu’un train est présent sur le tronçon. L'état global des 
trois tronçons est représenté par le triplet (g1, q2,g3) € Q où Q = {0,1}. Initialement, il 
n’y a aucun train sur les trois tronçons. Représenter le système par un Y automate À ayant 
pour ensemble d'états Q. 


b/ Le langage admissible K est l’ensemble des mots qui ne conduisent pas à la situation de 
blocage représentée à droite de la Figure C.1. 
Donner un automate raffinant À et reconnaissant le langage admissible K. 
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c/ On suppose que les entrées sur les tronçons e1, e2, e3 sont contrôlables et que les sorties ne 
le sont pas, soit la partition 


Ein = (51, 52,83), Ze = (e1, ©, 63} : 


Montrer que le langage K est contrôlable. Construire un superviseur. 


d/ On suppose en outre que e, est incontrôlable (la ligne 1 est une ligne rapide prioritaire que 
l’on ne souhaite pas interrompre). Le langage K est il contrôlable ? Calculer le sous langage 
contrôlable maximal K° (on donnera un automate reconnaissant KT). 


e/ En outre e n’est plus contrôlable. Calculer K°. 
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